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sites  were  studied  with  emphasis  on  the  time  dependence  of  the 
failure  process.  Three  subject  areasjeaxe^  identified: 

a.  Residual  stresses  due  to  changes  in  temperature  through 
the  glass  transition  range.  It  is  found  that  determina¬ 
tion  of  the  creep  compliance  or  of  the  relaxation  modulus 
is  the  most  important  material  property  for  accurate 
stress  determination  general  experimental  and  analytical 
agreement  prevails. 

b.  In  temperature  "accelerated"  crack  propagation  tests 
along  interfaces  it  is  found  that  crack  propagation  stops 
upon  raising  the  temperature.  The  reason  for  this  "unex¬ 
pected"  behavior  is  that  with  raising  the  temperature  the 
elimination  of  residual?  stresses  overcomes  the  reduction 
of  viscosity  so  that  crack  arrest  becomes  possible. 

c.  Fatigue  crack  propagation  is  studied  for  a  viscoelastic 
material.  In  contrast  to  metals  the  rate  of  crack  growth 
per  cycle  is  strongly  affected  by  the  frequency,  declin¬ 
ing  with  increasing  frequency.  However,  the  average  velo¬ 
city  (cm/sec)  per  cycle  increases  with  frequency.  Thus  it 
is  more  important  to  consider  the  time  under  stress  than 
merely  the  number  of  cycles.^ 
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1.  PURPOSE  AND  OBJECTIVE  OF  PROGRAM 


Polymers  are  increasing  in  prominence  in  the  design  of  aerospace  structures 
as  exemplified  by  the  use  of  epoxies  in  advanced  composite  materials  and  by 
adhesives  for  metal/metal  and  metal/composite  joining.  However,  only  in  very 
rare  instances  do  designers  worry  about  the  consequences  of  the  'time  depen¬ 
dent"  properties  which  these  materials  exhibit.  The  reasons  for  this  neglect  are 
manyfold;  however,  primary  amongst  them  is  the  general  lack  of  familiarity  with 
the  consequences  of  viscoelastic  material  behavior  on  the  one  hand,  and  the 
misjudgement  of  long-term  time  dependent  response  as  estimated  from  short 
term  laboratory  test,  on  the  other.  With  regard  to  the  latter  point  it  is  not 
uncommon  to  experience  laboratory  data  that  exhibit  time-  or  rate-sensitivity 
which  is  lost  in  the  data  scatter,  with  the  common  conclusion  that  time-  or  rate- 
sensitivity  is  absent  or  not  important.  Such  finding  is.  of  course,  justified  only  if 
the  application  of  the  polymer  is  likewise  limited  to  short  times  comparable  to 
those  encountered  in  the  laboratory  test,  but  is  not  justified  if  applications 
involve  time  spans  or  frequencies  several  orders  of  magnitude  different  than 
those  underlying  the  laboratory  test. 

From  a  structural  point  of  view  it  may  be  undesirable  to  design  structures  for 
extended  use  with  materials  that  possess  strongly  rate  sensitive  behavior. 
Nevertheless,  in  order  to  understand  the  limits  of  designs  incorporating  poly¬ 
mers  one  must  be  thoroughly  familiar  with  the  special  problems  derived  from 
viscoelastic  material  response. 

One  might  suggest  that  the  limits  of  viscoelastic  behavior  can  be  readily 
ascertained  by  simple  compliance  or  relaxation  testing.  For  some  deformation- 
limited  designs  that  may  be  true.  However,  our  recent  discovery  that  elevated 
stress  loads  accelerate  the  deformation  and  relaxation  processes  (non-linearly 


-  4  - 


viscoelastic  behavior)  points  out  that  standard  material  characterization  is 
marginally  adequate  at  best.  In  particular,  these  findings  point  out  that  because 
of  the  stress-induced  acceleration  of  creep  processes  the  fracture  behavior  of 
polymers  will  be  affected  significantly  Specifically,  the  high  stresses  at  the  tip  of 
a  crack  induce  time  dependent  phenomena  that  proceed  much  more  rapidly 
there  than  those  in  the  bulk  of  the  material  As  a  consequence  a  material  that 
might,  by  other  criteria,  appear  to  be  nearly  rigid  and  rate  insensitive  can 
indeed  exhibit  time  dependent  fracture  response.  Such  apparent  inconsistency 
is  observed  even  with  silicate  glasses.  It  is  thus  clearly  of  engineering  interest  to 
develop  and  disseminate  a  better  understanding  of  time  dependent  phenomena 
of  failure  in  structural  polymers. 

It  was  with  this  desire  and  purpose  in  mind  that  we  undertook  the  p.-esent 
program,  the  primary  intent  being  the  development  of  analysis  tools  for  the 
prediction  of  time  dependent  phenomena.  Because  the  physical  behavior  of 
these  materials  is  generally  poorly  understood  -  especially  with  regard  to  frac¬ 
ture  -  we  felt  that  an  experimental  program  with  some  analytic  undergirding 
was  an  appropriate  mode  of  operation. 

Specifically,  the  program  was  designed  to  elucidate  those  physical 
phenomena  that  control  the  time  dependence  of  the  failure  process.  Such 
failure  or  fracture  processes  cover  both  the  separation  of  a  viscoelastic  solid 
from  an  elastic  substrate  as  well  as  the  growth  of  cracks  in  rigid  polymers.  Both 
problems  are  germane  to  the  failure  of  bonded  joints  and  of  advanced  fiber  rein¬ 
forced  composites. 

2.  SUMMARY  OF  TECHNICAL  RESULTS 

We  feel  that  progress  has  been  made  on  all  fronts  with  one  exception.  We  del¬ 
ineate  the  progress  below,  but  comment  first  on  the  past  that  met  with 
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unforseen  difficulties. 

Part  of  an  earlier  program  on  which  the  present  effort  was  to  build  depended 
on  servomechamcal  control  by  an  8  bit  micro-processor.  It  turned  out  that  the 
Northstar  micro-processor,  while  only  about  four  years  old  at  the  time  a  prob¬ 
lem  arose,  was  claimed  to  be  too  old  to  be  worth  repairing.  The  processor  was 
finally  repaired  but  too  late  to  benefit  the  thesis  research  of  the  student  who 
was  to  use  it.  The  student  had  to  change  his  research  topic  and  has  since  com¬ 
pleted  his  thesis  (see  Section  2  3  below) 

We  have  learned  an  important  point  out  of  this  misfortune,  however.  Because 
we  were  totally  dependent  on  off-campus,  i.e.  outside  technical  service  and 
analysis  we  were  virtually  totally  subject  to  competition  for  repair  in  the  general 
market  place  but  with  a  relatively  trivial  problem  that  did  not  weigh  heavily  in 
the  general  competition.  Future  involvement  with  computer  will  entail  a)  more 
direct  involvement  in  understanding  and  designing  the  system,  plus  b)  working 
with  systems  for  which  expertise  exists  on  campus.  Thus  assistance  in  trouble 
shooting,  not  to  speak  of  cost  sharing,  will  be  of  direct  advantage  to  system 
maintenance. 

We  turn  now  to  outlining  briefly  the  three  major  areas  of  research  accom¬ 
plishments.  These  areas  are  delineated  in  two  theses,  one  a  Ph.D.  thesis  and  the 
other  a  thesis  for  the  degree  of  Aeronautical  Engineer.  The  details  are  well 
represented  in  these  documents,  which  are  appended  as  part  of  this  final  report. 
The  topics  contained  in  the  Ph  D  thesis  are 

a.  Analysis  of  residual  stresses  as  a  polymer  is  cooled  through  the  glass  trnsi- 
tion  temperture  and 

b.  Effect  of  residual  stresses  on  adhesive  failure. 


The  third  topic, 
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c.  relates  to  the  effect  of  viscoelastic  material  behavior  on  fatigue  crack  propa¬ 
gation.  The  latter  topic  is  the  subject  of  the  Engineer's  Thesis. 

Topics  a)  and  b)  are  subjects  that  were  motivated  by  bonding  problems  rela¬ 
tive  to  solar  cell  construction  but  are  of  a  much  more  fundamental  mature  than 
energy  engineering.  Although  that  reserch  was  started  under  DOE  sponsorship 
the  problems  considered  are  precisely  those  of  structural  adhesive  use,  in  fact, 
are  especially  germane  to  fiber  reinforced  composite  materials.  That  work  was 
substantially  done  and  completed  under  the  present  grant. 

2.1  Residual  Stresses  in  Polymers  Cooled  Through  the  Glass  TranantionTempcrature 

In  may  aerospace  structures  polymers  are  formed  together  with  non¬ 
polymers  (metal,  graphite  fibers)  while  the  polymer  is  in  the  semi-liquid  or  ’^oft" 
state  As  the  temperature  drops  material  shrinkeage  occurs  which  undergoes  a 
transition  at  the  glass  temperature.  While  this  is  a  qualitatively  well-known  fact 
the  quantitative  aspects  were  not  at  all  explored.  Because  of  the  great  sensitivity 
of  the  mechanical  properties  to  the  temperature  around  the  glass  transition 
this  was  always  judged  to  be  a  very  difficult  problem.  Nevertheless,  it  is  a  very 
important  engineering  problem  because  ail  polymer  using  structures  involve 
that  question.  Specifically,  if  it  turns  out  that  these  residual  stresses  are  large 
enough  to  almost  cause  fracture  without  much  additional  loading  it  would  be 
important  to  assess  the  residual  strength. 

The  study  revolved  thus  around  the  following  question:  Given  the  current 
knowledge  and  state  of  affairs  in  analyzing  thermoviscoelastic  behavior,  and 
given  material  characterization  techniques  that  considerably  exceed  in  preci¬ 
sion  what  is  practice  in  industry  and  moBt  engineering  schools,  how  well  can  one 
predict  the  residual  stresses  and/or  deformations  in  composite  structures  as 
the  temperature  passes  through  the  glass  transition  range.  One  might  call  this  a 
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careful  sensitivity  study. 

The  answer  to  this  question  is  represented  essentially  graphically  in  Figures 
24  through  27  in  Appendix  A  (starting  at  page  56).  These  figures  represent  com¬ 
puted  and  interferometrically  measured  radii  of  curvature  of  sandwiches  made 
of  glass  and  Polyvinylacetate.  It  turns  out,  that  the  most  influential  parameter 
in  this  comparison  is  the  creep  or  relaxation  behavior  (see  Figure  5  in  Appendix 
A)  Although  on  a  log-log  plot  such  data  may  'look  very  good"  it  turns  out  that 
most  of  the  deviation  in  the  Figures  24-27  in  Appendix  A  between  theory  and 
experiment  results  from  the  uncertainty  (variability)  of  the  time  dependence  of 
the  material  rheology.  The  data  in  Figure  5  was  measured  painstakingly  -  by  nor¬ 
mal  engineering  standards  -  and  repeatedly  on  high  purity  PVAc. 

The  net  result  of  this  study  is  thus  that  if  predictions  of  residual  stresses 
need  to  be  made  to  within  even  only  10-15  %  accuracy  such  an  estimation  will  be 
practically  speaking  impossible  within  the  normal  run  of  engineering  practice 
and  applications. 

2.2  Effect  of  Residual  Stresses  on  Adhesive  Failure 

One  of  the  most  important  subjects  in  the  assessment  of  durability  of  struc¬ 
tures  is  the  ability  to  either  predict  long  term  (10-30  years)  failure  or  to  possess 
accelerated  test  methods  with  the  help  of  which  one  can  determine  experimen¬ 
tally  whether  a  specific  design  will  endure.  The  most  common  means  of 
accelerating  tests  involving  polymers  is  to  test  under  elevated  temperatures. 
When  more  than  one  material  is  involved  such  temperature  changes  induce 
thermal  stresses  which  can,  by  themselves,  cause  fracture.  Again,  one  of  the 
most  often  used  tests  to  assess  adhesion  is  the  peel  test;  it  seemed  appropriate, 
therefore,  to  examine  the  effect  of  temperature  on  peel  rate. 


At  the  same  time  one  notes,  that  in  the  peel  test  there  is  rarely,  if  at  all,  men¬ 
tion  made  of  the  thickness  of  the  peeled  layer  other  than  a  mere  recording,  but 
not  as  a  recognition  that  peel  layer  thickness  would  or  could  affect  "peel 
strength." 

While  specimen-to-specimen  variation  renders  fairly  large  data  scatter  the 
answer  to  the  latter  question  is  that  there  is  a  dependence  of  peel  strength  in 
terms  of  energy  expended  per  unit  of  new  interface  surface.  The  experimental 
results,  while  not  the  best  in  data  clarity,  indicate  a  non-monotonic  dependence, 
which  makes  sense,  if  one  considers  that  the  viscoelastic  energy  expenditure 
depends  on  the  detail  of  the  strain  field  in  the  vicinity  of  the  disbond.  For  a 
thick  specimen  the  strain  rates  are,  on  the  average  low  because  the  local  bend 
radius  is  large  leading  to  "smaller"  strains  and  strain  rates,  hence  less  dissipa¬ 
tion.  At  the  same  time,  physical  observation  shows  that  more  volume  is  coming 
under  strain  the  larger  the  thickness  is,  not  only  because  of  the  increased  thick¬ 
ness  but  because  there  is  more  material  along  the  length  of  the  peel  layer  under 
strain  because  the  mean  bend  radius  is  larger  Moreover,  the  strain  field  tends 
to  become  more  non-uniform  as  the  thickness  increases  (see  page  87  of  Appen¬ 
dix  A).  The  net  effect  allows  then  for  a  maximum  of  viscoelastic  dissipation  at 
some  thickness  of  the  peel  layer. 

In  the  second  part  of  the  peel  study  it  was  examined  what  the  effect  of  resi¬ 
dual  thermal  stresses  was  on  the  peel  rate  at  different  temperatures.  It  turns 
out  that  in  a  certain  temperature  range  peel  does  not  take  place  while  under 
the  same  load.  The  explanation  is  found  in  the  fact  that  below  a  certain  tem¬ 
perature  the  thermally  induced  strain  energy  assists  the  fracture  process  When 
the  temperature  rises  to  a  given  value  the  thermal  and  cure  shrinkeage  induced 
strain  energy  no  longer  assists  fracture  and  the  crack  propagation  process 
stops  This  behavior  is  opposite  to  what  one  expects  normally  on  purely 
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rheological  grounds.  Crack  growth  can  then  be  promoted  only  by  increasing  the 
peel  force. 

2.3  Recent  Results  on  the  Validity  of  a  Non-Lin  early  ViscoelasticRepresentatian 

One  of  the  recurring  topics  in  the  failure  or  fracture  of  viscoelastic  solids  is 
the  question  of  non-linearly  viscoelastic  material  behavior.  In  particular  it  is  the 
goal  of  materials  behavior  research  to  find  that  mathematical  representation  of 
experimentally  accessible  states  of  stress  or  strain  histones  which  lend  them¬ 
selves  to  generalizations  to  less  readily  tested  states. 

Leaderman  has  suggested  many  years  ago  that  a  generalization  of  linearly 
viscoelastic  behavior  be  achieved  through  a  convolution  of  a  material  function, 
say  the  relaxation  modulus  with  a  non-linear  function  of  strain  instead  of  the 
strain  itself.  On  the  other  hand,  an  equally  acceptable  representation  in  terms 
of  the  creep  compliance  would  be  the  convolution  of  the  creep  compliance  with 
a  non-linear  function  of  the  (Cauchy)  stress.  The  latter  form  is  the  one  used  in 
Schapery's  "generalization"  of  the  J-integral  to  include  viscoelastic  material 
behavior. 

The  particular  form  of  material  representation  with  a  special  form  for  the 
non-linear  function  of  the  strain  was  used  by  Bloch,  Chang  and  Tschoegl  [i]  to 
generate  excellent  fit  of  experimental  data.  Let  e  be  the  normal  engineering 
strain  and  let  the  uniaxial  extension  ratio  be  X  =  1  +  e.  Then,  if  n  is  some  con¬ 
stant  at  a  temperature  T„  the  stress  can  be  represented  by 

a(t)  =  ^-/E( t  -u)  ^-}Xn(u)  -X  2(u)j  du  (1) 

The  corresponding  fit  to  experimental  data  is  shown  in  Figure  1,  details  as  to 
how  the  single  free  parameter  n  is  fitted  by  the  data  is  recorded  in  the  cited 
reference  Note  that  the  data  in  Figure  1  is  obtained  from  monotonic  loading 
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histories,  in  particular  from  uniaxial  constant  strain  rate  tests. 

In  many  situation  it  is  of  interest  to  deal  with  periodic  deformation  histones 
such  as  are  encountered  in  fatigue-type  problems  In  such  situations  it  is  of 
paramount  importance  to  be  able  to  estimate  the  amount  of  energy  dissipation 
in  a  material  in  order  to  assess  the  associated  temperature  rise  It  was  the 
intent,  therefore,  to  estimate  with  the  excellent  material  representation  in  mind 
the  amount  of  energy  dissipation  per  cycle  for  this  material  and  to  compare 
that  estimation  with  that  for  the  small  strain  behavior,  i.e.  for  linearly  viscoelas¬ 
tic  behavior 

In  order  to  compute  the  energy  dissipation  per  cycle  one  calculates,  for 
umaxial  deformation  the  loss  per  cycle 

T  T 

Wc  =  /a(t)  ^-dt  =  /(j(t)edt.  (2) 

O  o 

Suppose  we  consider  only  a  cyclic  tensile  strain  be  letting 

£  =  s0  (1  +  sincjt)  (3) 

_IL 

Let  the  function  in  curly  brackets  be  F*(u)  =  \n(u)  -  X  2  (u)  so  tha^  with  (2) 

_  n_ 

F*(u)  =  [  1  +  e0(l  +  sinuu)  ]n  -  [  1  +  e„(l  +  sinuu)  ]  2  (4) 

which  has  the  Fourier  series  representation 

bn 

V  - sinnwu - cosna>n-b0 

^  nu  ncj 

F*(u)  =  =  — 1 —  Va„  +•  b2  sin(ncju  +  dj  (5) 

ncj 

tan<5„  = - — 

Substitution  of  (5)  into  (1)  yileds  for  steady  state  condition  and  by  definition  of 
the  complex  modulus  E*(i>)  =  E(u)el4(omega'  and  the  Loss  tangent  tan  A  (omega), 
the  periodic  stress  is 
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PREFACE 

PoiVrner.c  coatings  are  often  applied  in  order  to  protect  a  substrate  from  the 
.nf. -er.ee?  of  tne  environment  If  polymer  films  are  to  adhere  for  many  years  to 
-'■-i'  r.  ■  -os’  rates  .t  s  necessary  to  understand  how  material  and  system  parame- 
•■‘*e :•  :r.e  process  of  separation  between  coating  and  substrate.  The 
one  no  me-  i  r.vo.vec  piay  a  recurring  role  in  all  problems  of  separation  at  an 
r.  ter'  ice  net  ween  dissimilar  materials. 

The  present  work  deals  with  the  long  time  integrity  of  polymer  coatings.  The 
goal  of  this  study  is  the  long  time  prediction  of  the  initiation  and  propagation  of 
an  interfacial  crack,  that  can  lead  to  the  debonding  of  a  protective  layer  The 
immediate  motivation  for  this  study  arose  in  the  need  for  encapsulation  of  solar 
photovoltaic  cells.  These  have  become  important  as  a  possible  alternative 
source  of  energy  However,  the  work  presented  here  should  be  viewed  m  a  more 
general  manner  because  it  applies  to  other  coatings  subject  to  similar  environ¬ 
mental  loading  conditions.  For  example,  the  peeling  of  paint  from  its  substrate 
falls  into  this  category  Due  to  a  changing  atmosphere  or  possibly  deficient 
manufacturing  procedures  the  coating  materials  might  no  longer  adhere  per¬ 
fectly.  opening  the  path  to  further  degradation  of  the  adhesive  bond  Problems 
of  this  kind  occur  in  the  encapsulation  of  solar  cells  and  are  exemplified  in  Fig¬ 
ure  1  which  illustrates  a  delammation  in  a  solar  ceil  array  ;  the  coating  has 
separated  from  its  underlying  substrate  after  lengthy  exposure  to  environmen¬ 
tal  loading  conditions.  The  appearance  of  such  structural  defects  indicates  the 
need  to  determine  under  what  conditions  the  adhesive  bond  will  continue  to  fail. 
The  present  study  deals  with  two  distinct  aspects  related  to  the  long  time 
integrity  of  a  viscoelastic  matenal  layer  bonded  to  an  elastic,  rigid  substrate. 


The  first  part  concerns  the  magnitude  of  the  stresses  present  in  layered 
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erties  of  the  adhesives  needs  to  be  carefully  analysed.  This  time  dependence  is 
also  reflected  in  the  energy  required  to  create  new  surfaces  as  interfacial 
debonding  proceeds  ;  the  adhesive  fracture  energy  is  one  of  the  dominant 
parameters  in  time  dependent  adhesive  failure.  In  our  investigation  it  is  charac¬ 
terized  through  peel  testing. 

With  the  knowledge  of  the  pertinent  material  properties  as  well  as  of  the 
adhesive  fracture  energy,  we  then  proceed  to  formulate  a  criterion  for  continu¬ 
ing  interfacial  crack  propagation.  The  analysis  is  carried  out  for  elastic  solids, 
with  the  effect  of  viscoelastic  behavior  incorporated  later  on.  Debond  tests  pro¬ 
vide  a  way  to  check  how  well  the  theoretical  predictions  correspond  to  experi¬ 
mental  debond  results. 
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ABSTRACT 

An  important  source  of  interface  fracture  contributing  to  adhesive  failure  in 
a  bunaterial  sandwich,  consisting  of  a  rigid  substrate  and  a  viscoelastic  encap- 
sulant  material,  arises  from  residual  stresses.  The  encapsulant  is  often  depo¬ 
sited  on  the  substrate  above  its  glass  transition  temperature  region  but  used 
below  this  temperature  range.  In  order  to  determine  the  magnitude  of  the  resi¬ 
dual  stresses  a  viscoelastic  stress  analysis  of  a  bimaterial  sandwich  is  carried 
out,  taking  into  account  the  time-dependent  material  properties  of  the  polym¬ 
eric  layer  and  the  environmental  "loading"  conditions.  The  theoretical  analysis  is 
paralleled  by  an  experimental  examination  of  the  time-dependent  out-of-plane 
deformation  of  thin,  circular  sandwiches. 

Polyvinyl  acetate  was  chosen  as  a  model  material  exhibiting  significant 
viscoelastic  effects  under  room  test  conditions.  Therefore  the  pertinent  physical 
and  mechanical  properties  of  PVA<.  are  determined  ;  these  include  the  thermal 
coefficient  of  expansion,  the  shear  creep  compliance  and  the  relaxation 
modulus.  In  the  experimental  work  BK-7  glass  is  taken  as  the  "rigid"  substrate. 
The  measurements  connected  to  the  stress  analysis  are  monitored  with  laser 
interferometry  (Newton’s  rings).  A  comparison  between  theory  and  experiment 
completes  the  viscoelastic  stress  analysis. 

In  the  second  part  of  this  study  time  dependent  adhesive  failure  of  rubbery 
materials  is  investigated.  Polymeric  materials  are  being  used  increasingly  for  a 
wide  variety  of  applications.  Some  of  these  materials  are  applied  as  protective 
layers  to  isolate  their  substrates  from  a  hostile  environment.  Others  achieve 
remarkable  structural  bond  strengths  thereby  displacing  the  traditional 
mechanical  fasteners  like  bolts  and  rivets.  If  one  wants  to  investigate  the  long 
time  integrity  of  a  layer  assembly  the  time  dependence  of  the  material  prop- 
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2.4  Effect  of  Viscoelastic  Behavior  on  Fatigue 

It  is  a  commonly  acceptedfact  that  crack  growth  rates  (da/dn)  under  fatigue 
type  loading  does  not  depend  on  test  frequency  (in  metals).  As  a  result  one  can 
accumulate  fatigue  data  quickly  at  high  frequencies.  In  viscoelastic  materials, 
we  have  calculated1  that  under  cyclic  loading  cracks  should  propagate  more 
rapidly  on  a  per  unit  time  basis  the  higher  the  frequency  That  prediction  is 
verified  in  the  following  sense:  Crack  growth  per  cycle  decreases  as  the  fre¬ 
quency  increases.  Yet  the  decrease  is  such  that  when  one  considers  the  average 
crack  speed  (per  cycle )  that  latter  quantity  increases  with  frequency.  Thus  it  is 
found  that  it  is  essentially  the  time  under  load  which  controls  crack  growth  and 
not  the  deformations  per  unit  cycle  as  in  metals.  These  results  are  summarized 
in  Figures  4.8  and  4. 12  in  Appendix  B. 
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Ph  D.  and  one  Professional  Engineer's  Degree.  Dr.  Luc  Heymans  now  works  for 
the  Eastman  Kodak  Compay  in  Rochester,  N.Y.,  and  Engr.  Samuel  Chang  is 
employed  by  American  Telecom  in  Anaheim,  CA.  Furthermore  the  following  stu¬ 
dents  have  contributed  to  the  program  and  have  derived  educational  benefits 
from  it:  Mile  Moser  and  Sotirios  Natsiavas. 

4.  COUPLING  ACTIVITIES 

Under  this  grant  the  above  mentioned  degree  students  produced  thesis  which 
are  included  here  as  Appendices  A  and  B.  These  theses  will  be  written  into  three 
publications  according  to  the  outlines  in  Section  2. 1-2.3. 

1.  W.G.  Knauss  and  H.  Dietmann,  'Crack  Propagation  Under  Variable  Load  Histones  in  Linearly 

Viscoelastic  Solids,"  fni.  J.  Ehgng  Set,  Vol  8,  pp.  843-658  (1870). 
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the  linearly  viscoelastic  solid  which  has  the  same  relaxation  modulus.  It  is  of 
interest  to  note  that  for  the  material  represented  in  the  measurements  in 
Bloch,  Chang  and  Tschoegl's  work,  i.e.  when  n*  1.2  the  energy  dissipated  by  the 
non-linear  body  is  about  90  %  of  that  for  the  linearly  viscoelastic  one,  and  in  a 
way,  that  is  rather  insensitive  to  the  level  of  the  maximum  strain  (in  the  range 
examined) 

This  result  would  seem  to  be  of  considerable  interest  inasmuch  as  it  would 
represent  the  first  estimate  of  the  effect  of  non-linearly  viscoelastic  material 
behavior  on  hysteresis.  There  is,  however,  one  serious  catch  to  this  result  which 
questions  the  validity  of  the  representation  in  equation  (l)  and  thus  questions 
also  the  validity  of  Schapery's  non-linear  J-integral  extension.  The  latter  state¬ 
ment  derives  from  the  fact  that  Schapery  uses  the  same  type  of  representation. 


Noting  that  integrals  of  the  type 


Zn  u 

f  sinnut  sinmwt  dt 


0  n  /  m 
rr/2n  =  m 


are  involved  in  evaluating  (2)  we  recognize  that  in  the  evaluation  of  (2)  with  the 
use  of  (7)  only  the  fundamental  frequency  in  the  stress  expression  contributes 
to  the  dissipation  integral.  Specifically,  all  the  terms  in  the  series  representa¬ 
tion  for  the  stress  higher  than  the  fundamental  mode  do  not  contribute  to  dissi¬ 
pation.  That  result,  being  entirely  a  consequence  of  the  mathematical  represen¬ 
tation  of  the  stress,  is  highly  suspect  in  its  physical  correctness.  In  fact,  each 
term  in  the  series  representation  contains  a  term  that  involves  a  phase  shift 
due  to  the  loss  tangent.  Yet,  those  terms  do  not  contribute  to  the  loss.  We  can 
only  conclude  that,  inspite  of  the  excellent  representation  for  montomc  load 
histories,  that  even  this  non-linear  formulation  of  only  the  unieaxial  stress- 
strain  behavior  is  not  universally  valid. 
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h 


a 


o  ■»  ^/a  ^  x  J»2 

ff(t)  =  - —  £  E(mcj)  — 1:2 - — sln[  mwt  +  <5m  +  A(mcj)  ]  (5) 

It  turns  out  for  all  the  parameters  considered  below  b0  is  very  small  and  may 
be  taken  to  be  zero,  and  b!  is  identically  zero.  We  have  thus,  more  explicitly, 

a  =  j^-E(u)  ^-sin[  wt  +  A(u)  ] 

o  99  '\/a  ^  4-  h  ^ 

+  3^-  £  E(mu)  — ”  -  —  sin[  mat  +  Am  +  A  {mu)  ]  (7) 

Substitution  of  (7)  into  (2)  yields  then  the  energy  dissipated  per  cycle  for  the 
non-linearly  viscoelastic  material.  This  computation  has  been  carried  out  for  a 
range  of  parameters  "n"  as  well  as  for  a  range  of  extension  rations. 

One  finds  for  the  energy  dissipated  per  cycle  Wc 

Wc  =  ^p-(— )e0E"(cj)  (E"  =  imaginary  part  of  the  dynamic  modulus) 

•3  CJ 

One  can  show  that 

n+2 

ara  =  Kl  +  e0)n_1  +  £<1  +  E0)  2  i  £0Amcj 

with  A„  a  sequence  of  numbers,  the  first  of  which,  At,  does  not  vary  much  from 
unity  in  magnitude.  In  considering  the  limit  e„  -*  0  one  finds  for  the  dissipation 
at  infinitesimal  strain  per  cycle  Wum 

WUm  =  rre|E"(«)  . 

So  that  the  ratio  for  ’hon-linear"  to  ‘linear"  histories  is 

R  s  ^-=  §*<1  +  £o)n'*  +  J<1  +  O'^iA, 


The  results  are  presented  in  Figure  2.  Here  I  have  defined  as  "R"  the  ratio  of 
the  energy  dissipated  by  the  non-linearly  viscoelastic  solid  to  that  dissipated  by 


Figure  1.  Interfacial  failure  between  solar  cells. 
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structures  involving  viscoelastic  encapsuiant  materials.  Some  polymers  or 
encapsulant  materials  are  likely  to  be  deposited  on  their  substrates  at  higher 
temperatures  followed  by  a  cooldown  to  the  working  temperature  range.  The 
first  part  of  this  work  focuses  therefore  on  an  assessment  of  how  well  we  are 
presently  able  to  predict  such  stresses.  At  this  stage  interfacial  fracture  is  not 
considered.  For  purposes  of  predicting  long  time  stability  of  an  adhesive  bond  it 
may  become  necessary  to  compute  and  determine  the  magnitude  of  these  resi¬ 
dual  stresses.  In  particular,  if  the  temperature  at  which  the  coating  is  bonded 
to  its  substrate  lies  above  the  glass  transition  temperature  of  the  viscoelastic 
material  a  cooling  through  the  glass  transition  range  may  not  be  avoidable.  This 
cooling  is  accompanied  by  relaxational  behavior  as  well  as  by  changes  in  the 
thermal  expansion  characteristics  of  polymeric  materials.  The  resulting  ther¬ 
mal  stresses  or  strains  induced  by  the  mismatch  in  thermal  expansion  between 
adherend  and  substrate  are  then  important.  For  most  encapsulant  materials 
adhering  to  a  rigid  sublayer  the  difference  between  thermal  expansion 
coefficients  typically  varies  around  5  x  10-4/°C  ;  a  cool  down  of  100°C  (180°F) 
may  readily  introduce  several  percent  strain  and  cannot  thus  be  disregarded 
unequivocally.  With  very  rigid  substrates  the  thermally  stored  energy  will  not  be 
disposed  of  through  bending.  Hence  more  driving  force  remains  available  for 
interfacial  fracture. 

2.  The  second  type  of  problem  under  study  here  concerns  the  separation  of  a 
rubbery  material  from  its  substrate,  using  a  fracture  mechanics  approach.  The 
behavior  of  rubbery  materials  differs  markedly  from  the  one  exhibited  by  the 
materials  studied  in  the  first  part  of  this  work  ;  there  the  cooldown  from  higher 
temperatures  through  the  glass  transition  region  makes  the  polymeric  materi¬ 
als  behave  increasingly  like  a  hard  or  even  brittle  material.  The  investigation  of 
debonding  of  such  rigid  or  semi-rigid  polymers  is  very  difficult  ,  as  a  precursor 
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to  such  studies  it  seems  prudent  to  deed  first  with  debonding  of  rubbery  materi¬ 
als  which  allow  more  reasonable  assumptions  in  the  stress  analysis  of  the 
failure  or  crack  propagation  process.  In  the  study  of  the  debonding  of  rubbery 
materials  the  viscoelastic  character  of  the  coating  or  encapsulating  material 
can  be  accounted  for  more  readily.  Especially  the  high  stresses  and  strains 
around  an  interfacial  crack  tip  generate  the  need  to  incorporate  the  relevant 
viscoelastic  phenomena  in  any  fracture  analysis  for  adhesive  bonds.  The  former 
are  difficult  to  determine  accurately  and  one  may  have  to  be  satisfied  with  rea¬ 
sonable  estimates  for  their  amplitude. 

Besides  material  properties  and  the  evaluation  of  stresses  and  strains,  other 
parameters  involved  in  time  dependent  adhesive  failure  need  investigation.  One 
of  these  parameters  arises  from  the  need  to  quantify  the  surface  energy  associ¬ 
ated  with  the  creation  of  new  surfaces.  For  crack  growth  in  viscoelastic  elasto¬ 
mers  it  is  possible  to  define  fracture  progression  in  terms  of  a  rate  dependent 
fracture  energy.  Such  a  formulation  is  possible  as  long  as  there  is  no  significant 
creep  or  relaxation  in  the  body  far  from  the  crack  tip  so  that  an  elastically 
stored  energy  can  be  defined,  and  viscous  dissipation  due  to  crack  propagation 
is  confined  to  a  small  region  in  the  crack  tip  vicinity.1  For  rubbery  adhesives 
there  is  an  analogous  function  for  interface  separation  that  needs  to  be  deter¬ 
mined  in  order  to  assess  the  failure  rate  of  bond  interfaces.  Inasmuch  as  an 
elastic  stress  analysis  is  conducted,  the  time  dependent  failure  process  can  be 
described  solely  through  this  material  function  of  the  adhesive  and  adherend. 
In  the  kind  of  vanishingly  small  interfacial  crack  growth  rates  the  problem 
reduces  to  one  of  purely  elastic  fracture  mechanics.  In  the  second  portion  of 
this  thesis,  which  deals  with  the  fracture  analysis,  we  are  thus  concerned  with  a 

quasi-elastic  energy  consideration  in  which  the  forces  driving  the  interfacial 

1.  H.X.  Muller  and  W  G.  rCnauss,  The  Fracture  Energy  and  Some  Mechanical  Properties  of  a 

Polyurethane  Elastomer,  Trans.  Soc.  Rheoi.,  15  (1971)  pp.  217-233. 
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fracture  process  are  derived  from  temperature  changes  as  well  as  residual 
stresses  induced  by  cure  shrinkage  of  the  elastomer.  Additional  steady  state 
stresses  can  be  induced  through  volume  changes  caused  by  the  outgassing  of 
low  molecular  weight  components  of  the  polymer  or  by  ingression  of  water  or 
organic  solvents. 

The  sources  of  forces  driving  interface  fracture  derive  from  a  variety  of 
environmental  conditions,  mainly  diurnal  and  seasonal  changes  to  which  protec¬ 
tive  polymer  Aims  may  be  exposed.  The  cyclic  loading  pattern  of  these  forces 
gives  rise  to  a  fatigue  type  load  history  for  the  adhesive  bond.  Small  values  of 
the  amplitude  of  the  temperature  differential,  whether  between  day  and  night  or 
between  seasonal  averages  can  contribute  to  sizeable  strains.  In  an  analogous 
manner  possible  material  sensitivity  to  humidity  changes  the  integrity  of  the 
bond  interface.  This  combined  effect  of  temperature  and  moisture,  in  addition 
to  the  cumulative  fatigue  damage  over  longer  periods  of  time,  significantly 
affects  the  lifespan  of  an  encapsulated  solid. 

Interfacial  failure  is  not  a  new  problem.  It  has  been  investigated  from  a 
variety  of  viewpoints  most  of  which  are  connected  with  surface  chemistry  or 
surface  physics.  In  that  connection  it  is  useful  for  later  reference  to  recall  here 
the  relative  size  scales  associated  with  the  technical  problem  that  motivates  this 
work.  They  include  the  in-plane  dimensions  and  sizes  of  the  various  layers  in  the 
composite  structure.  Where  the  panel  containing  the  solar  cells  is  on  the  order 
of  meters  (or  yards),  the  thicknesses  of  the  protective  coatings  are  several  ord¬ 
ers  of  magnitude  smaller.  For  later  reference  we  note  that  : 

1.  A  silicon  cell  typically  measures  1  to  2  mm  In  thickness. 

2.  Protective  layers  are  thinner,  namely  on  the  order  of  10-1  to  10-4  mm. 


3.  The  interphase  or  transition  between  dissimilar  solids  is  about  10~°  mm. 

It  is  also  appropriate  to  record  here  the  time  scales  of  interest.  If  an  interfa¬ 
cial  crack  propagates  at  a  rate  as  low  as  1  x  10-7  centimeter  per  minute,  such  a 
debond  rate  would  produce  a  separation  on  the  order  of  1  centimeter  after  20 
years,  which  is  a  typical  projected  lifespan  for  an  encapsulated  solar  cell.  On 
the  other  end  of  the  time  spectrum  one  has  to  consider  the  material  relaxation 
time  of  the  encapsulant.  In  order  to  exclude  viscoelastic  phenomena,  the  glass 
transition  temperature  range  has  to  remain  sufficiently  far  below  the  operating 
temperatures  if  one  wants  to  deal  with  rubbery  materials.  At  the  crack  tip  how¬ 
ever  material  relaxation  and  viscous  dissipation  demand  a  more  careful 
analysis.  From  these  considerations  the  importance  of  the  size  and  time  scales 
appears  more  clearly. 

Summarizing,  it  is  useful  to  keep  in  mind  that  the  two  problems  under  study 
in  this  engineering  analysis  of  polymer  film  adhesion  are  necessarily  interre¬ 
lated  through  the  common  viscoelastic  phenomena  occurring  in  both  cases. 
Nevertheless  the  present  investigation  into  the  debonding  of  rubbery  materials 
is  limited  to  the  formulation  of  a  crack  propagation  criterion  using  an  energy 
approach. 
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1.  INTRODUCTION 

We  consider  now  the  9tress  analysis  for  a  viscoelastic  layer  bonded  to  a  thin 
substrate  and  cooled  through  its  glass  transition  temperature  range.  It  is  our 
intention  to  quantify  the  residual  stresses  in  the  adhesive  bond  through 
theoretical  calculations  as  well  as  experimental  measurements.  The  validity  of 
the  assumptions  made  in  the  theoretical  analysis  can  thus  be  checked  against 
the  experimental  results.  At  the  same  time,  the  determination  of  the  accuracy 
of  the  experimental  techniques,  in  particular  of  the  physical  property  determi¬ 
nation  is  a  corollary  purpose. 

It  is  clear  that  this  stress  analysis  requires  a  ''careful”  material  characteriza¬ 
tion.  The  need  to  include  the  viscoelastic  character  of  the  polymer  points  the 
choice  towards  materials  exhibiting  their  glass  transition  within  easy  experi¬ 
mental  reach,  i.e.  close  to  ambient  temperatures.  For  this  reason  polyvinyl  ace¬ 
tate  (PVAc),  an  uncrosslinked  amorphous  polymer  and  of  narrow  molecular  dis¬ 
tribution  is  chosen  as  the  model  material,  since  it  possesses  a  glass  transition 
temperature  zone  around  30  °C. 

Since  the  bulk  modulus  does  not  change  by  far  as  much  as  the  shear 
modulus  throughout  the  transition  region,  it  is  proposed  for  engineering  pur¬ 
poses  that  a  constant  value  is  a  reasonable  approximation.  This  assumption  is  a 
part  of  the  present  check  on  our  ability  to  predict  thermoviscoelastic  stresses  in 
the  bimaterial  sandwich.  The  same  assumption  however  does  not  hold  with 
respect  to  the  relaxation  modulus  or  creep  compliance  in  shear.  These  material 
properties  need  to  be  determined  accurately  over  a  wide  time  interval,  Even 
though  the  shear  relaxation  modulus  is  used  later  on  in  the  stress  analysis,  the 
creep  compliance  is  measured  because  a  torsion  apparatus  measuring  the  latter 
was  available  The  relaxation  modulus  could  thus  be  calculated  from  the  creep 
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data. 

The  temperature  dependence  of  the  relaxation  mechanisms  in  amorphous 
polymers  is  well  known  [l].  Because  of  the  great  sensitivity  of  the  material 
behavior  to  temperature  variations  the  temperature  control  during  the  meas¬ 
urement  of  the  viscoelastic  material  properties  must  be  very  fine  in  order  to 
avoid  disproportionately  sizeable  changes  in  the  time  scale  of  the  material 
response.  In  a  similar  manner  the  sensitivity  of  the  polymer  to  moisture  [2] 
makes  it  necessary  to  keep  humidity  levels  under  control  so  as  not  to  introduce 
any  stray  effects  into  the  measurements. 

Because  stresses  are  introduced  in  the  bimaterial  sandwich  by  thermally 
induced  volume  changes  the  pertinent  expansion  characteristics  need  to  be 
determined,  in  a  range  including  the  glass  transition.  Volume  changes  result 
from  thermal  loading  histories  or  curing  of  the  polymeric  material.  In  order  to 
avoid  grappling  with  thermal  transients  the  thickness  of  the  specimens  was 
chosen  to  be  so  small  so  as  to  ensure  a  uniform  temperature  throughout  the 
sample  at  all  times. 

Next  we  review  several  analytical  treatments  of  problems  related  to  the 
current  interest.  Several  decades  ago  work  on  bimaterial  laminates  concen¬ 
trated  on  bimetallic  strips.  Timoshenko  [3]  reported  on  his  investigation  in  con¬ 
nection  with  bimetal  strip  thermostats  usi^g  a  general  theory  for  bending 
deformation  of  a  bimetal  strip,  submitted  to  uniform  heating. 

An  early  analysis  of  problems  related  to  failure  in  bonded  joints  was  carried 
out  by  Goland  and  Reissner  [4].  In  their  investigation  of  cemented  lap  joints, 
they  obtained  explicit  solutions  for  the  stresses  for  two  limiting  cases,  i.e the 
case  of  a  very  thin  cement  layer,  and  the  case  where  the  joint  flexibility  is 
mainly  due  to  that  of  the  cement  layer,  More  recently  the  growing  interest  in 
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adhesive  bonding  for  automotive  and  aerospace  applications  has  triggered  an 
extensive  research  effort.  The  advent  of  the  finite  element  method  opened  up 
new  ways  to  analyze  adhesive  joints.  It  is  however  beyond  the  scope  of  this  study 
to  include  a  broad  review  of  the  work  related  to  the  knowledge  of  the  stress  dis¬ 
tribution  in  different  kinds  of  adhesive  joints  For  this  the  reader  is  referred  to  a 
paper  by  Kinloch  [5].  That  author  also  compiled  an  elaborate  review  of  surface 
and  interfacial  phenomena  relevant  to  the  adhesive  joining  of  materials  [6].  In 
the  following  we  limit  ourselves  to  review  research  efforts  more  closely  related  to 
the  topic  of  this  investigation. 

Durelli  and  coworkers  [7,8,9]  conducted  an  extensive  study  of  stresses  and 
strains  in  rectangular  strips  of  epoxy  material,  bonded  on  one  face  to  a  rigid 
plate  and  then  shrunk.  Conducted  by  means  of  photoelasticity  techniques,  they 
included  slabs  with  different  edge  geometries  in  their  investigation  In  their  work 
on  interfacial  crack  propagation  Mulville  and  Vaishnav  [10]  determined  via  pho¬ 
toelastic  techniques  the  magnitude  of  the  residual  stresses  due  to  casting  and 
curing  of  an  epoxy  bonded  to  aluminum.  These  stresses  were  found  to  contri¬ 
bute  up  to  20  percent  of  the  strain  required  for  crack  initiation  at  the  interface 
Weitsman  [ll]  calculated  approximately  residual  thermal  stresses  that  develop 
within  the  resin  of  a  fiber  reinforced  composite  laminate  Comparisons  with 
linear  elasticity  indicate  that  viscoelastic  relaxation  may  reduce  the  residual 
stresses  by  about  20  percent.  He  also  presented  an  elastic  analysis  of  the 
stresses  in  adhesive  joints  due  to  moisture  and  temperature  using  variational 
principles  [12].  It  is  observed  that  the  significant  stresses  at  the  interface  are 
essentially  confined  to  the  neighborhood  of  the  outer  edge. 

Weitsman  has  also  analyzed  the  effect  of  moisture  sorption  on  the  mechanical 
behavior  of  viscoelastic  adhesive  layers  [  1 3 ]  .  for  steady  ambient  humidity  the 
viscoelastic  stresses  are  smaller  than  their  elastic  counterparts  In  their  work 


i 


-  5  - 


Abbott  and  Bumpton  [14]  concluded  that  the  physical  properties  of  a 
polyurethane  adhesive  polymer  are  affected  by  exposure  to  moist  environments 
particularly  at  higher  temperatures.  Sargent  and  Ashbee  [15]  carried  out  an 
optical  interference  investigation  of  swelling  due  to  water  uptake  by  adhesives. 
They  found  that  the  swelling  in  their  specimens  was  strongly  inhomogeneous, 
resulting  in  the  presence  of  both  compressive  and  tensile  stresses  normal  to  the 
bond 

In  this  part  of  this  thesis  we  consider  a  model  specimen  which  consists  of 
only  two  layers  of  dissimilar  materials  with  free  edges  (PVA*  and  BK-7  glass). 
The  material  characterization  is  delineated  in  Chapter  2  Next  the  calculation  of 
stresses  and  strains  in  a  bimaterial  specimen  emphasizing  the  time  dependent 
material  properties  of  the  adhesive  is  developed  in  Chapter  3  The  fundamental 
field  equations  for  isothermally  varying  temperature  fields  in  viscoelastic 
materials  have  been  stated  by  Muki  and  Sternberg  '  ’.6],  and  by  Morland  and  Lee 
'.7]  .  these  equations  are  used  to  formulate  the  model  presented  here.  The  final 
step  in  this  analysis  consists  in  the  experimental  examination,  which  parallels 
the  model,  the  measurements  consist  of  the  time-dependent  out-of-plane  defor¬ 
mations  of  thin  sandwiches  composed  of  glass  and  PVAc  through  the  use  of  opti¬ 
cal  interferometry  ;  that  work  is  presented  in  Chapter  4 
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2.  VISCOELASTIC  MATERIAL  CHARACTERIZATION 

This  chapter  is  concerned  with  the  determination  of  mechanical  and  physical 
properties  of  polyvinyl  acetate  pertinent  to  the  stress  analysis.1  Section  2.1  stu¬ 
dies  the  thermal  expansion  behavior  of  PVAe  leading  to  the  measurement  of  its 
temperature  dependent  thermal  coefficient  of  expansion.  In  Section  2.2  the 
creep  behavior  of  polyvinyl  acetate  is  investigated  and  the  shear  relaxation 
modulus  is  then  calculated  numerically  via  the  appropriate  relations  between 
the  material  parameters. 

2.1  Thermal  Expansion  Behavior 

The  thermaL  expansion  characteristics  of  polyvinyl  acetate  were  determined 
using  a  glass  dilatometer  containing  mercury  as  the  confining  liquid  (18, 19]. 8 
This  expansion  contains  a  reversible  component  on  which  other  changes  in 
volume  due  to  variations  in  moisture  content,  to  curing  and  outgassing  of  low 
molecular  weight  components  are  superimposed.  The  method  used  here 
excludes  these  spurious  factors. 

Polymers  display  significant  changes  in  their  mechanical  and  physical  prop¬ 
erties  in  the  glass  transition  range  ;  these  include  the  coefficient  of  thermal 
expansion,  the  heat  capacity  and  the  compressibility.  Outside  the  glass  transi¬ 
tion  region  the  relation  between  volume  change  and  temperature  remains 
nearly  linear  .  the  intersection  of  the  extension  of  these  linear  portions  from 
below  and  above  the  transition  domain  determines  the  position  of  the  glass- 
transition  temperature  Tg.  The  value  of  Tg  will  vary  slightly  with  the  rate  of  cool¬ 
ing  or  heating 

1.  The  average  weight  molecular  weight  of  polyvinyl  acetate  was  indicated  by  the  manufacturer 
(Poiysciencea  inc.,  Warrington.  P^  to  be  500.000. 

2.  in  the  early  stages  of  the  material  characterization  it  was  attempted  to  use  a  10-38  Fotomc 
Sensor  (Mechanical  Technology  Lncj  to  record  the  expansion  of  a  PVAC  sample.  This  non-contact 
optical  displacement  detector  forms  an  application  of  a  fiber  optics  principle.  Due  to  non- Linear 
and  non-repeata'ole  response  to  temperature  changes  in  the  sensor,  its  further  use  was 
discontinued  and  a  mercury  dilatometer  was  constructed  suosequently 


The  algorithm  of  Hopkins  and  Hamming  [29]  is  used  to  determine  the  relaxa¬ 
tion  modulus  of  PVAC  with  the  master  creep  curve  of  Figure  7  as  input  for  the 
calculations.  The  results  are  presented  in  Figure  8. 


It  should  be  noted  however  that  algorithms  available  for  numerically  invert¬ 
ing  viscoelastic  functions  are  usually  illustrated  by  generating  creep  functions 
from  relaxation  data,  i.e.  the  converse  procedure  of  the  one  attempted  here. 
These  calculations,  like  the  one  illustrating  the  procedure  in  [29],  are  successful 
because  errors  in  calculating  the  creep  function  do  not  propagate  and  amplify 
in  subsequent  stages  of  the  computation  We  document  this  assertion  by  con¬ 
sidering  the  recursion  formula  used  in  reference  [29]  ;  values  of  the  creep  com¬ 
pliance  at  successive  time  intervals  are  calculated  using  the  expression 


'•fi  (tn  +  1/2)  = 


t-rn-l  ^ (ti+ 1/2)  |  f(tn+l  tj)  f(tn+]  tj+i) 

_ i=0 _ [ _ 

f(tn«-l  —  tn) 


In  equation  (6)  p(t)  refers  to  normalized  values  of  the  creep  function,  and 

t 

f(t)  =  f  ^  (t)  dr  (7) 

0 

in  which  Y'(t)  represents  normalized  values  of  the  shear  modulus.  In  order  to 
study  potential  sources  of  error  propagation  we  look  at  the  coefficient  of 
?(t„-1/e)  in  the  expression  for  <p{in+\/z)  This  coefficient,  after  application  of  the 
mean  value  theorem,  is 


V*  (tn»l  jO  tn-1 

V'  (Win  ~  t-n+l  —  tn 


where  t,  and  are  values  of  time,  such  that  tn  >  >  tn_1  and  tn+l  >  -iSfe  >  tn, 

and  where  tj  correspond  to  equidistant  time  intervals  on  the  logarithmic  scale 
such  that  tn_j  <  tn  <  tnM  Hence  each  bracket  has  an  absolute  value  less  than 
unity  so  that  the  potential  error  in  V'Un-i/a)  is  multiplied  by  a  number  less  than 
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log  aT  = 


Ci  (T-Tr) 
C2  +  (T-Tr) 


(4) 


The  constants  and  C2  in  equation  (4)  are  found  to  be 
Cx  =  9.49 

C2  =  36.12  for  T  >  Tg  and  Tr  =  40. 

These  values  seem  acceptable  in  light  of  the  data  presented  by  Ferry  [27], 

As  indicated  before,  our  results  complement  imtial  results  obtained  by 
Knauss  and  Kenner  [2],  Comparison  of  their  data  as  well  as  the  present  ones 
with  those  by  Plazek  indicate  a  slight  shift  for  these  results  toward  shorter  time 
Residual  moisture  present  in  the  PVA,.  seems  the  most  reasonable  explanation 
for  this  difference. 

Plazek’s  data  [25]  are  used  to  estimate  the  steady  state  compliance0  J°  the 
limiting  value  is  taken  to  be  log  J°  =  -  3.60  (in  m2/N).  Figure  7  shows  the  master 
curve  for  the  recoverable  creep  compliance  of  PVAC  at  40  °C,  including  the  long 
time  values  from  reference  [25]  ;  the  dashed  line  beyond  log  t  =  9.50  indicates 
this  extrapolation. 


2.2.2  Relaxation  Modulus  of  Polyvinyl  Acetate  As  the  final  step  in  the  viscoelas¬ 
tic  characterization  of  the  model  polymer,  the  shear  relaxation  modulus  /n(t)  is 
calculated  from  the  experimental  data  for  J(t).  The  creep  compliance  and  the 
relaxation  modulus  are  connected  by  the  relations  [27], 

f  J(t-f)M(Odf  =  t 
0 

or  (5) 


/ 


M  t-£)J(0d* 


t 


0.  The  steady  state  compliance  is  deflned  as 

lim  [  J(t)  -  t/77  ]  =  J°  . 

t-*  - 
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heating  degrades  the  material,  it  becomes  nearly  impossible  to  obtain  experi¬ 
mental  data  in  this  temperature  domain  in  a  "finite"  time  interval.  Therefore 
the  upper  temperature  data  are  given  a  ''natural"  extension  towards  a  slope 
equal  to  unity.  This  procedure  allows  an  estimate  of  the  viscous  contribution  in 
the  temperature  range  covered  by  the  measurements. 

Subsequently  horizontal  shifting  is  attempted  relative  to  the  data  for  24  °C, 
obtained  by  Knauss  and  Kenner  [2],  A  smooth  composite  curve  is  obtained  from 
the  individual  curves,  as  shown  in  Figure  5,  where  the  results  for  the  four  test 
runs  involving  several  specimens  are  presented.  The  above  mentioned  "natural" 
extension  is  indicated  in  dashed  line  segments,  namely  beyond  log  t  =  9.50.  It  is 
also  found  that  the  viscous  contribution  in  the  range  of  test  temperatures  is 
essentially  negligible5 

It  is  apparent  from  the  quality  of  the  master  curve  that  PCA*.  can  be  con¬ 
sidered  to  behave  like  a  thermorheologically  simple  material  whose  temperature 
dependence  can  be  described  by  a  shift  factor  ayo(T),  such  that 

J(T.t)  =  (T0/T)  J  T0,[t/hy9(T)]  .  (3) 

This  time-temperature  shift  factor  ax  is  shown  in  Figure  6  for  dry  PVAc,  includ¬ 
ing  error  bars  representing  the  uncertainty  in  shifting  the  individual  curves  of 
Figure  4.  The  different  symbols  correspond  to  the  various  test  runs,  including 
the  results  of  reference  [2].  One  should  note  a  marked  change  in  slope  in  the 
glass  transition  temperature  range.  For  temperatures  above  Tg  the  shift  factor 
ax  can  be  represented  by  the  WLF  equation  [1],  namely 

5.  Recovery  experiments,  associated  with  one  of  the  four  test  runs,  had  failed  to  yield  satisfactory 
results  for  the  value  of  the  viscosity.  The  recovery  times  turned  out  to  be  extremely  long  in  the 
temperature  range  under  consideration.  Another  approach  for  determining  the  value  of  r), 
suggested  by  Ninomiya  [28],  also  fell  short  in  producing  useful  results  ,  this  was  due  to  a  poorly 
defined  extrapolation  of  (dJ(t)/dt)  versus  1/t. 


log  Jp(t) 
(m2/N ) 
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the  geometry  of  the  specimens  was  changed  significantly  with  respect  to  the  one 
used  in  reference  [2].  Short  cylindrical  PVAc-samples  are  glued  with  an  Eastman 
910  adhesive  to  aluminum  holders.  The  dimensions  of  the  cylindrical  PVA* 
specimens  range  from  7.94  to  17.72  mm  in  diameter,  and  from  5.0  to  5.96  mm  in 
length.  For  these  composite  specimens  the  amount  of  strain  in  the  metal  parts 
is  negligible  compared  to  the  strain  in  the  polymeric  central  part. 

Figure  4  displays  the  recorded  creep  compliance  for  dry4  PVAc  as  a  function 
of  time  with  the  temperature  as  a  parameter,  for  one  of  the  four  test  runs.  The 
error  level  for  these  curves  is  evaluated  to  be  approximately  2  %.  From  these 
data  a  master  creep  curve  for  PVA*  is  generated.  At  first  the  curves  are  shifted 
vertically  by  an  adjustment  factor  T/T0.  where  T„  is  taken  to  be  40  °C  ,  hence  the 
temperature  dependence  of  the  rubberlike  nature  of  the  response  is  accounted 
for  [27],  This  shift  is  included  in  Figure  4  where  Jp(t)  represents  the  vertically 
shifted  creep  compliance.  In  order  for  subsequent  horizontal  shifting  to  be  suc¬ 
cessful.  the  importance  of  the  long  time  viscous  contribution  needs  to  be  deter¬ 
mined. 

According  to  the  theory  of  linear  viscoelasticity  the  creep  compliance  of  an 
amorphous  uncrosslinked  polymer  can  be  described  by, 

J(t)  =  Jg  +  Jd(t)  +  -  (2) 

V 

where  Jg,Jd(t)  and  rj  are  the  glassy  compliance,  the  delayed  creep  compliance 
and  the  '  .scosity,  respectively.  In  the  long  time,  high  temperature  domain  J(t)  is 
largely  determined  by  the  viscous  contribution  t  /  r),  indicated  by  ^  linear  por¬ 
tion  in  a  plot  displaying  log  J(t)  versus  log  t  with  a  slope  equal  to  unity.  Because 
it  takes  an  inordinate  amount  of  test  time  to  reach  this  45°  slope  and  long  time 


4,  The  term  "dry"  refers  to  the  in  vacuum  storage  of  the  specimens.  This  does  not  exclude  some 
possible  residual  moisture  which  could  not  be  extracted  from  the  samples  in  our  laboratory. 
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contraction  of  the  polymeric  samples  does  not  yield  repeatable  values  for  Tg. 
The  values  of  the  cubical  thermal  coefficient  of  expansion  are  found  to  be 
2  50x10 ~*/°C  in  the  glassy  domain  and  5.98xlQ"4/°C  above  Tg.  The  error  for 
these  data  was  calculated  to  be  about  3  %.  The  test  results  from  different  runs, 
indicated  by  distinct  symbols  in  Figure  3,  proved  to  be  very  repeatable.  The 
numerical  data  for  this  particular  type  of  PVA*  are  in  good  agreement  with 
results  previously  obtained  for  PVA*  by  Kovacs  [23]. 

2J2  Shear  Creep  Behavior 

In  order  to  complete  the  necessary  material  characterization  of  the  model 
material  PVA,.,  its  shear  creep  compliance  is  determined.  The  testing  was  car¬ 
ried  out  on  a  creep  torsiometer,  as  outlined  in  reference  [24].  These  authors 
also  report  creep  data  for  PVAC  in  the  temperature  region  around  the  glass  tran¬ 
sition  temperature  [2]. 3  These  initial  results  are  here  extended  in  the  upper 
temperature  region,  so  as  to  account  also  for  the  possible  viscous  flow  contribu¬ 
tion  of  this  uncrosslinked  polymer.  This  latter  determination  is  necessary  to 
assess  the  residual  stresses  due  to  cooldown  from  the  melt. 

2.2.1  Experimental  Results  and  Analysis  The  creep  compliance  J(t)  can  be 
derived  from  the  experimentally  measured  twist  angle  i?  (t)  as  a  function  of 
time,  through  the  relation  for  torsion  of  circular  cross-sections,  given  by 


J(t)  = 


TTd4  tf(t) 


where  d,  L  and  T  are  the  specimen  diameter,  specimen  gauge  length  and  applied 
torque  respectively.  The  derivation  of  equation  (1)  is  trivial  and  is  found  in 
Appendix  1.  In  order  to  guarantee  the  accuracy  of  the  measurements  at  higher 
temperatures  and  due  to  the  limitations  of  the  torsiometer  described  in  [24], 

3.  Complementary  work  on  PVA,,  has  been  reported  by  Plazek  [25,26]  in  a  broader  temperature 
range. 
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tubing  through  which  liquid  nitrogen  can  flow  for  cooling  purposes  ;  a  thermo¬ 
couple  monitors  the  temperature  of  the  specimen  environment.  The  uniformity 
of  the  temperature  throughout  the  bath  is  maintained  by  two  continuously  mov¬ 
ing  propellers.  Specimen  volume  changes  result  in  changes  in  mercury  level,  in 
the  capillary,  which  are  read  using  the  ruling  on  the  capillary. 

For  accurate  results  a  calibration  is  necessary  to  account  for  both  the  non¬ 
uniformity  of  the  capillary  bore  and  the  intrinsic  thermal  expansivity  of  the 
apparatus.  The  former  calibration  was  made  prior  to  making  a  U-bend  in  the 
capillary.  With  a  stop-cock  affixed  at  the  lower  end  of  the  capillary,  mercury  was 
drawn  by  a  vacuum  into  the  capillary.  The  mercury  level  was  then  lowered  in 
small  increments  by  opening  the  stop-cock  ;  the  corresponding  volume  of  the 
capillary  was  determined  by  weighing  the  removed  mercury  with  a  Mettler  HL32 
analytical  balance  with  a  readibility  range  of  0.1  mg,  and  then  using  standard 
tabulated  data  for  the  density  of  mercury.  The  average  volume  of  the  tube  per 
height  increment  was  found  to  be  0.00013  ml/mm  ;  any  nominal  volume  above 
reference  level  computed  using  this  average  can  thus  be  corrected. 

The  calibration  for  thermal  expansion  effects  involves  measuring  the 
apparent  volume  change  with  temperature  while  the  specimen  chamber  is  com¬ 
pletely  filled  with  mercury.  The  change  in  volume  with  temperature  of  a  volume 
of  mercury  equal  to  the  nominal  specimen  volume,  yields  the  volume  change 
that  has  to  be  subtracted  from  the  recorded  data  when  a  specimen  is  present. 

2.1.2  Experimental  Results  The  data  of  the  dilatometric  measurements  are 
presented  in  Figure  3.  The  thermal  expansion  characteristics  for  dry  PVA* 
display  the  change  of  mercury  level  in  the  capillary  versus  the  temperature  for  a 
cooling  rate  of  approximately  4  °C  per  hour  from  which  the  glass  transition  tem¬ 
perature  Tg  is  determined  at  29.0  °C  At  faster  cooling  rates  the  non-equilibrium 
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An  extensive  phenomenological  study  of  the  glass  transition  in  amorphous 
polymers  and  in  particular  in  PVAc  has  been  carried  out  by  Kovacs  [20].  Other 
glass  dilatometry  measurements  yielding  very  satisfactory  results  include  the 
findings  of  Rush  [21,22].  Both  these  authors  follow  the  method  of  quenching  the 
dilatometer  from  a  reference  temperature  into  baths  at  some  other  tempera¬ 
ture.  From  the  equilibrium  volumes  at  the  bath  temperature  the  thermal  expan¬ 
sion  data  are  then  calculated. 

2.1.1  Specimen  Preparation  and  Experimental  Method  The  polyvinyl  acetate 
obtained  from  the  manufacturer  in  the  form  of  clear  spherical  beads,  is  molded 
at  a  temperature  of  125  °C  in  a  steel  mold.  The  details  of  the  preparation  are 
reported  in  reference  [2].  Cylindrical  specimens  about  35  mm  long  and  6.35  mm 
m  diameter,  are  used  in  the  dilatometric  measurements. 

The  design  of  the  apparatus  was  guided  to  a  large  extent  by  the  description 
outlined  by  Bekkedahl  [18].  The  essential  parts  of  the  dilatometer  are  shown  in 
Figure  1.  A  specimen  chamber  is  formed  from  the  laterally  expanded  ends  of 
two  Pyrex  glass  tubes  which  are  clamped  together  via  an  0-ring  seal.  The  upper 
tube  proceeds  to  a  stopcock  and  subsequently  terminates  at  a  reservoir.  The 
lower  tube,  a  capillary  of  nominally  0  5  mm  inside  diameter,  undergoes  a  U-bend 
and  then  ascends  vertically  for  approximately  700  mm  along  a  rule  with  a  least 
count  of  0.5  mm.  These  components  are  supported  on  a  lucite  plate,  as  shown  in 
Figure  2. 

The  specimen  is  placed  in  the  specimen  chamber  and  the  chamber  clamp 
secured.  The  remaining  volume  of  the  dilatometer  is  filled  with  mercury  from 
the  reservoir  via  the  stop-cock,  which  is  then  closed  to  exclude  any  air  bubbles. 
The  specimen  chamber  is  submerged  in  a  silicone  bath  contained  in  a  Dewar 
flask  The  latter  contains  an  inductance  heating  element  and  a  coiled  copper 


Figure  B.  Master  relaxation  curve  for  dry  PVAc  at  40  °C. 
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unity  in  size  if  the  input  function  Y'(t)  is  a  monotonically  decreasing  function. 
Hence  the  effect  of  the  error  dies  out. 

In  the  present  case,  where  it  is  attempted  to  generate  the  shear  modulus 
from  creep  data,  the  error  propagation  becomes  unstable  because  the  input 
function  is  a  monotonically  increasing  function  and  because  of  the  steep  slope 
of  the  master  creep  curve  of  PVA,.  in  the  transition  range.  This  consideration 
becomes  apparent  by  evaluating  the  ratio  between  the  first  pair  of  brackets  in 
expression  (8).  With  the  creep  compliance  as  input  function,  the  product  of  the 
first  term  is  no  longer  less  than  unity.  Nevertheless  using  time  intervals  whose 
amplitude  was  chosen  on  the  basis  of  [ 29 J  it  was  tried  to  invert  the  creep  com¬ 
pliance.  In  the  single  precision  mode  the  error  growth  was  unstable  ;  therefore 
the  computations  were  executed  in  double  precision.  The  error  oscillations  did 
not  completely  die  out,  but  a  solution  throughout  the  whole  time  spectrum  was 
nevertheless  obtained.  Consequently  an  essentially  trial  and  error  inversion 
procedure  was  finally  developed.  The  domain  where  the  inital  solution  for  fi( t) 
did  not  appear  smooth  is  approximated  by  a  smooth  curve.  Then  as  a  check, 
the  numerically  obtained  shear  relaxation  modulus  was  itself  again  inverted  to 
insure  that  the  experimentally  measured  creep  compliance  was  recovered.  In 
these  calculations  the  choice  of  the  time  intervals  seemed  the  critical  factor  : 
with  an  appropriate  choice  of  this  parameter,  very  few  iterations  were  neces¬ 
sary.  At  long  time  values,  i.e,  beyond  log  t  =  9 .50,  the  data  are  extended  in 
accord  with  the  procedure  used  earlier  for  the  creep  data. 


3.  VISCOELASTIC  STRESS  ANALYSIS 

Having  obtained  the  pertinent  material  properties  for  the  viscoelastic 
material  we  proceed  to  the  stress  analysis  of  a  bimaterial  specimen,  composed 
of  an  elastic  stiff  substrate  and  an  attached  viscoelastic  layer.  The  theoretical 
results  obtained  next  will  be  compared  with  experimental  results  in  Chapter  4. 
Section  3.1  outlines  the  relations  derived  to  compute  the  stresses  as  a  function 
of  time.  The  second  part  of  this  chapter,  namely  Section  3.2,  presents  the 
numerical  results  of  the  modelling. 

3.1  Formulation  of  the  Problem 

It  is  assumed  that  a  perfect  bond  exists  between  the  two  dissimilar  materials 
involved.  Furthermore  we  consider  only  slow  thermal  changes  such  that  the 
sample  does  ot  experience  a  temperature  gradient  through  its  thickness  nor 
incurs  any  non-equilibrium  volume  changes.  The  analysis  allows  only  for  small 
deformations. 

The  geometry  of  the  specimen  is  shown  in  Figure  9.  The  viscoelastic  material 
is  modelled  by  a  thermorheologically  simple  material  ;  its  mechanical  properties 
include  the  shear  relaxation  modulus  /i(t)  and  its  temperature  dependent  shift 
factor  aTo(T),  the  bulk  modulus  K„,  taken  to  be  a  constant,  and  the  thermal 
coefficient  of  expansion  a(T).  The  latter  is  assumed  to  be  simply  a  function  of 
the  temperature,  i.e.  no  explicit  time  or  temperature  rate  dependence  is  allowed 
[30].  Comparison  with  experimental  results  later  will  examine  the  validity  of 
this  assumption.  The  mechanical  behavior  of  the  adhering  substrate  is  charac¬ 
terized  by  corresponding  elastic  and  physical  parameters. 

Let  a  displacement  field  be  defined  for  both  material  domains  as  : 


Adhesive  (viscoelastic) 


Substrate  (elastic) 
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3w 

3r 


v  =  0 


w  =  w(r) 


(9) 


where  u.  v.  and  w  are  the  displacements  in  the  r,  and  z-direction  respectively  ; 
the  coordinates  are  defined  in  Figure  9  and  e0  is  the  radial  strain  of  the  bottom 
fiber  of  the  bimaterial  element.  Let  the  curvature  be  given  by  : 


1 

R 


3zw 

ar2 


(10) 


where  w(r)  denotes  a  spherical  deflection,  whose  magnitude  is  supposed  to  be 
small  enough  to  have  equation  (10)  hold  near  r  =  0.  Our  immediate  purpose 


here  is  to  calculate  the  strain  e0  as  well  as  the  curvature 


result  of  a  changing  temperature. 


3zw 

3rz 


occurring  as  a 


The  stress  in  the  z-direction  is  zero  so  that  the  only  non-vanishing  stresses 
are  ar  =  ~  a  ,  while  az  =  0.  Furthermore  all  shear  strains  vanish  (except  near 

an  edge).  Under  non-isothermal  conditions  : 


O'ldc  =  3Ko 


ericic  -  3ai3 


(11) 


where  indicates  the  temperature  difference  with  respect  to  the  reference  tem¬ 
perature.  Also, 

Sy  =  2fj.  •  dejj  $  (12) 

where  Sy,  e^  are  deviator  stresses  and  strains,  respectively.  For  the  present 


t  Repeated  indices  denote  the  summation  convention. 

t  The  star  notation  (♦)  refers  to  Stieltjes  convolution  ;  equation  (12)  should  read  in  full  as  : 

Sjj  =  2 /  p(t-{)deqtt)  ■ 
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geometry,  we  have 


ffkk  =  ar  +  a*  -  2or 


Sr  =  far 


1  ,  2 aT 

e'  3  T(  HCt3w>) 


Hence  equation  (12)  becomes: 


ar(z,t)  =  6/x(t)  •  d  [  er(t)  - 


2crr(z,t) 


-a(T(t))*(t)] 


gr(z,t)  4-  •  d<xr(z.t)  =  8/u(t)  •  d[er(t)  -  a(T(t))  (t)].  (15) 

If  the  temperature  varies  during  the  stress  history  one  has,  upon  applying 
instantaneous  temperature  reduction  [16], 

df  =  f  16} 


C-  f  4£— 
?  J0  aT[T(C)] 


(17a) 


t  '  =  / 


aT[T(0] 


(17b) 


Thus  equation  (15)  can  be  written  as, 

4  r  8*r(*.f)  „r  .,x  a[er(0-a(T(0)tf(f)]  ,.„x 

CTr(Z,t)  3k7-(  M(t  -f) — ^ - d£  =  6 J  M t'-f) - ^ - 3*  (18) 


Because  by  assumption  the  temperature  T  does  not  depend  on  z  (nor  on  r  and  i$) 
and  since  the  temperature  reduced  time  t'  is  a  function  of  t,  any  explicit 


and  consequently 
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• 

eTi  (z.t’)  =  18K„  M(t')  •  dPi(t')  •  d[fr(z,t*)  -ff(t')] 

(26) 

'a 

% 

where  the  subscript  1  refers  to  the  range  of  the  viscoelastic  material 

in  which 

\ 

this  expression  for  the  aT  stress  is  valid. 

In  a  similar  fashion  the  stress  in  the  elastically  behaving  substrate  can  be 

• 

computed.  The  elastic  constitutive  equations  equivalent  to  equations 

(11)  and 

(12)  are  : 

oitit  =  3K  [sjd,  -  3onJ] 

(27a) 

• 

Si]  =  2/zeij. 

(27b) 

Hence, 

• 

(7r(z,t)  =  Qfj.  [er(t)  QK  a(T(t))  tf(t)] 

(28) 

or 

» 

]-i 

<7r(z,t)  =  SjLt  1  +  [er(t)  -a(T(t))  *(t)] 

(29) 

Using  the  appropriate  relations  between  the  elastic  moduli,  and  substituting  the 

• . J 

expression  for  the  strain,  namely 

■« 

fr(t)  =  e,(t)  -  z^-f-(t) 

3r* 

(30) 

» 

into  equation  (29)  yields  analogously  to  equation  (26), 

»rt(z.t)  =  7^T-[£o(t)  -z~-  (t)  -aa(T(t))tf(t)] 

*  1  -i/2  dr* 

(31) 

• 

where  the  subscript  2  refers  to  the  range  of  the  elastic  substrate. 

Next,  equilibrium  of  forces  and  moments  can  be  considered.  This 

amounts  to  the  requirements  : 

• 

— ■  -J.  •-»  La  ■  -■ 
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2(hj  *  h,) 

/  <fr(z.t')dz  =  0  (32a) 

O 

and 

2 Oil  *■  hg) 

J  z  <7r(z.t')dz  =  0  (32b) 

0 

Executing  the  indicated  integration  over  the  thickness  of  the  bimaterial  speci¬ 
men  yields  the  following  pair  of  equations  : 

E, 

18K.h,Mi(t’)  *  dP,(t')  •  de0(t')  + -r— -h2  f0(t’) 

I  —1/2 

-  l8K0(hf  +  2h1h*)Mi(f)  •  d^(f)  *  d0-(f)  +  0-^’) 

E 

=  lSKoh^itt')  •  dPj(t’)  .  dn  j(t')  +  — ^ —  hz  a2S(t‘)  (33a) 

1  ~l/2 


|18Ko(hf +  2h1hg)#*i(t')  •  dPi(t')  *  df0(f)  +  ^_h|e0(t') 


!  a  9  o  a 

-I  24K0(hi  +  3h!ha  +  3hih8)A4i(t')  *  dPi(t’)  *  d  — —  (t') 

dr 


E2  4,3  d2vt  ,  ,x 


Ep  .2  , 


=  10Ko(hr  +  2ht  haWt’)  .  dP,(f)  *  dtf  i(t’)  +  h2a2  3(f)  .  (33b) 

1  — I/a 


This  system  of  two  equations  (33a)  and  (33b),  in  two  unknowns  e0  (t')  and 


2——  (t‘).  yields  solutions  in  the  temperature  reduced  time  domain.  They  are 
dr2 


related  to  their  counterparts  in  real  time  through  expressions  similar  to  equa¬ 
tions  (19),  namely, 
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fo(t')  =  e0(t) 


d2w 

Sr2 


(t')  = 


a2w 

Sr2 


(t). 


(34a) 

(34b) 


Through  relations  (17)  the  solution  in  the  real  time  domain  t  can  be  computed. 

In  the  limit  of  two  Linearly  elastic  solids  the  analysis  simplifies.  Those 
limit  computations  are  outlined  in  Appendix  2.  In  temperature  domains  far 
removed  from  the  glass  transition  temperature  range  the  viscoelastic  coating 
acts  essentially  like  a  purely  elastic  medium.  Hence  the  glassy  and  long  time 
properties  of  polyvinyl  acetate  can  be  used  to  calculate  the  elastic  asymptotes, 
characterizing  the  behavior  of  the  bimaterial  layer  in  the  corresponding  tem¬ 
perature  domains. 


3.2  Numerical  Results 


Two  cases  are  considered  ;  the  first  concerns  the  heating  of  a  bima- 
terial  sample  from  temperatures  well  below  the  glass  transition  of  the  polymer 
up  to  the  temperature  region  where  the  viscoelastic  material  essentially  acts 
like  a  rubbery  solid.  The  second  case  deals  with  the  cooling  of  a  similar  speci¬ 
men  between  the  same  extreme  temperature  domains,  this  time  however  start¬ 
ing  in  the  melt  region.  For  each  of  the  two  examples  three  graphs  are  presented 
below  The  first  displays  the  graphical  representation  of  the  relation  between  the 
real  time  t  and  the  temperature  reduced  time  t'  as  given  by  equation  (17).  The 
next  two  figures  represent  the  solutions  of  equations  (33)  for  the  radial  strain  in 
the  bottom  fiber  of  the  specimen  ci(t),  where 


Sl  =  *o(t)  -  a2(T-T0) 


(35) 


and  the  radius  of  curvature  — —  (t) 

dr4 


It  should  be  noted  that  the  radial  strain  of 


the  specimen  is  referred  to  the  bottom  fiber  (z=Q)  and  not  to  the  neutral  axis  or 


» 


plane  of  the  bimatertal  layer  assembly  ;  this  formulation  was  preferred  because 
the  location  of  the  neutral  axis  or  plane  changes  continuously  -  though  not 
much  -  due  to  the  varying  material  properties  as  a  function  of  time  throughout 
the  thermal  loading  history. 

Some  mechanical  and  physical  properties  of  polyvinyl  acetate  still  need 
to  be  defined  From  the  glassy  shear  modulus  and  a  Poisson’s  ratio  v  -  0.41 
32],  the  bulk  modulus,  taken  to  be  a  constant,  equals  Ko  =  8.28x10®  N/m2  This 
value  of  K„  combined  with  the  experimental  rubbery  shear  modulus  /U.  yields  v 
-  0.50  in  the  rubbery  limit.  The  properties  of  BK-?  glass  included  in  these  two 
sample  calculations  are  [31]  : 

a  =  7  IxlO'VC 

E  =  B.lxlO9  N/m2 

u=  0  208 

In  both  cases  the  temperature  changes  with  a  constant  rate  of  ±  0.050  °C/min. 
Finally  the  thicknesses  are  hj  =  0.5  mm  for  the  polymer,  hz  =  0.3  mm  for  the 
glass  substrate  respectively. 

3  2.1  Numerical  Evaluation  Procedure  The  solution  of  the  system  of  equations 
(33)  is  carried  out  in  several  steps  listed  below  : 

1.  First  the  relation  between  the  real  time  t  and  the  temperature  reduced  time 
scale  t‘  is  computed. 

2  Subsequently  the  system  of  two  equations  (33)  is  solved  in  the  t'  domain.  This 
computation  involves  determining  the  various  material  properties  at  the  succes¬ 
sive  time  steps  at  which  the  convolution  integrals  need  to  be  evaluated. 

3  The  final  step  consists  in  solving  a  system  of  equations  using  Cramer’s  rule  ; 
the  solution  in  the  real  time  domain  t  is  then  obtained  through  equations  (17) 


3.2  2  Heating  Case  The  results  pertinent  to  the  heating  case  are  shown  in  Fig¬ 
ures  10,11  and  12.  The  relation  between  t  and  t‘  yields  a  monotonically  increas¬ 
ing  function  resulting  in  an  unambiguous  one  to  one  correspondence  between 
the  two  time  scales  involved.  Hence  the  solution  of  equations  (33)  can  be  readily 
reduced  to  real  time,  once  it  is  known  on  the  temperature  reduced  time  scale  t’. 

fiZyf 

The  solutions  ejt)  and  — r-(t)  exhibit  a  very  similar  qualitative  behavior  and 

dr 

can  thus  be  discussed  concurrently. 

(1)  At  the  onset  of  the  heating  at  T  =  -  25  °C  (Tg  =  29  °C)  the  bimaterial  speci¬ 
men  acts  essentially  like  a  sandwich  with  two  elastic  media.  For  purposes  of 
comparison  the  elastic  asymptotes,  computed  using  equations  (A2.2)  and  (A2.3), 
are  also  drawn  in  Figures  1 1  and  12  ;  in  these  the  polymer  is  characterized  by  its 
glassy  properties. 

(2)  Approaching  the  glass  transition  region  one  first  notices  a  slight  increase  in 
the  strain  and  the  curvature  over  the  values  predicted  by  the  glassy  elastic 
analysis  ;  this  feature  reflects  the  increase  in  the  coefficient  of  thermal  expan¬ 
sion  of  PVAe  before  the  material  relaxation  occurs. 

(3)  A  few  degrees  (5  to  7)  above  the  glass  transition  temperature  (Tg  =  29  °C) 
relaxation  the  polymer  becomes  the  dominant  phenomenon  causing  a  sharp 
reversal  in  the  slopes  of  both  curves. 

(4)  Once  the  glass  transition  range  has  been  crossed  the  specimen  gradually 
starts  acting  again  like  a  purely  elastic  bimaterial  layer  with  the  polymer  now 
behaving  in  a  rubberlike  fashion.  Eventually  the  solutions  to  converge  to  the 
elastic  asymptotes  calculated  using  the  results  outlined  in  Appendix  2.  This 
correspondence  is  quite  rapidly  achieved  for  the  axial  strain  rL(t)  ;  for  the  cur- 
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vature  — 5-  (t)  the  convergence  takes  a  little  longer. 


Temperature  V 


ceeding  through  the  glass  transition  region  of  PVA*  to  higher  temperatures 
il  the  radius  of  curvature  seems  to  reach  a  limit  value,  i.e.  the  fringe  spacing 
erved  on  the  TV  monitor  appears  to  remain  unchanging.  Some  examples  of 
uced  experimental  data  a*re  shown  in  Figures  20  and  22  which  display  the 
ius  of  curvature  R,  calculated  from  the  observed  fringe  patterns,  as  a  func- 
1  of  time  and  temperature.  The  calculations  leading  to  the  displayed  values 
R  are  performed  using  a  least  square  fit.  Figures  21  and  23  display  the 
responding  curvature  1/R.  The  initial  values  of  R  and  1/R  are  "negative"  in 
ord  with  our  definition  of  the  "sign"  of  the  displayed  curvature. 

It  should  be  pointed  out  at  this  time  that  the  magnitude  of  R  or  1/R  in 
Figures  20  through  27  represent  averaged  "instantaneous"  values  for  the 
ius  of  curvature  or  for  the  curvature  at  the  corresponding  time.  One  would 
iect  that  with  continuously  and  slowly  varying  temperature  the  curvature 
inges  also  correspondingly  slowly  and  continuously.  Surprisingly  this  was  not 
case,  but  in  a  certain  temperature  range  the  curvature  would  increase  or 
:rease  steadily  with  a  superposed  cyclic  excursion  about  the  slowly  changing 
an. 

In  each  of  the  test  runs  corresponding  to  Figures  20  through  23,  a 
niftcant  change  in  R  or  1/R  can  be  observed  when  the  glass  transition  zone  is 
issed.  The  samples  are  eventually  heated  up  to  65  °C  or  58  °C  respectively  ;  at 
ise  temperatures  the  recorded  fringe  pattern  no  longer  seems  to  vary  as  a 
iction  of  the  prescribed  temperature  history. 

After  the  polymer  relaxes  at  these  high  temperatures  from  6.5  to  10 
irs  the  samples  are  cooled  down  to  room  temperature  conditions.  It  turns 
.  that  the  radius  of  curvature  does  not  drop  its  magnitude  substantially  with 
:reasing  temperature.  When  approaching  the  glass  transition  zone  from 
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recorded  with  two  Chromel-Alumel  thermocouples.  One  of  them  is  fixed  in  the 
chamber  right  above  the  sample  ;  the  other  one  is  embedded  in  a  bimatenal 
specimen,  made  out  of  BK-7  glass  and  Solithane  113,  and  mounted  right  next  to 
the  glass  -  PVAC  specimen  being  tested.  This  procedure  allowed  monitoring  the 
temperature  in  the  chamber  and  in  a  specimen  similar  to  the  one  being  studied 
simultaneously  ;  the  difference  in  reading  between  the  two  thermocouples  aver¬ 
aged  1  °C  ;  consequently  the  output  from  the  embedded  thermocouple  is  taken 
as  the  indicator  of  the  temperature  in  the  neighborhood  of  the  tested  specimen 
unless  specified  otherwise.  The  moisture  content  inside  the  environmental 
chamber  is  measured  by  the  wet  bulb  method.  In  addition,  dry  air  is  continu¬ 
ously  pumped  into  the  chamber  to  prevent  moisture  from  building  up. 

4.3  Experimental  Results 

One  series  of  experiments  involves  samples  composed  of  a  0.25  mm 
(0  010")  thick  glass  slide  and  a  polymer  layer  about  0.50  mm  (0.025”)  thick.  The 
other  specimens  consist  of  a  0.50  mm  (0.020")  thick  glass  slide  and  a  0.35  mm 
(0  014")  polymer  coating. 

At  this  point  we  need  to  specify  the  "sign"  of  the  initial  curvature  of  the 
specimen  at  the  beginning  of  an  experiment.  We  define  a  curvature  to  be  "posi¬ 
tive"  if  has  its  center  on  the  positive  z-axis  as  defined  in  Figure  9.  It  was  verified 
experimentally  that  the  pristine  glass  slides  and  the  bimaterial  specimens  actu¬ 
ally  display  a  "negative"  curvature  at  room  temperature  this  observation 
explains  the  shape  of  the  sample  shown  in  Figure  18  The  "sign"  of  curvature  can 
be  checked  by  considering  the  reflected  parallel  light  beam  coming  of!  the  sur¬ 
face  under  investigation.  If  we  can  locate  a  focal  point,  the  specimen  displays  a 
"positive"  curvature  with  respect  to  the  incoming  light  beam. 

Bimaterial  specimens  are  heated  starting  from  ambient  conditions, 
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Figure  18.  Measurement  system  for  observing  Newton's  rings. 


rings  are  then  given  by  : 


r  =  Vm  X  f?  m  =  0,1,2,...  (38) 

A  remaining  point  of  interest  concerns  the  intensity  of  the  central 
spot,  corresponding  to  the  point  of  contact  between  the  convex  surface  and  the 
reference  plane.  It  is  known  that  a  relative  phase  change  of  rr  occurs  under 
reflection  at  the  air-to-glass  surface  [33]  .  If  there  were  no  such  phase  changes, 
the  rays  reflected  from  the  two  surfaces  in  contact  should  be  in  the  same  phase 
and  produce  a  bright  spot  a  the  center.  Hence  when  the  contact  is  perfect,  the 
central  spot  in  the  reflected  light  is  black. 

4.2.2  Experimental  Setup  A  schematic  of  the  experimental  setup  is  shown  in 
Figure  18  The  curved  bimaterial  sample  rests  on  an  optically  flat  surface.  This 
wedge  shaped  support  ls  made  out  of  BK-7  glass  and  has  a  wedge  angle  of  10° 
which  feature  prevents  interference  with  reflections  from  the  front  surface  of 
the  wedge.  The  monochromatic  light  beam  of  a  He-Ne  laser  (X  =  6328  A)  is 
directed  towards  the  specimen  through  a  beam-splitter,  which  also  collects  the 
resulting  interference  pattern  from  the  reflected  beam.  The  Newton's  rings  are 
observed  through  a  Nikon  Microscope  SMZ-10  connected  to  a  Hitachi  CCTV  Cam¬ 
era  (Model  HV-16SU)  throughout  the  experiment.  The  recorded  data  can  be 
viewed  continuously  on  a  Conrad  TV  Monitor  while  they  are  also  stored  on  mag¬ 
netic  tape  on  a  Sanyo  Video  Tape  Recorder  for  later  reference.  An  example  of  a 
recorded  fringe  pattern  is  shown  in  Figure  19. 

The  specimen  and  its  support  are  enclosed  in  an  environmental 
chamber,  the  temperature  of  which  can  be  controlled  to  ±  0.5  °C.  Nearly  con¬ 
stant  heating  rates  could  be  imposed  fairly  readily  ,  cooling  to  temperatures 
below  ambient  conditions  in  a  steady  fashion  posed  more  difficulties  due  to 
problems  in  regulating  the  flow  of  nitrogen  to  the  chamber.  The  temperature  is 
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Figure  17.  Interference  method. 


squeeze  out  any  excess  polymer  material.  After  another  20  minutes  the  press  is 
deactivated  and  the  specimen  is  allowed  to  cool  in  its  mold  under  ambient  con¬ 
ditions  for  several  hours. 

Due  to  this  cooldown  residual  thermal  stresses  are  introduced  in  the 
bimaterial  sample.  Hence  at  room  temperature  the  composite  specimen  displays 
a  curvature.  Apart  from  that  the  glass  slides  themselves  display  an  initial  cur¬ 
vature,  which  amplitude  must  be  accounted  for  in  the  data  reduction. 

4.2  Experimental  Method 

4.2.1  Interference  Method  The  optical  method  used  to  measure  the  out-of-plane 
displacement  of  the  bimaterial  specimen,  is  based  on  interference  of  light  waves 
(Newton’s  rings).  A  schematic  of  the  optical  setup  used  in  the  interference 
method  is  shown  in  Figure  17.  The  reflected  light  contains  a  set  of  fringes, 
resulting  from  either  constructive  (bright  fringes)  or  destructive  (dark  fringes) 
interference  ;  for  the  geometry  shown  in  the  schematic  the  condition  for  a 
minimum  intensity  is  : 

d  =  m  cos  (36) 

where  m  =  0,1,2,...  denotes  the  order  of  the  dark  fringes  ;  X  and  are  the 
wavelength  of  the  used  light  and  the  incident  angle  respectively.  For  normal 
incidence  each  consecutive  dark  fringe  corresponds  to  a  change  in  the  out-of¬ 
plane  displacement  equal  to  half  a  wavelength  of  the  incident  light.  At  a  given 
fringe  of  order  m,  counted  from  the  center  fringe,  the  air  film  thickness  along 
that  fringe  is  given  by  : 

d  =  m  7  m  =  0, 1,2,3,...  (37) 

If  the  fringes  are  produced  in  the  air  film  between  a  spherical  surface  and  an 
optically  flat  surface,  the  resulting  fringes  will  be  circular.  The  radii  of  the  dark 


4.  EXPERIMENTAL  EXAMINATION  OF  TME  DEPENDENT  DEFORMATION 

This  last  chapter  of  Part  1  concerns  measurements  related  to  the 
verification  or  examination  of  the  viscoelastic  stress  analysis  of  the  bimaterial 
specimen  consisting  of  a  glass  slide  and  a  PVAc  layer.  Section  4.1.  describes  the 
preparation  of  the  specimens.  In  Section  4.2  the  experimental  method  is  out¬ 
lined  ;  in  the  subsequent  sections  the  results  are  presented,  as  well  as  a  com¬ 
parison  between  the  experimental  data  and  the  theoretical  calculations  of  the 
previous  chapter. 

4.1  Test  Materials  and  Specimen  Preparation 

In  view  of  our  desire  to  measure  deformations  optically  it  is  preferable 
to  work  with  glass  adherends,  although  that  is  a  not  a  necessary  requirement 
However,  a  clear  adherend  allows  direct  examination  of  the  interface  integrity. 
Circular  slides  of  BK-7  glass  were  then  used  having  a  diameter  of  25.4  mm  (1.0") 
and  are  either  0.25  mm  (0.010")  or  0.50  mm  (0  020”)  thick. 

A  bimaterial  specimen  and  the  polymer  film  deposition  fixture  are 
shown  in  Figure  16.  In  the  latter  the  thickness  of  the  film  is  controlled  by  metal 
spacers.  In  order  to  prevent  the  polymer  from  sticking  to  the  spacers,  the  glass 
slides  are  enclosed  by  a  Teflon  ring  of  the  same  thickness  as  the  spacers.  In 
addition  the  top  and  bottom  sides  of  the  mold  are  sprayed  with  KRAXO  1711 
Release  Agent  (Contour  Chemical  Company)  to  avoid  potential  adhesion  prob¬ 
lems  when  removing  the  specimen  from  the  mold.  The  PVAC7  grains  obtained  as 
such  from  the  manufacturer,  are  poured  on  to  glass  slide,  before  th  -  top  plate 
of  the  mold  is  secured.  With  the  mold  still  not  entirely  closed,  it  is  placed 
between  the  two  plates  of  a  hot  press  at  about  150  °C.  After  heating  for  about  20 

minutes,  enough  pressure  is  applied  to  close  the  film  deposition  fixture  and  to 

7.  The  average  weight  molecular  weight  of  polyvinyl  acetate  was  Jidicated  by  the  manufacturer 
(Poiysciences  Lnc.,  Warrington,  PA)  to  be  500.000. 


».  Curvature  1/Rasa  function  of  time  (cooling) 


Figure  13.  Relation  between  t  and  t’  in  the  cooling  case. 


3.2.3  Cooling  Case  Figures  13,14  and  15  portray  the  results  from  the  case  of  a 
cooldown  from  a  temperature  T  =  150  °C  towards  the  glassy  region  of  the  viscoe¬ 
lastic  layer.  The  relation  between  the  two  different  time  scales  t  and  t‘.  indi¬ 
cated  in  Figure  13,  shows  a  noticeable  convergence  to  an  almost  constant  value 
of  t'.  corresponding  to  higher  values  of  t  and  Low  temperatures.  This  particular 
behavior  clearly  exemplifies  the  "freezing-in"  of  the  temperature  reduced  time 
t‘.  In  order  to  ensure  the  one-to-one  reciprocity  between  t  and  t',  the  calcula¬ 
tion  of  expressions  (17)  has  to  be  carried  out  in  the  double  precision  mode  ; 
even  then  a  rescaling  of  the  values  of  t'  close  to  the  asymptotic  limit  of  t' 
remains  necessary  in  order  to  have  a  readily  applicable  correspondence  between 

t  and  t‘  available.  As  in  the  previous  example,  the  solutions  ®i,(t)  and  — r-  (t) 

dr * 


behave  very  similarly. 

(1)  During  the  initial  cooling  in  the  rubbery  region,  the  bimaterial  specimen  acts 
like  an  elastic  composite  specimen  ;  the  behavior  matches  the  elastic  asymp¬ 
totes  as  indicated  in  Figures  14  and  15. 

(2)  Upon  nearing  the  glass  transition  zone  from  above  the  freezing-in  of  the 
material  properties  starts  far  above  Tgl  in  this  case  around  75  °C.  This  consti¬ 
tutes  a  rather  surprising  result,  as  it  occurs  about  45  °C  above  Tg  for  this  partic¬ 
ular  cooling  rate. 

(3)  After  passing  through  the  glass  transition  zone  of  PVA,,  around  29  °C.  the 

d^yf 

rate  of  increase  in  e0(t)  and  — r-  slows  down  to  converge  to  the  elastic  limits, 

dr* 


signaling  glassy  behavior  of  the  polyvinyl  acetate  layer. 


Figure  12.  Curvature  1/Rasa  function  of  time  (heating). 
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Figure  20.  Radius  of  curvature  of  a  bimaterial  specimen  as  a  function  of  time 
and  temperature  (thin  glass  slide). 


Temparalur*  (*C) 


Timm  (mm  ) 


Timm  ( min  ) 


Figure  21.  Curvature  of  a  bimaterial  specimen  as  a  function  of  time  .and  tem 


perature  (thin  glass  slide). 
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Figure  23.  Curvature  of  a  bimaterial  specimen,  as  a  function  of  time  and  tem 
perature  (thick  glass  slide). 
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above  a  complicated  fringe  pattern  appears  which  is  definitely  not  axisym- 
metric.®  The  complicated,  if  not  incoherent,  fringe  pattern  which  can  be 
observed  seems  to  transform  into  one  commensurate  with  cylindrical  bending. 
Several  degrees  below  the  glass  transition  temperature  (Tg  =  29  °C)  the  ’’regular” 
Newton's  rings  reappear, That  fact  is  indicated  by  the  data  points  on  the  far 
right  of  the  time  scale  of  Figures  20  through  23.  In  this  temperature  domain 
the  "sign”  of  R  and  1/R  in  Figures  20  through  23  is  indicated  as  "positive”.  An 
explanation  is  presented  in  the  analysis  of  the  experimental  results.  It  is  worth 
stating  that  this  "unstable"  behavior  of  the  sandwich  around  the  glass  transition 
domain  during  cooldown  has  been  observed  in  every  test  run  ;  the  "irregular 
fringe"  pattern  developed  consistently  and  independently  of  the  cooling  rate 
imposed  on  the  specimen. 

4.4  Analysis  of  the  Experimental  Results 

From  the  results  of  the  first  experiments  involving  the  glass-PVA<. 
bimaterial  specimens  it  became  clear  that  it  would  be  necessary  to  define  a 
reference  state  to  be  used  as  a  standard  against  which  the  induced  strains  could 
be  compared  through  the  measurement  of  the  radius  of  curvature.  The  time 
dependent  character  of  the  material  properties  raises  this  question  continu¬ 
ously  ;  this  can  be  shown  very  easily  by  observing  the  change  in  curvature  of  a 
sample  leaving  the  temperature  unchanged.  After  four  hours  at  room  conditions 
(24  °C)  the  radius  of  curvature  R  was  found  to  have  nearly  doubled  its  initial 
value,  clearly  demonstrating  the  relaxational  behavior  of  the  polymer.  From 
these  initial  findings  resulted  the  procedure  followed  during  the  measurements 
leading  to  the  results  shown  in  Figures  20  through  23.  The  plateau  reached  by  R 
at  about  30  °C  above  the  glass  transition  temperature  of  PVA*.  seems  to  meet  the 

8.  The  presence  of  this  irregular  fringe  pattern  is  responsible  for  the  lack  of  data  points  in  the 
corresponding  temperature  domain  of  Figures  20  through  23. 
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requirements  for  the  definition  of  a  reference  state.  As  mentioned  previously 
any  further  increase  in  temperature  does  not  appear  to  change  R  any  more, 
thereby  not  varying  the  existing  stress  levels  anymore.  At  these  temperatures 
the  material  is  approaching  its  viscous  flow-region  where  additional  stresses 
cannot  build  up  any  longer.  Therefore  it  is  postulated  that  this  state  represents 
a  "stress  free"  reference  state  for  the  specimens,  even  though  circular  fringes 
cam  still  be  seen  on  the  TV  monitor.  It  was  suspected  that  these  remaining 
Newton’s  rings  were  caused  mainly  by  the  initial  curvature  of  the  glass  slides. 
The  validity  of  this  assumption  can  be  checked  by  measuring  the  curvature  of 
the  virgin  glass  slides  ;  for  this  purpose  the  polymer  coating  was  removed  from 
its  substrate  by  submerging  the  samples  in  acetone  in  which  polyvinyl  acetate 
dissolves.  The  interference  fringes  generated  in  the  air  gap  between  the  pristine 
glass  slide  and  the  optical  flat  accounted  in  part  for  the  pattern  observed  at  the 
higher  temperatures,  but  also  indicated  that  at  the  above  mentioned  reference 
state  some  residual  curvature  of  the  bimaterial  specimen  was  still  present.  For 
the  thin  glass  slides  the  residual  curvature  between  our  defined  reference  state 
and  the  curvature  of  the  virgin  glass  slide  equals  0.1xl0~3/mm.  In  the  case  of 
the  thicker  slides  the  amplitude  of  this  residual  curvature  could  not  be  deter¬ 
mined  because  the  experimental  error  is  of  the  same  order  as  this  difference 
between  the  reference  state  and  the  curvature  of  a  pristine  glass  slide. 

We  now  proceed  to  make  a  comparison  between  the  experimental 
results  and  the  theoretical  calculations  outlined  in  Chapter  3  ;  the  heating  and 
cooling  thermal  loading  histories  are  analyzed  separately.  Because  it  is  assumed 
that  the  theory  only  accounts  for  small  deformations,  linear  superposition  is 
applied  in  order  to  subtract  the  contribution  of  the  reference  state  from  the 
experimentally  recorded  fringe  patterns.  In  accord  with  the  nomenclature  of 
Figure  17,  the  out-of-plane  displacement  d  assuming  spherical  deflection  is  given 
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by  : 

d  =  R2  -  vRz  -  r2  (39) 

The  thermally  induced  curvature  to  be  compared  with  the  theoretical  results,  is 
calculated  from 


1  _  2dT 

Rt  r2  +  af 


(40) 


where  dr  =  d  -  d<„ 

in  which  dt,  d  and  d<,  represent  the  thermally  induced,  experimentally  recorded 
and  reference  out-of-plane  displacement  at  a  distance  r  from  the  point  of  con¬ 
tact  between  the  sample  and  the  reference  plane  respectively. 

1  In  Figures  24  through  27  experimentally  obtained  data  for  heating  histories 
are  plotted  on  a  comparative  basis  with  the  theoretical  curves,  resulting  from 
the  viscoelastic  stress  analysis.  The  first  two  graphs  refer  to  specimens  with  a 
thin  glass  slide.  The  main  difference  between  theory  and  experiment  lies  in  the 
maximum  amplitude  of  the  change  in  curvature  ;  in  the  four  cases  presented 
here  the  specimens  undergo  a  maximum  change  in  curvature  with  respect  to 
their  reference  state  ranging  from  60%  to  95%  of  the  theoretically  predicted 
value.  In  Figure  27  the  influence  of  the  value  of  the  bulk  modulus  K,,  is  also 
included  ;  two  theoretical  curves  corresponding  to  the  indicated  values  of  K„  are 
drawn.  A  30%  reduction  in  the  amplitude  of  K0  does  not  constitute  a  sizeable 
difference  and  supports  our  assumption  of  a  constant  bulk  modulus  K„. 

Some  characteristics  common  to  the  four  cases  are  worth  mentioning. 
In  every  case  test  data  indicate  that  the  specimen  relaxes  initially,  indicated  by 
a  decline  in  the  value  of  the  curvature  at  an  almost  constant  or  very  slowly 
increasing  temperature.  This  phenomenon  had  been  observed  earlier  when 
observing  the  relaxation  of  a  specimen  at  constant  temperature.  Therefore  the 
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Figure  24.  Change  in  curvature  of  a  bimaterial  specimen  as  a  function  of  tem 
perature  (thin  glass  slide). 


Figure 


o 


-I*02 


ro 


-0.3r 


■0.4 


-0.5h 


- Theory:  K0=8.28xlO  N/m 

v  =0.41 


- Theory  :  K0=  5.5  x  10  N/m‘ 

v=  0.37 


o  Experiment 


55 


50  o 
0  „ 


45  © 
40  o 

L. 

© 

35  a 


30 

25 


E 

.© 


Figure  27.  Change  in  curvature  of  a  bimaterial  specimen  as  a  function  of  tem¬ 
perature  (thick  glass  slide). 
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temperature  profile  used  in  the  theoretical  calculations  has  been  approximated 
toy  a  piece  wise  linear  graph  as  indicated  in  Figures  24  through  27.  In  most  of 
the  cases  the  experimental  data  points  lag  behind  the  corresponding  theoretical 
values,  although  the  peak  value  is  concurrently  reached  somewhere  between 
35  0  and  37  °C.  The  relaxation  occurs  somewhat  faster  than  the  theory  predicts. 
This  is  a  significant  result  of  this  investigation  as  it  appears  that  full  relaxation 
of  the  PVA,.  does  not  start  until  about  7  °C  after  the  glass  transition  has  been 
passed. 

It  will  be  noted  that  some  data  points  indicate  apparent  negative 
values  of  1/R.  This  ’’behavior"  results  from  taking  a  reference  state  equal  to  the 
one  attained  during  the  long  time  relaxation  at  the  higher  temperature.  As  can 
be  seen  in  Figures  20  and  22,  the  radius  of  curvature  corresponding  to  this 
reference  state  is  smaller  in  absolute  value  than  the  ones  reached  immediately 
after  pronounced  relaxation  occurs.  This  phenomenon  is  more  pronounced  in 
the  specimens  containing  the  thinner  glass  slide. 

Two  more  assumptions  made  in  the  theoretical  calculations  need  to  be 
discussed.  They  involve  the  time  dependence  of  the  thermal  coefficient  of  expan¬ 
sion  as  well  as  the  uniformity  of  tue  temperature  of  the  specimen.  In  our 
analysis  no  time  dependence  of  the  expansion  coefficient  is  assumed.  If  one  were 
to  introduce  some  rate  dependence  of  the  thermal  coefficient  of  expansion  a  in 
the  analysis,  the  resulting  effect  would  be  a  shift  of  the  theoretical  curves  in  Fig¬ 
ures  24  through  27  parallel  to  the  (horizontal)  time  or  temperature  axis.  This 
time  dependence  of  a  could  account  for  some  of  the  discrepancies  between 
theory  and  experiment,  but  its  influence  cannot  be  clearly  identified  from  our 
measurements.  We  need  to  point  out  the  measurements  were  performed  at  a 
heating  rate  similar  to  the  one  used  in  the  determination  of  the  thermal  expan¬ 
sion  characteristics  of  PVAe  reported  on  earlier  in  Chapter  2. 
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Concerning  the  uniformity  of  the  temperature  throughout  a  bimateriai 
specimen,  we  need  to  mention  that  the  size  of  the  specimen  is  chosen  such  that 
the  bimaterial  sample  can  respond  "instantaneously"  to  temperature  changes.  It 
turns  out  however  that  the  specimens  are  extremely  sensitive  to  small  tempera¬ 
ture  variations,  which  behavior  results  in  the  previously  mentioned  oscillatory 
behavior  of  the  interference  fringes  in  a  certain  temperature  range.  Neverthe¬ 
less  if  the  temperature  is  kept  constant,  the  fringe  motion  ceases  as  was  experi¬ 
mentally  verified. 

2.  As  indicated  earlier  in  this  chapter  the  cooling  of  the  specimens  from  tem¬ 
peratures  well  above  the  glass  transition  temperature  to  ambient  conditions 
resulted  in  unresolved  experimental  problems,  namely  the  appearance  of  a  com¬ 
plicated.  non-circular  fringe  pattern.  The  repeatability  of  this  phenomenon 
leads  us  to  believe  that  some  kind  of  instability  develops  when  cooling  through 
the  glass  transition  range  of  polyvinyl  acetate.  It  is  believed  that  the  bimateriai 
samples  switch  the  "sign"  of  their  curvature  when  cooled  through  the  glass  tran¬ 
sition  region.  This  phenomenon  was  also  observed  experimentally  when  heating 
bimateriai  samples  from  temperatures  well  below  ambient  conditions  (0  °C)  to 
the  glass  transition  temperature  range. 
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{SION  OF  LINEARLY  VISCOELASTIC  MATERIAL  SPECIMENS  WITH  CIRCULAR  CROSS- 
mONS 

The  relation  between  the  creep  compliance  J(t)  and  the  experimentally 
asured  twist  angle  i?(t)  as  a  function  of  time  is  derived  for  a  linearly  viscoe- 
tic  material.  The  sample  possesses  a  circular  cross-section  ;  the  geometry  is 
>wn  in  Figure  A1  : 


z 


Figure  Al.  Torsion  specimen 

>m  kinematic  considerations  the  shear  strain  is  given  by, 

2e,«(t)  =  £-*(t)  (Al.l) 

i  the  corresponding  shear  stress  is 
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r,4(t)  =  2(jl  *  det«  =  t  (A1.2) 

From  statics  the  twisting  moment  M  can  be  computed  by  : 

d/2  ,  v 

M(t)  =  Znf  Tt«r2dr  =  ^-dV*d(  (A1.3) 

so  that  from  equation  (A1.3)  the  twist  angle  can  be  calculated  as  : 

d(t)  =  ^J(t)  ’dMCt).  (A1.4) 

rrd4 

If  the  moment  M(t)  is  applied  in  a  step  load  fashion,  such  that 

M(t)  =  M0h(t)  (A1.5) 


we  can  write: 


J(t)  4dM(t)  =  /  J(t-0  dM(f) 


=  /  J(t-*)d(K0h(f)) 


=  M0  J(t). 


(A1.8) 


Combining  equations  (A  14)  and  (A1.5)  the  creep  compliance  is  then  given  by 


J(t)  = 


TTd4  d(t) 
32L  M0 


(A1.7) 


which  yields  equation  (1).  It  is  worth  pointing  out  that  this  formula  has  exactly 
the  same  form  as  the  one  used  for  torsion  of  elastic  specimens  with  circular 
cross-sections. 


t  The  star  notation  (•)  refers  to  Stieltjes  convolution  ;  equation  (Al  .2)  should  read  in  full  as  follows 

T  (t)  =  f  ^  M t-$)  d  (  )• 
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APPENDIX2 


STRESS  ANALYSIS  CHECK  FOR  A  BTHATEKIAL  INVOLVING  ELASTIC  HEDIA  ONLY 

If  both  materials  forming  the  bimaterial  specimen  under  consideration 
are  elastic  solids,  the  system  consisting  of  equations  (33),  can  be  solved  expli¬ 
citly  in  terms  of  the  material  and  geometric  parameters. 

The  stress  in  either  material  is  given  by, 


<Tr,(z,t)  =  T—  e0(t)  -  Z”7-(t)  -  04(1(0)  (t) 

1  1  -i/i  dr * 


(A2.1) 


where  i  =  1,2  refers  to  the  range  either  one  of  the  materials.  Substitution  of 

(A2.1)  into  the  equilibrium  equations  and  integrating  over  the  thickness  of  the 

02w 

specimen,  yields  a  system  of  two  equations  and  two  unknowns.  £„(t)  and  — -=-(t), 

dir 

which  can  be  solved  analytically.  The  solution  of  this  system  is  : 

|pr(t)  =  -  — ~-hZ'  [T(t)  -  T0]  [a,(T(t))  -  az(T(t))]  (ht  +  h8)  (A2.2) 


e0(t)  =  [T-(-0— To]-{[^fal(T(t))  h*  +  il  ««(T(t))  hz 


(A2.3) 


-  ^1^2[3h|(hf  +  2h,he)(al(T(t))  -a8(T(t))) 


where 


-  4h,hz(aI(T(t))  h2Z  +  ag(T(t))hf) 


A  =  $jhi  +  $ghg  +  4$i$ehihg(hi  +  hg)+6hihg$i$g 
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1.  INTRODUCTION 

In  this  part  of  the  thesis  we  are  concerned  with  time  dependent  adhesive 
failure  and  long  time  integrity  of  adhesive  systems  ;  emphasis  is  placed  on  the 
adhesion  of  rubbery  materials  to  relatively  rigid  substrates.  Besides  the  viscoe¬ 
lastic  character  of  the  polymeric  coating,  the  amount  of  energy  needed  to  gen¬ 
erate  new  debond  surfaces  needs  to  be  described  quantitatively  When  studying 
adhesive  failure  of  rate  sensitive  materials  one  has  to  recognize  the  influence  of 
the  time  dependent  properties  of  the  polymer  on  the  fracture  energy.  Therefore, 
if  one  attempts  to  determine  this  surface  energy  experimentally,  it  is  necessary 
to  separate  an  intrinsic  and  a  rate  dependent  interface  energy  from  the  meas¬ 
ured  fracture  energy.  The  rate  independent  energy  then  corresponds  to  that 
obtained  for  infinitesimally  slow  interfacial  crack  propagation  rates,  and  conse¬ 
quently  should  not  vary  if  different  test  methods  are  used  to  measure  it.  The 
rate  dependence  needs  to  be  clarified  further.  Let  the  fracture  or  adhesive 
energy  measured  at  some  rate  of  unbonding  c  be  denoted  by  7(6).  Then  as  origi¬ 
nally  formulated  by  Muller  and  Knauss  [  1]  one  may  write 

7  =  r  .  g(c)  (i) 

where  T  is  the  intrinsic  fracture  energy  and  g(c)  is  a  nondimensional  function  of 
c.  In  [l]  the  theory  states  specifically  that  the  fracture  energy  can  be 
represented  by  the  product  of  an  intrinsic  adhesive  failure  energy  I\  which 
depends  solply  on  the  physical  and  chemical  nature  of  the  fracture  interface, 
and  a  function  g(c)  accounting  for  the  rheological  losses  based  on  the  descrip¬ 
tion  governed  by  the  WLF  equation.  Efforts  leading  to  establishing  a  characteris¬ 
tic  measure  of  joint  strength,  independent  of  test  geometry  and  loading  condi¬ 
tions,  were  carried  out  by  Gent  and  coworkers  [2,3,4].  Their  studies  involved 
three  types  of  testpieces,  i.e,  simple  extension,  pure  shear  and  peeling  samples. 
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The  formulation  exemplified  by  equation  (1),  was  Later  confirmed  by  Andrews 
and  Kinloch  [5].  The  details  of  their  experimental  work  can  be  found  in  refer¬ 
ences  [0,7]. 

In  the  determination  of  the  intrinsic  energy  of  adhesion  F  and  the  rate 
dependent  adhesive  energy  y  many  different  test  methods  have  been  used  over 
the  years  and  newer  ones  are  under  study.  The  former  ones  include  blister 
tests,  cone  tests,  peel  tests,  lap  shear  and  cantilever  beam  tests  ;  they  have  been 
extensively  described  in  the  literature  [8].  A  fairly  complete  overview  of  the  vari¬ 
ous  test  methods  for  adhesive  joints  has  been  completed  recently  by  Kinloch  [9]. 
The  choice  of  specific  method  to  measure  the  energy  associated  with  interfa¬ 
cial  fracture,  whether  it  concerns  T  or  y  .  depends  to  a  large  extent  on  the  func¬ 
tion  of  the  adhesive  system. 

Of  particular  interest  in  a  fracture  mechanics  approach  to  adhesive  failure  is 
the  interaction  of  the  local  crack  tip  deformation  modes.  In  the  theory  of  frac¬ 
ture  mechanics  this  notion  refers  to  the  different  types  of  deformations  occur¬ 
ring  in  elastic  media  or  adhesive  systems  with  an  initial  flaw.  In  general.  Mode  I 
refers  to  a  displacement  during  which  the  crack  surfaces  in  the  crack  tip  vicin¬ 
ity  move  directly  apart,  perpendicular  to  the  plane  of  the  crack.  Mode  II  is  asso¬ 
ciated  with  a  displacement  in  which  the  crack  surfaces  slide  over  one  another, 
i.e.  perpendicular  to  the  leading  edge  of  the  crack.  In  Mode  III  a  similar  motion 
occurs  parallel  to  the  leading  edge  of  the  crack.  In  monolithic  materials  dis¬ 
placements  applied  normal  to  and  tangential  to  the  crack  plane  result  only  in 
normal  and  tangential  displacements  of  the  crack  tip  respectively  In  adhesive 
systems  however  combinations  of  various  deformation  modes  appear  in  general 
under  all  loading  conditions.  Their  existence  has  been  demonstrated  experi¬ 
mentally  by  Liechti  [10],  leading  to  the  formulation  of  a  vectorial  crack  opening 
displacement  criterion  in  adhesive  systems  accounting  for  the  presence  of  the 
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two  fracture  modes.  This  observation  raises  the  question  of  the  mode  depen¬ 
dence  of  the  adhesive  energy.  Several  authors  have  reported  a  dependence  of 
the  adhesive  fracture  energy  on  the  loading  mode  [11,12,13]  resulting  from  a 
variety  of  test  methods.  Several  reasons  for  this  apparent  variation  can  be  men¬ 
tioned  however  [12]  ;  among  them  we  mention  the  different  rates  of  viscous  dis¬ 
sipation  and  the  actual  fracture  surface  area  for  each  loading  mode. 

The  viscoelastic  character  of  the  adhesive  material  is  reflected  in  the  rate 
dependence  of  the  adhesive  fracture  energy.  This  characteristic  also  provides  a 
way  of  checking  whether  a  material  is  thermorheologically  simple,  if  the  ampli¬ 
tude  of  the  residual  thermal  stresses  is  negligible.  Under  those  circumstances, 
one  can  perform  debonding  experiments  to  measure  the  adhesive  energy  at 
different  temperatures  necessarily  including  the  glass  transition  range  of  the 
polymer.  If  a  master  curve,  displaying  the  adhesive  energy  versus  the  crack 
speed,  can  be  constructed  from  the  data  by  the  time-temperature  superposition 
principle,  the  material  is  said  to  be  thermorheologically  simple. 

One  needs  to  keep  in  mind,  however,  that  the  temperature  itself  remains  a 
significant  experimental  parameter  because  it  can  be  associated  with  one  of 
several  other  potential  sources  of  energy  available  to  the  debonding  process,  in 
addition  to  the  mechanical  load  being  applied  during  these  tests.  One  energy 
source  arises  from  the  strain  energy  due  to  thermal  stresses  caused  by  the 
mismatch  between  the  coefficient  of  expansion  of  adherend  and  substrate,  if  the 
tests  are  performed  at  temperatures  far  different  from  the  one  at  which  the 
polymer  was  deposited  on  its  substrate.  Furthermore  the  curing  of  the  polymer 
and  the  attendant  shrinkage  causes  straining  and  hence  form  an  additional 
source  of  strain  energy.  The  presence  of  these  residual  stresses  has  also  to  be 
accounted  for  in  a  quantitative  manner. 


Much  attention  has  been  directed  towards  the  determination  of  the  complex 
stress  field  around  the  interface  crack  tip.  While  investigating  some  problems  of 
geophysical  interest,  Williams  [14]  established  in  the  framework  of  the  linearized 
theory  of  elasticity  that  the  singular  behavior  of  the  stress  is  proportional  to  the 
inverse  square  root  of  the  distance  from  the  tip  of  the  crack,  but  that  the 
stresses  possess  a  sharp  oscillatory  character  confined  to  the  immediate  vicinity 
of  the  crack  tip.  Sih  and  Rice  [15]  confirmed  Williams’  finding  of  the  oscillatory 
character  of  the  stresses  ;  these  authors  also  showed  how  the  complex-variable 
method  combined  with  eigenfunction  expansion  used  in  references  [14,15]  can 
be  applied  to  formulate  the  problem  of  bonded  dissimilar  elastic  planes  contain¬ 
ing  cracks  along  the  bond  [l8]  England  [17]  examined  the  problem  of  a  crack 
opened  under  constant  pressure  between  two  dissimilar  materials.  Erdogan  con¬ 
sidered  the  case  of  two  bonded  semi-infinite  planes  containing  cracks  along  the 
bond  [18].  All  these  analyses  possess  the  oscillatory  singularities,  which  require 
the  upper  and  lower  surfaces  of  the  crack  to  overlap  near  the  ends  of  the 
cracks.  In  his  work  Erdogan  concluded  that  in  practical  applications  the 
phenomenon  of  stress  oscillation  may  be  ignored  [18],  However  more  recently 
Knowles  and  Sternberg  [19]  showed  that  the  oscillatory  singularities  arising  in 
mterface-crack  problems  stem  from  the  linearization  of  such  problems,  as 
alluded  to  by  England  [17],  Knowles  and  Sternberg  found  the  crack  to  open 
smoothly  near  its  ends,  regardless  of  the  specific  loading  at  infinity. 

In  their  continuing  effort  to  obtain  a  better  understanding  of  time-dependent 
fracture  in  adhesive  systems,  researchers  have  tried  to  formulate  equations 
governing  fracture  in  viscoelastic  solids.  Williams  [20]  discussed  the  essential 
similarity  between  cohesive  Griffith  fracture  and  adhesive  bond  fracture,  both  of 
which  in  their  elastic  model  possess  singular  stresses  near  the  crack  tip  It  then 
becomes  natural  to  attempt  to  use  similar  fracture  criteria  for  configurations 
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including  adhesive  bonds.  In  [21]  Knauss  formulates  two  fracture  criteria  appli- 
able  for  cohesive  crack  propagation  in  viscoelastic  materials.  One  is  the  crack¬ 
opening  displacement  criterion  which  states  that  continued  fracture  is  possible 
when  the  displacement  u  at  the  trailing  end  of  the  cohesive  zone  attains  some 
constant  and  critical  value  u0.  The  other  criterion  equates  the  work  done  by  the 
boundaries  of  the  cohesive  material  on  the  disintegrating  material  which  is 
equal  to  a  fracture  energy.  In  contrast  to  the  description  of  fracture  in  rate- 
insensitive  materials  it  appears  that  the  size  of  the  damage  or  process  zone  is  a 
necessary  parameter  in  the  discussion  on  fracture  in  (linearly)  viscoelastic 
solids.  For  many  applications  it  is  sufficient  to  assume  that  the  growth  law  for 
steady  crack  propagation  applies  instantaneously.  In  [21]  Knauss  points  out 
that  if  the  stress  intensity  factor  does  not  change  much  during  the  time  interval 
in  which  the  crack  travels  the  length  of  the  cohesive  zone,  the  steady  crack  pro¬ 
pagation  equations  still  apply  In  another  paper  Knauss  [22]  presented  an 
expression  for  the  rate  of  interfacial  unbonding  c,  if  two  poorly  compressible 
viscoelastic  solids  are  joined  ;  for  plane  strain  conditions  the  crack  speed  is 
given  by, 

Kz  [c,l,2]  =  ~~jT  (2) 


D,  (S-)  +  Da  (-) 
c  c 


where 
Dt  and  D? 
a 

c.c 


K  [  c,  1,  2  ] 

The  stress  intensity 


=  the  creep  compliance  of  the  two  joined  solids, 

=  length  of  the  process  zone, 

=  crack  size  and  crack  tip  velocity, 

=  intrinsic  surface  energy  required  for  unit 

crack  extension. 

=  crack  tip  stress  intensity  factor 

factor  K  depends  in  general  on  the  material  properties  of 


the  two  solids  which  are  indicated  by  the  numerals  1  and  2.  Numerical  efforts  to 
evaluate  the  stress  intensity  factors  for  bimaterial  bodies  have  been  carried  out 
by  Hong  and  Stern  [23],  and  Smelser  [24],  In  the  present  study  the  strain  energy 
release  rate  is  chosen  as  the  governing  parameter  for  mterfaciai  crack  propaga¬ 
tion. 

Before  we  can  attempt  to  carry  out  a  failure  analysis  it  follows  from  equation 
(2)  that  some  other  pertinent  data  need  to  be  determined.  Of  special  interest  is 
the  surface  energy  associated  with  the  generation  of  new  surfaces  as  indicated 
earlier.  The  measurement  of  the  fracture  energy  forms  the  subject  of  the  study 
in  Chapter  2  where  the  influence  of  geometrical  and  environmental  conditions  is 
also  investigated.  Subsequently,  in  Chapter  3,  two  delamination  models  only 
involving  elastic  adherends  and  adhesives,  are  developed  based  on  the  strain 
energy  release  rate  criterion  ;  this  formulation  leads  to  crack  propagation  cri¬ 
teria  for  elastic  bimaterial  specimens  subject  to  thermal  and  moisture  loading 
conditions.  The  next  phase  of  this  work  also  includes  the  influence  of  the 
viscoelastic  character  of  the  polymer  layer  in  a  bimaterial  sandwich.  Finally, 
results  of  debond  tests  are  reported  in  order  to  verify  the  validity  of  the  delami- 
nation  criteria.  In  the  experimental  work  a  crosslinked  polyurethane  elastomer, 
Solithane  113,  is  used  as  the  model  adhesive.  Its  properties  have  already  been 
extensively  characterized  in  the  past  [1,25],  and  its  weak  bond  to  glass 
adherends,  ensuring  interfacial  fracture,  makes  it  a  desirable  material  for  these 
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considerably  smaller.  These  considerations  are  represented  qualitatively  in  Fig¬ 
ure  6  The  combination  of  these  two  effects  therefore  shows  a  peak  value  for 
some  intermediary  value  of  the  thickness  of  polymer  films.  Numerical  values  for 
the  dissipation  D  calculated  for  the  different  film  thicknesses  in  accord  with 
equation  (10)  show  a  behavior  similar  to  the  one  exhibited  in  the  inset  of  Figure 
4  m  a  qualitative  manner. 

An  important  remaining  question,  and  unanswered  in  this  study,  is  the  defor¬ 
mation  mode  interaction  in  the  peel  test  One  can  expect  the  peel  angle  to  play  a 
role  in  determining  the  amount  of  Mode  11  loading  at  the  crack  tip.  Anderson 
and  coworkers  contend  [12]  that  in  the  90°  peel  test,  the  major  portion  of  the 
load  should  be  in  Mode  I.  However,  it  is  reasonable  to  assume  that  some  Mode  II 
deformation  is  present,  due  to  bending  near  the  debond  tip.  In  more  recent 
work  Crocombe  and  Adams  [13]  used  the  finite  element  method  to  show  that  the 
proportion  of  Mode  II  loading  is  essentially  independent  of  the  peel  load  and 
angle,  and  the  adhesive  or  adherend  modulus.  The  proportion  appears 
significant,  being  about  30  %  for  compressible  solids  ( v  =  0  3)  and  decreasing  to 
about  .2  Z  as  the  adhesive  becomes  incompressible.  The  results  in  the  litera¬ 
ture  on  this  matter  are  still  very  limited.  One  can  however  conclude  that  an 
appreciable  amount  of  Mode  II  deformation  at  the  crack  tip  is  present  in  the 
peel  test.4 

2  4  2  Temperature  Effect  in  Fracture  Energy  Determination  The  temperature 
is  expected  to  have  two  distinct  effects  on  the  debonding  process 

1  First  it  influences  the  thermorheological  behavior  of  the  polymer 

2  The  temperature  introduces  thermal  stresses  and  corresponding  elastically 


4 


4  was  considered  to  investigate  the  interaction  of  the  deformation  modes  in  our  test  geometry 
inrough  a  finite  element  code,  containing  a  bimateriei  cracx  element.  The  TFXGA?  computer 
code  not  uemg  -eadiiy  avauaole  we  were  unable  to  include  these  resuits  m  our  study. 
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ment  with  other  investigations  [ 27,32]  at  least  at  first  sight.  In  these  referenced 
studies  an  increase  of  y  with  thickness  is  reported,  a  conclusion  that  cannot  be 
drawn  from  our  data.  Especially  the  results  for  the  two  thickest  films  need 
further  explanation. 

When  plotting  values  of  7  as  a  function  of  the  thickness  h  for  a  particular 
crack  speed  c,  as  indicated  in  the  inset  in  Figure  4,  one  observes  a  peak  for 
thickness  value  of  h  =  1.6  mm  (1/16").  Let  us  attempt  to  explain  this  observa¬ 
tion  in  terms  of  certain  processes  not  exphcitely  accounted  for  in  equation  (3). 
The  influence  of  the  high  stresses  and  strains  as  well  as  of  the  viscous  dissipa¬ 
tion  around  the  moving  crack  tip  is  hard  to  determine  directly.  One  can  try. 
however,  to  estimate  the  bulk  viscous  dissipation  D  in  the  polymer  sheet  There¬ 
fore  we  determine  first  the  bending  strain  t  and  the  associated  strain  rate  t  in 
the  plate  as 

*  -  fr  :b) 


and 


c  y 


ds 


(9) 


where  y  denotes  a  point  location  from  the  midpiane  of  the  polymer.  The  bulk 
viscous  dissipation  D  is  proportional  to  i2  and  the  material  volume  V,  according 
to 

D  ~  s2V  -  c 2f  ^  i(d..y^gl)z  ds  (10) 

J0  ds 

In  this  representation  two  counteracting  effects  are  present  In  the  thinner 
films,  which  always  form  a  full  90°  bend,  the  strains  and  thus  the  strain  rates 
are  higher  ;  in  the  thicker  Alms,  which  never  display  a  complete  90°  bend  in 
these  tests,  the  volume  of  the  material  remains  higher,  but  the  strain  rates  are 
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Figure  5.  Polynomial  (of  degree  2)  least  squares  fit  through  the  experimental 
data  of  the  peel  tests. 
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restriction  on  the  accuracy  of  these  measurements  derives  from  the  limited 
resolution  with  which  the  debond  position  could  be  measured.  This  problem 
results  from  the  need  to  monitor  the  tests  from  outside  of  the  thermal  chamber 
in  which  the  experiments  were  performed.  The  resolution  of  the  ruler  indicating 
the  crack  propagation  is  around  0.5  mm.  The  accuracy  of  the  various  geometri¬ 
cal  parameters  is  also  subject  to  limitations,  as  they  have  to  be  obtained  from 
photographic  recordings.  Figure  5  shows  the  results  of  a  polynomial  (of  degree 
2)  least  squares  fit  through  the  experimental  data. 

From  Figures  4  and  5  it  is  apparent  that  the  energy  consumption  in  the 
separation  process  is  a  function  of  the  adhesive  thickness  As  pointed  out  in  the 
Introduction  to  Part  II,  one  of  the  main  points  of  interest  remains  the  limiting 
value  of  the  adhesive  fracture  energy  for  vanishing  crack  speeds.  The  data 
recorded  for  the  four  different  thicknesses  seem  to  converge  to  a  common  lower 
limit  of  about  f  =  .02  g/mm  ±  30%,  which  one  identifies  as  the  intrinsic  energy 
of  adhesion  F  ;  for  the  limit  value  F  of  the  adhesive  energy  the  thickness  effect 
seems  to  disappear. 

We  now  discuss  the  reasons  for  the  increase  in  the  adhesive  energy  7  as  the 
crack  propagation  increases.  The  deviation  from  the  limiting  value  of  the  intrin¬ 
sic  energy  of  adhesion  F  is  due  to  a  number  a  processes  which  are  not  included 
in  the  energy  balance  equation  (3).  High  stresses  occur  at  the  tip  of  the  inter¬ 
face  crack  ;  in  addition  to  the  strain  energy  associated  witn  these  stresses, 
viscous  dissipation  occurs  in  the  viscoelastic  material  in  the  vicinity  of  the  mov¬ 
ing  crack  tip.  This  dissipation  is  clearly  a  function  of  the  crack  tip  local  strain 
rates  which  in  turn  are  proportional  to  the  crack  speed  c  An  additional  energy 
absorbing  process  is  the  bulk  viscous  dissipation  in  the  polymeric  adhesive. 


Even  with  these  considerations  in  mind,  these  results  are  not  in  full  agree- 
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Figure  4.  Adhesive  energy  7  versus  crack  propagation  rate  c  for  various  film 
thicknesses  at  ambient  temperature. 


check  whether  the  solution  y(x)  of  equation  (4)  generates  the  experimentally 
recorded  deflection. 

In  this  comparison  it  is  tried  to  match  the  experimental  and  calculated 
deflection  and  slope  of  the  free  end.  Using  the  boundary  conditions  (5)  does  not 
yield  the  experimentally  recorded  deflection  curve.  Therefore  the  boundary  con¬ 
dition  describing  the  slope  at  the  assumed  "built-in”  end  (x  =  0)  needs  some 
further  clarification.  The  existing  moment  causes  a  rigid  body  like  rotation  of 
the  debonded  polymer  film.  Thus  the  second  boundary  condition  in  equation  (5) 
is  changed,  allowing  small  initial  angles  varying  (from  0°  to  7°)  at  the  clamped 
end  ;  this  correction  can  be  thought  of  as  aui  elastic  rotation  correction  at  the 
built-in  end.  A  similar  consideration  was  formulated  by  Ripling  and  Mostovoy 
[31]. 

From  equation  (3)  the  adhesive  peel  energy  y  can  be  calculated  as  a  function 
of  the  crack  propagation  speed  c.  During  the  tests  with  the  thicker  adherends 
when  a  constant  90°  bend  angle  did  not  result,  the  effective  loading  varied  and 
the  debond  speed  seldom  stayed  constant,  due  to  continuous  variations  in  the 
geometry  and  thus  the  moment,  hence  the  parameters  L  and  l  were  also  slowly 
varying  functions  of  time. 

2.4  Experimental  Results  and  Discussion 

2-4.1  Dependence  of  the  Fracture  Energy  on  Film.  Thickness  Initial  tests  were 
performed  at  room  temperature,  using  weights  ranging  from  1  g  to  15  g,  depend¬ 
ing  on  the  film  thickness.  The  four  different  "film''  thicknesses  employed  were 
0.8  mm  (1/32"),  1.6  mm  (1/16”),  3.2  mm  (1/8”)  and  6.4  mm  (1/4").  The  results 
are  shown  in  Figure  4  displaying  the  peel  energy  y  per  unit  crack  extension 
versus  the  crack  speed  c.  The  data  exhibit  distinct  scatter,  but  that  scatter  is 
considerably  smaller  than  that  in  many  data  reported  in  the  literature  [4].  One 


soning  the  following  relation  results  . 
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where  F  and  y  are  the  peel  force  and  the  adhesive  energy  respectively  . 
v,  a.  U.  y.E.  1  and  b  indicate  the  downward  velocity  of  the  free  end  of  the  cantil¬ 
ever  plate,  the  crack  speed,  the  length  of  the  film  (along  its  midplane),  Poisson's 
ratio.  Young's  modulus,  the  section  moment  of  inertia,  and  the  width  respec¬ 
tively  ;  M  is  the  moment  resulting  from  the  applied  peel  force  and  the  body 
forces.  The  energy  associated  with  the  longitudinal  extension  of  the  polymer 
was  found  to  be  negligible. 


In  order  to  calculate  the  work  done  by  the  body  forces  as  well  as  its  bending 
strain  energy,  the  deflection  shape  of  the  cantilever  plate  needs  to  be  known.  In 
accordance  with  the  theory  of  large  deflections,  the  following  differential  equa¬ 
tion  governs  the  deflection  y(x)  (R  =  radius  of  curvature) 
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subject  to  the  boundary  conditions  : 

y(o)  =  0 
y'(o)  =  o. 

Furthermore, 


(5) 


M  =  F  (L  -  x)  +  ^7 - for  0  <  x  <  l 

2cosp 

M  =  F  (L  -  x)  for  1  <  x  <  L 


(6) 

(7) 


in  agreement  with  the  nomenclature  of  Figure  3  ;  q  denotes  the  weight  of  the 
polymer  sheet  per  unit  length.  The  results  of  the  numerical  calculations  are 
compared  with  the  data  compiled  from  the  photographic  records,  in  order  to 
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At  higher  temperatures  the  dissolved  gas  in  the  saline  solution  formed  bub¬ 
bles  on  the  specimen  thus  affecting  the  net  weight  on  the  specimen.  In  order  to 
prevent  the  formation  of  air  bubbles  a  0.5?!  volume  of  Kodak  Photo-Flo  200  solu¬ 
tion  ls  added  to  the  fluid.  In  addition  the  latter  was  boiled  before  the  actual  test¬ 
ing  to  drive  off  the  dissolved  air. 

2.3  Analytical  Developments 

As  indicated  in  the  Introduction  to  this  chapter,  one  of  the  primary  goals  of 
this  investigation  remains  the  influence  of  the  thickness  of  the  adhesive  on  the 
adhesive  energy  y.  In  addition  we  also  want  to  include  the  viscoelastic  character 
of  the  polymer  film  in  our  study. 

When  performing  peel  tests  it  is  often  desirable  to  obtain  a  90°  bend  in  the 
debonding  adhesive  film  ;  if  the  geometry  remains  constant  during  the  time  of 
the  experiment,  it  makes  the  analysis  necessary  to  evaluate  the  adhesive  frac¬ 
ture  energy  relatively  straightforward.  For  the  range  of  thicknesses  involved  in 
our  experiments  it  was  found  that  the  thicker  films  could  not  always  form  a  full 
90°  bend  ;  this  situation  occurred  for  the  range  of  peel  weights  used  if  the  thick¬ 
ness  exceeded  2  mm. 

Next  we  proceed  to  formulate  an  energy  balance  to  determine  the  fracture 
energy  y.  The  geometry  of  the  peel  test  is  shown  in  Figure  3.  For  analysis  pur¬ 
poses  we  assume  that  the  glass  plate  remains  mechanically  rigid.  Work  is  being 
done  by  the  descending  peel  force  as  well  as  by  the  body  forces  of  the  hanging 
portion  of  debonded  polymer.  The  latter  part  includes  work  as  a  result  of  the 
translation  of  the  "film"  as  well  as  of  the  rotation  of  the  free  hanging  polymer  in 
the  case  of  thicker  films  when  the  peel  angle  does  not  remain  constant.  The 
energy  thus  generated  serves  to  create  new  surfaces  and  is  also  stored  as  bend¬ 
ing  strain  energy  in  the  debonded  film.  Hence,  using  "strength  of  materials"  rea- 
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a  prescribed  load  orientation  to  the  adhesive  and  aimed  at  gradually  removing  it 
from  its  adherend  with  a  more  or  less  constant  rate.  This  can  be  achieved  by 
either  prescribing  a  constant  displacement  rate  or  by  applying  a  constant  load 
to  the  adhesive  joint.  In  the  former  configuration  the  peel  testing  apparatus 
needs  to  include  a  system  that  allows  continuous  adjustment  of  the  peel  angle  in 
order  to  be  performed  in  a  tensile  testing  apparatus  [13,28].  In  the  present  work 
preference  is  given  to  the  second  way  of  loading,  i.e.  a  constant  weight  is  applied. 
The  experimental  set-up  in  the  latter  case  also  is  simpler  to  assemble.  In  addi¬ 
tion  the  need  to  enclose  the  testing  apparatus  in  an  environmental  chamber 
allowing  control  of  the  test  temperature  pointed  towards  the  use  of  a  simple 
method. 

In  addition  to  the  dead  weight,  we  have  to  account  for  the  increasing  weight 
of  the  debonded  polymer  film  as  peeling  progresses.  The  large  thickness  of  the 
Solithane  113  films  in  these  measurements  (from  0.8  mm  to  6.4  mm)  gives  rise 
to  body  forces  on  the  order  of  the  applied  peel  force.  Therefore  it  was  decided  to 
submerge  the  debonded  portion  of  the  adhesive  in  a  fluid  with  a  density  that 
matches  that  of  Solithane  113.  A  5%  by  weight  saline  solution  [30]  satisfies  this 
requirement.3  The  peel  testing  apparatus  is  shown  schematically  in  Figure  2.  It 
consists  of  a  stainless  steel  frame  which  supports  the  specimen  and  rests  on  the 
edges  of  a  water  tank.  The  latter  is  made  out  of  glass  so  that  the  test  geometry 
can  be  monitored  during  peeling. 

The  motion  of  the  crack  tip  is  monitored  via  a  mirror  and  using  a  tran¬ 
sparent  ruler  mounted  on  the  glass  adherend.  At  regular  time  intervals  pictures 
of  the  peel  test  geometry  are  taken  in  order  to  keep  track  of  the  debonding  pro¬ 
cess.  Kodak  technical  pan  film  2415  provided  high  contrast  records, 

3.  The  potential  effect  of  salt  deposition  at  the  interface  on  the  values  of  the  adhesive  energy 

remains  an  unknown  and  is  unaccounted  for  in  this  work. 
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weight  is  thoroughly  degassed  and  mixed.  The  liquid  polymer  is  then  pressure 
injected  into  the  gap  between  the  center  plate  of  the  mold  and  the  glass 
adherends,  through  feeding  channels  drilled  in  the  center  plate.  The  mold  gap  is 
then  sealed  and  the  polymer  cured  at  room  temperature  for  one  week  ;  the  peel 
specimens  are  used  within  a  month  after  the  cure  has  been  completed.  The 
symmetric  design  of  the  fixture  permits  the  simultaneous  casting  of  two  speci¬ 
mens  of  equal  thickness.  This  procedure  reduces  somewhat  the  possibility  of 
differences  in  material  adhesive  properties  between  specimens.  In  order  to 
prevent  sticking  of  the  Solithane  113  to  the  metal  parts  of  the  mold,  these  were 
sprayed  with  a  thin  layer  of  KRAXO  1711  Release  Agent  (Contour  Chemical  Com¬ 
pany).  The  Solithane  113  specimens  are  approximately  290  mm  long,  50  8  mm 
wide,  the  thickness  varying  according  to  the  spacer  used,  bonded  to  the  glass  , 
the  dimensions  of  the  latter  exceed  those  of  the  Solithane  113  by  about  1  cm  all 
around. 

The  cure  at  room  temperature  eliminates  the  residual  stresses,  resulting 
from  the  different  thermal  expansion  properties  of  adherend  and  adhesive. 
However,  the  cure  shrinkage  of  the  latter  introduces  residual  stresses.  Their 
magnitude  has  been  estimated  by  Liechti  [10]  ;  the  cure  shrinkage  was  deter¬ 
mined  to  be  AV /V  =  1  34xl0_z  which  corresponds  to  a  thermal  cool  down  of 
25  °C.Z 

2.2  Experimental  Method  for  Peel  Testing 

The  adhesive  peel  test  can  be  performed  in  different  ways,  including  stripping, 
floating  roller  and  climbing  drum  tests.  All  these  methods  are  designed  to  apply 


1.  The  terms  "Resin"  and  "Catalyst"  are  commercial  designations.  Chemically  the  "Resin"  is  a 
trifunctional  isocyanate  which  is  the  product  of  a  reaction  between  Castor  Oil  and 
Toiyienediisocyanate.  The  "Catalysf’is  a  triol  and  consists  essentially  of  Castor  Oil. 

2.  The  thermal  expansion  characteristics  of  Solithane  1 13  were  determined  previously  m  the 
dilatometer  described  in  Chapter  2  of  Part  I.  Above  the  giass  transition  the  linear  coefficient  of 
thermal  expansion  equals  1  79xlO-4/°C 


Figure  1.  Peel  test  specimen  and  casting  fixture. 


2.  MEASUREMENT  OF  THE  ADHESIVE  FRACTURE  ENERGY 

In  this  chapter  the  determination  of  the  adhesive  energy  is  presented.  The 
peel  method  has  been  widely  used,  because  it  is  an  easily  applied  techmque 
[4,7.13.26,28,29].  However,  as  will  be  seen  below,  the  peel  test  is  not  free  of 
interpretational  difficulties.  It  is  in  part  for  this  reason  that  we  investigate  the 
interfacial  failure  process  in  terms  of  the  peel  test.  The  specific  goals  of  the  peel 
tests  are  : 

1.  To  determine  a  value  for  the  intrinsic  energy  of  adhesion  f,  as  well  as  to 
study  the  rate  dependence  of  the  fracture  energy  y  at  ambient  temperature. 

2.  To  examine  the  effect  of  the  thickness  of  the  polymer  as  well  as  of  its  viscoe¬ 
lastic  nature  on  the  value  of  the  fracture  energy  y  The  latter  influence 
manifests  itself  through  relaxational  behavior  and  viscous  dissipation 
phenomena. 

3.  To  investigate  the  influence  of  residual  (cure)  and  thermal  stresses  on  the 
basis  of  the  thermorheological  behavior  of  the  polymer. 

The  next  two  sections  describe  the  materials  used  and  the  experimental  method. 
Section  2.3  presents  the  energy  balance  considerations  leading  to  the  value  of  y 
The  experimental  results  and  the  discussion  form  the  object  of  Section  2.4. 

2.1  Test  Materials  and  Specimen  Preparation 

In  order  to  allow  visual  monitoring  of  the  peel  test,  glass  was  chosen  as  the 
transparent  adherend.  In  the  Introduction  the  reasons  for  choosing  Solithane 
113  as  the  model  adhesive  have  already  been  given. 

A  specimen  and  the  casting  mold  are  shown  in  Figure  1.  By  means  of  set¬ 
screws  the  glass  plate  is  forced  against  a  spacer  which  determines  the  thickness 
of  the  polymer  film.  A  mixture  of  50%  Solithane  '  Resin"  and  50%  "Catalyst"  1  by 


Figure  6.  Qualitative  representation  of  the  volume  V  of  the  polymer  film,  the 
strain  rate  t,  and  the  dissipation  D  as  a  function  of  the  film  thickness 
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stored  energy. 

In  all  testing  data  available  in  the  literature  at  this  time  only  the  first  of  these 
effects  is  considered,  and  the  effect  of  the  second,  the  temperature  influence  is 
not  considered.  Therefore,  it  is  appropriate  to  examine  whether  this  second 
effect  can  be  neglected,  and  if  so,  under  what  conditions. 

The  enclosure  of  the  set-up  in  the  thermal  chamber  allowed  measurements  in 
the  temperature  interval  ranging  from  0  °C  to  70  °C.  The  lower  limit  was 
imposed  by  the  freezing  point  of  the  saline  solution  (around  -3  °C)  [30]  ;  at  tem¬ 
peratures  above  70  °C  excessive  evaporation  of  the  solution  made  it  more 
difficult  to  keep  the  fluid  level  constant  over  protracted  periods  of  time.  In  order 
to  keep  the  crack  propagation  speeds  at  the  low  end  of  the  velocity  domain,  test 
periods  of  12  hours  were  common. 

The  test  results  for  the  0.8  mm  (1/32”)  and  16  mm  (1/16")  thick  Aims,  are 
displayed  in  Figure  7.  These  data  show  the  interfacial  crack  propagation  rate  c 
as  a  function  of  the  test  temperature,  the  arrows  pointing  downward  in  Figure  7 
represent  zero  crack  speed  data  at  the  corresponding  temperatures.  The  peel 
tests  were  performed  using  the  same  constant  weight  at  all  temperatures.  The 
measurements  were  limited  to  the  thinner  films,  which  keep  their  geometry 
essentially  constant  during  the  peeling  process  In  both  cases  the  crack  speed 
seems  to  reach  a  maximum  value  around  10  0  to  20  °C  above  room  temperature 
At  higher  temperatures  the  crack  speed  becomes  vanishingly  small.  This  latter 
observation  is  in  direct  contradiction  to  the  thermorheological  effect 

For  the  purpose  of  an  explanation  we  note  that  equation  (2)  can  be  written 
as 


Figure  7.  Crack  propagation  rate  c  as  a  function  of  the  test  temperature. 
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where  E„  is  Young  moduies  for  Solithane  1 13  in  its  rubbery  domain  and  G  is  the 
energy  strain  release  rate.  The  contribution  of  the  glass  adherend  on  the  left 
hand  side  of  equation  (2)  has  been  neglected.  In  this  analysis  the  sources  of 
energy  driving  the  crack  propagation  are  twofold.  One  arises  from  the  strain 
energy  stored  in  the  polymer  as  a  result  of  the  cure  shrinkage.  The  strain 
energy  release  rate  due  to  curing  amounts  to  : 


Gc 


Eh  ,  AV  ,g  To  ~T 
1  -  i/  1  3V  '  Tc 


(12) 


where 

Av 

v 


cure  shrinkage  (numerical  values  from  ref.  [10]) 


Tc  =  cure  temperature  (room  temperature)  taken  equal  to  25  °C 

T0  =  Tc  +  25  °C.  T0  represent  temperature  at  which 

the  curing  induced  strains  would  vanish  as  determined 
by  Liechti  [10]. 

For  reasons  of  simplicity  the  strain  energy  release  rate  due  to  the  peel  load,  is 
taken  to  be  : 


Gp  = 


F 

b 


(13) 


Introducing  equations  (12)  and  (13)  into  equation  (11)  yields  : 

2 


D  ( — )  = 
c  ax 


Eh.  1  f  *vc  >2  .  _F_ 
r  i  - 1/ v  3v  '  r  i 


(14) 


From  equation  (14)  the  crack  velocity  c  can  be  calculated  as  a  function  of  the 
temperature  T  with  the  intrinsic  energy  of  adhesion  T  as  a  parameter  ;  the 
results  of  such  a  computation  are  presented  in  Figure  6.  For  lower  values  of  T, 
i.e.  up  to  F  =2  0.020  g/mm.  a  continuous  relation  is  obtained  over  the  complete 
temperature  interval.  For  higher  values  of  T,  equation  (14)  is  no  longer  satisfied 
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over  the  whole  temperature  domain  and  hence  the  criterion  for  crack  propaga¬ 
tion  is  no  longer  met  for  all  values  of  T.  In  order  to  compare  experimental  and 
theoretical  results,  one  tries  to  match  the  results  of  Figure  7  with  one  of  the 
theoretical  curves  of  Figure  8  for  a  constant  value  of  F.  The  theoretical  curve 
corresponding  to  T  =  0.02  g/mm  displays  a  behavior  very  similar  to  the  one 
exhibited  by  our  experimental  data  as  exemplified  in  Figure  9  for  a  film  thick¬ 
ness  h  =  0.8  mm.  It  should  be  borne  in  mind  that  the  cure  shrinkage  strain 
energy  vanishes  around  50  °C,  thus  reflecting  the  correspondence  of  the  cure 
shrinkage  to  a  thermal  cooldown  of  25  °C  [10],  At  temperatures  removed  from 
50  °C.  the  strain  energy  stored  in  the  material  volume  as  a  result  of  curing 
becomes  much  more  important.  From  our  experimental  results  it  appears  that 
this  source  of  strain  energy  plays  a  significant  role  in  the  lower  end  of  the  test 
temperature  range.  We  can  also  conclude  that  the  value  of  the  intrinsic  frac¬ 
ture  energy  is  somewhat  less  than  0  02  g/mm. 

Two  more  considerations  concerning  equation  (14)  deserve  attention.  The 
first  one  concerns  the  value  of  F  to  be  substituted  in  equation  (14).  Besides  the 
peel  weight,  the  body  forces  connected  with  the  unbonded  polymer  not  yet  sub¬ 
merged  in  the  saline  solution,  also  pull  on  the  still  adhering  portion  of  the  Sol- 
ithane  113  sheet.  Therefore  estimates  were  made  to  account  for  the  body  forces 
in  some  sample  calculations  ;  no  significant  change  from  the  previous  calcula¬ 
tions  was  obtained  however. 

Another  point  of  interest  is  the  length  of  the  damage  zone  a.  In  order  to 
effect  a  fit  of  the  analysis  represented  in  Figure  8  with  the  experimental  data  of 
Figure  7  a  value  of  a  =  2k  was  necessary.  In  connection  with  Figure  8  it  should 
be  pointed  out  however,  that  small  variations  in  T  result  in  rather  significant 
changes  in  the  value  of  c.  Hence  variations  in  F,  corresponding  to  vertical  shift¬ 
ing  in  Figure  8,  can  easily  raise  the  value  of  a  up  to  lOoA.  For  comparison  the 
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9.  Comparison  between  calculated  and  experimentally  obtained  crack 
propagation  rate  as  a  function  of  temperature. 


value  of  a  in  cohesive  fracture  in  Solithane  113  obtained  by  Knauss  and  Muller 
can  be  referred  to,  i.e.  a  =  150A  [  1].  Hence  our  theoretical  results  seem  to  con¬ 
stitute  a  reasonable  source  for  comparison  with  the  experimental  data,  and 
underline  the  importance  of  the  residual  (cure)  stresses 


3.  DELAHINATION  MODELS 

In  this  section  two  examples  of  failure  progression  are  considered  ;  they 
exemplify  the  fracture  mechanics  approach  to  two  typical  adhesive  failure 
processes  germane  to  protective  film  debonding.  An  energy  criterion  is  applied 
in  order  to  predict  crack  propagation,  for  the  case  that  initial  imperfections 
(debond)  are  present  at  the  polymer-substrate  interfaces. 

In  one  example  case  shown  in  Figure  10a,  one  is  concerned  with  the  growth  of 
a  debond  at  the  edge  of  a  interface  between  a  polymer  and  a  rigid  substrate. 
Straining,  induced  by,  say,  a  decreasing  temperature  or  by  outgassing  of  the 
polymer,  initiates  and  drives  the  crack  propagation.  The  second  example  out¬ 
lined  in  Section  3.2,  involves  the  buckling  of  the  adhesive  layer,  after  a  portion 
of  the  coating  has  separated  from  the  substrate.  Subsequent  heating  or  solvent 
water  infusion  into  the  polymer  coating  of  the  bimaterial  structure  may  cause 
the  debond  to  grow. 

In  both  cases  the  circular  geometry  of  the  specimen  favors  axisymmetry  of 
the  problem.  It  is  also  appropriate  to  summarize  the  assumptions  made  in  the 
sample  calculations,  namely  : 

1.  The  coating  material  behaves  like  a  linearly  elastic  material  except  possibly 
in  the  highly  stressed  region  at  the  crack  tip.  This  assumption  describes 
well  the  general  behavior  of  a  polymer  far  away  from  its  glass  transition. 
Viscoelastic  effects  will  be  discussed  later  on. 

2.  The  adhering  layer  possesses  a  uniform  temperature  throughout  its  volume. 

In  each  example  the  elastic  strain  energy  release  rate  G  is  calculated  which 
equals  to  the  derivative  of  the  elastic  energy  with  respect  to  the  crack  surface 
increase.  The  crack  propagation  criterion  used  in  order  to  predict  the  onset  of 
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Figure  10.  Delamination  configurations. 
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additional  debonding,  was  originally  formulated  by  Griffith  ;  it  states  that  crack 
propagation  will  occur  if  the  energy  released  upon  crack  growth  is  sufficient  to 
provide  all  the  energy  required  for  crack  growth.  For  an  axisymmetric  geometry 
the  criterion  can  be  written  as, 


£  _  9U_  _  i  au 

3A  2ttR  3R 


C-5) 


where  dA  is  the  area  to  be  debonded,  while  R  and  F  is  the  radial  coordinate  and 
the  energy  required  to  create  a  new  surface,  respectively.  The  study  of  the  edge 
delamination  is  complemented  with  experimental  results  of  debond  tests,  del¬ 
ineated  in  Section  3  4 


3. 1  Edge  Delamination 

The  geometry  of  the  edge  delamination  case  is  shown  in  Figure  1 1  A  circular 
polymer  film  of  thickness  h  and  radius  R0  is  bonded  to  a  rigid  substrate  over  a 
circle  of  radius  R  so  that  an  annulus  R0-R  is  unbonded,  where  R>>  Rc-R  and 
R0-R>>  h  The  latter  assumption  allows  neglecting  the  more  complex  stress 
analysis  necessary  if  the  size  of  the  debond  were  comparable  to  the  film  thick¬ 
ness  Upon  changing  the  temperature,  the  film  tends  to  contract  or  expand,  but 
because  of  the  partial  bond,  an  equibiaxiai  homogeneous  stress  field  is  esta¬ 
blished  inside  a  circle  of  radius  rc  ,  where  rc  >  >  h.  We  aim  at  calculating  the 
temperature  T  at  which  unbonding  will  occur  if  T0  is  the  temperature  at  which 
bonding  was  effected 

In  the  circular  region,  where  r  <  rc,  one  finds  the  stress  field 


aT  =  a*  =  a 


Due  to  adhesion  to  the  rigid  substrate,  the  in-plane  mechanical  strains  vanish, 


namely, 
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Figure  11.  Geometry  of  the  edge  delamination. 
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and  hence, 


£r  —  £4  —  Q 


-E  eT 

a  =  - 

1  -v 


where  *>r  =  a  (T-T0  )  ;  E  and  v  are  Young’s  modulus  and  Poisson's  ratio  respec¬ 
tively  :  A  a  is  the  difference  between  the  coefficients  of  thermal  expansion  of  the 
polymer  and  the  substrate  ;  the  value  of  the  latter  can  be  neglected  with  respect 
to  the  one  for  the  polymer,  and  hence  A  a  -  a.  The  equibiaxial  character  of  the 
stress  field  yields  the  strain  energy  per  unit  volume 


u  =  crz  (1  -  v) 


so  that  the  total  stored  energy  U  is  given  by 

~  n  rf  h( !-(/)-  ~~  ttR2  h(  1  -  v)  ■  (20) 

In  order  to  create  a  new  annular  surface  Aa  (of  radius  R  with  width  AR)  as  the 
debonding  continues,  the  energy  AS  required  is 

AS  =  2tvR  AR  7  .  (21) 

During  the  unbonding  process,  energy  stored  elastically  In  the  sheet  is  released. 
The  rate  at  which  this  occurs  per  unit  bond  radius  equals 

W  '  aR=  —  "Rh(1  -“l  :22) 

Upon  substituting  equation  (18),  one  obtai  -s  : 

G  =  _u.|t  =  Eh^Ini)l  .») 


According  to  equation  (15)  unbonding  can  occur  when 
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T  -  Tg  >  ^ 
a 


ZiLzA  i  y/z 

Eh 


(24) 


Note  that  the  strain  energy  release  rate  G  is  independent  of  the  bond  size  R, 
namely 

G  =—-(!-«/)•  (25) 


As  already  suggested  earlier,  a  more  detailed  analysis  is  needed  for  small  values 
of  (R0-R)/h  as  indicated  in  Figure  12.  Numerical  methods  seem  more  appropri¬ 
ate  to  solve  problems  relating  to  the  early  stages  of  thts  particular  kind  of  crack 
growth,  namely  crack  initiation  and/or  small  initial  flaws.  In  the  region  where 
R/h  >  R  */h,  debonding  will  proceed  continuously,  and  presumably,  at  a  con¬ 
stant  rate  if  enough  energy  is  being  provided  to  create  new  crack  surfaces,  i  e  if 
G  >  7  If  G  <  7,  no  separation  will  take  place. 


Finally  one  notes  that  strain  energy  release  rale  G  may  result  from  causes 
other  than  temperature  induced  straining,  namely  solvent  (water)  swelling  and 
mass  loss  (outgassing  migration  of  low  molecular  weight  components).  Then  the 
non-dimensionalized  strain  energy  release  rate  is  given  by 


as  long  as  e0  is  the  strain  due  to  thermal  and/or  moisture  loading  and/or  due  to 
outgassing  of  the  polymer,  i.e. 


c0  -  A  a  (T  -  T0)  +  otHzo  C  -  0 


[27) 


where  an2o.  C,  f?  are  the  coefficient  of  expansion  due  to  water  (or  other  solvent) 
absorption,  the  water  (or  solvent)  concentration  and  the  shrinkage  due  to  out- 
gassing  respectively 
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we  account  for  the  cure  shrinkage  stresses,  the  strain  energy  release  rate  G  can 
be  written  as  : 

G  =  a  (T0-T  +  ATC)  2  (37) 

where  ATC  =  25  °C  represents  a  thermal  cooldown  equivalent  to  the  cure  shrink¬ 
age  [  10]. 

According  to  equation  (37)  debonding  could  start  for  a  AT  =  35  °C  and 
AT  =  24  °C  for  the  two  thicknesses  involved  in  the  debond  tests.  From  equation 
(37)  one  concludes  that 

G  ~  h  (T0  -T  +  ATC  )  2  (38a) 

~  h  (AT)2  (3Bb) 

When  evaluating  equation  (38b)  the  ratio  of  the  average  strain  energy  release 
rate  G  associated  with  each  polymer  film  thickness,  is  equal  to 


G. 

Ga 


h,(ATi)2 

h2(ATa)2 


0  8(1 00)2 
1  6(90)2 


0.62 


(39) 


One  has  to  consider  once  more  the  contribution  of  the  viscoelastic  character  of 
the  coating  tn  order  to  explain  this  result.  As  was  also  experimentally  observed 
m  the  peel  tests  involving  several  different  thicknesses,  it  appears  that  the 
thicker  films  require  a  higher  adhesive  fracture  energy  7  and  consequently  a 
higher  strain  energy  release  rate  G  in  order  for  interfacial  crack  propagation  to 
occur  It  remains  dfficult  to  estimate  the  absolute  value  of  the  temperature  at 
which  debonding  will  proceed  on  the  basis  of  the  elastic  (rubbery)  properties  of 
Solithane  113.  One  should  note  that  no  debond  could  be  observed  experimen¬ 
tally  unless  at  temperatures  well  below  the  glass  transition  range  of  the  poly¬ 
mer.  which  induces  a  pronounced  change  in  material  properties.  In  addition  the 
observed  debond  does  not  occur  in  the  crack  speed  range  of  the  intrinsic  frac- 
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Figure  17.  Debond  p 


eight  times  the  thickness  of  the  adhesive.  Some  specimens  were  allowed  to  cool 
down  as  the  chamber  lowers  its  temperature  ;  others  were  quench  cooled  once 
the  test  environment  has  reached  its  desired  temperature.  This  different  cooling 
history  is  reflected  in  the  initial  debond  rates,  as  shown  in  Figure  17  The  latter 
represents  the  crack  progression  as  a  function  of  time.  The  specimens  cooled  at 
a  lower  rate  start  debonding  at  a  slower  pace  than  the  quenched  specimens. 
Eventually  the  crack  propagation  arrests,  as  the  built-in  stresses  are  continu¬ 
ously  relieved  due  to  relaxation  in  the  viscoelastic  material  for  both  the  slow 
and  fast  cooling  histories.  In  addition  the  viscosity  of  Solithane  113  increases 
significantly  below  its  glass  transition  temperature,  previously  determined  to  be 
±  25  °C  5  ;  this  certainly  influences  the  value  of  the  fracture  energy  7,  a 
phenomenon  illustrated  by  the  experimental  data  of  the  previous  chapter  where 
the  influence  of  the  temperature  on  7  was  investigated.  It  must  be  pointed  out 
that  the  debond  does  not  progress  in  a  perfect  axisymmetric  fashion. 

Some  of  the  immediate,  practical  results  of  the  tests  are  the  temperatures  at 
which  debonding  could  first  be  noticed  after  a  waiting  period  of  about  12  hrs  for 
the  slowly  cooled  specimens  and  about  half  an  hour  for  the  quench  cooled  sam¬ 
ples.  For  the  larger  thickness  this  occurred  around  -42  °C  ,  for  the  thinner  film 
T=-51  °C  to  see  continuing  crack  propagation,  which  seems  consistent  with  the 
smaller  amount  of  strain  energy  stored  in  the  thinner  coating.  For  the  quench 
cooled  specimens  debond  started  to  be  noticeable  a  few  degrees  below  the 
corresponding  temperature  for  the  slowly  cooled  samples  We  now  proceed  to 
calculate  a  lower  bound  estimate  for  the  temperature  at  which  debonding 
should  start,  using  the  (elastic)  material  properties  of  Solithane  1  13  in  its  rub¬ 
bery  domain  It  should  be  recalled  that  the  strain  energy  release  rate  G  associ¬ 
ated  with  the  propagation  of  an  edge  delamination  is  given  by  equation  (26)  ;  if 

:>■  1'ne  therm ai  expansion  characteristics  of  Soiithane  1  13  were  determined  :.n  the  dilatometer, 

_ised  _n  me  material  cnaracterizauon  of  PVAe  m  Chapter  2  of  Part  1. 


Figure  16.  Specimen,  for  debond  testing. 
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3.4  Debond  Testing 

In  the  final  section  of  this  chapter  the  experimented  results  of  debond  tests 
are  presented.  The  tests  are  performed  in  order  to  make  a  comparison  with  the 
fadure  analysis  outlined  previously  in  Section  3.1  of  this  Chapter. 

3.4  1  Experimental  Method.  In  order  to  relate  the  fracture  energy  from  the  peel 
tests  to  these  measurements  Solithane  113  is  used  as  the  model  adhesive  To 
monitor  the  debond  progression  glass  is  also  again  selected  as  the  adherend 
Because  of  the  wide  test  temperature  range  involved  it  appears  necessary  to  use 
Pyrex  glass  which  can  sustain  large  temperature  differentials.  Before  the  poly¬ 
mer  is  cast  on  to  the  circular  glass  substrate,  the  latter  is  coated  with  a  thin 
(0.075  mm  or  0.003”)  self-adhesive  Teflon  film,  except  for  a  centered  circular 
portion  with  a  diameter  of  25.4  mm  (1.0").  The  Solithane  113  polymer  is  then 
cast  on  the  glass  slab  following  a  procedure  similar  to  one  used  to  make  the  peel 
test  specimens.  The  cure  at  room  temperature  eliminates  the  residual  thermal 
stresses,  but  cannot  avoid  the  cure  shrinkage  stresses  whose  magnitude  has 
been  determined  previously  [  10]. 

The  specimens,  an  example  of  which  is  shown  in  Figure  16,  were  tested  in  a 
Missimers  temperature  test  chamber  whose  temperature  can  be  controlled  to 
within  t  1  °C.  The  temperature  is  recorded  with  the  help  of  a  Chromel-Alumel 
thermocouple. 

3.4.3  Experimental  Results 

1.  The  first  series  of  tests  concentrated  on  cooling  specimens  from  room  tem¬ 
perature  to  temperatures  well  below  ambient  conditions  The  thermal  loading 
history  is  motivated  by  the  analysis  outlined  before  Two  different  thicknesses 
are  involved  in  the  measurements,  namely  0.8  mm  (1/32")  and  16  mm 
(1/16").  The  thickness  is  chosen  such  that  the  initial  debond  is  equal  to  at  least 
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the  adhesive  energy  y  from  its  lower  limit  T  reflects  the  rheological  or  dissipa¬ 
tive  properties  of  the  polymer.  This  observation  underlines  the  distinction 
between  a  rate  independent,  intrinsic  energy  of  adhesion  F  and  its  rate  depen¬ 
dent  counterpart  y  as  stated  by  Muller  and  Knauss  [l]. 

For  very  slow  crack  propagation  one  deals  essentially  with  the  intrinsic 
energy  of  adhesion  F  ;  in  this  case  the  bulk  of  the  polymeric  coating  does  not 
experience  pronounced  viscoelastic  effects.  It  seems  prudent  however  to  allow 
for  viscous  dissipation  which  occurs  in  any  case  near  the  crack  tip,  by  increas¬ 
ing  somewhat  the  value  of  the  required  adhesive  energy  y.  At  higher  crack 
speeds  the  rate-sensitivity  of  the  viscoelastic  material  becomes  much  more 
significant  and  consequently  one  needs  to  use  values  of  y  higher  than  F  m  an 
analysis  in  order  to  account  for  some  viscous  dissipation  at  non-zero  crack 
growth  rates. 

If  the  layered  structure  is  subject  to  temperature  or  moisture  variations 
which  mduce  viscoelastic  behavior  in  the  bulk  of  the  polymer  one  needs  to  con¬ 
sider  severe  changes  in  material  properties  as  well  as  the  influence  of  the 
environment  on  the  adhesive  fracture  energy.  The  former  can  be  avoided  by 
working  in  the  temperature  range  far  above  the  glass  transition  region  of  the 
polymer  coating.  The  sensitivity  of  the  fracture  energy  to  temperature  and 
moisture  makes  it  necessary  to  determine  the  adhesive  fracture  energy  y  under 
a  wide  range  of  environmental  conditions.  The  influence  of  the  temperature  can 
be  significant  as  experimentally  verified  in  Chapter  2.  By  allowing  variations  on 
the  value  on  y.  accounting  for  viscous  effects  as  well  as  environmental  condi¬ 
tions,  the  results  of  the  elastic  analyses  enable  us  to  make  engineering  esti¬ 
mates  of  the  environmental  conditions  under  which  time-dependent  adhesive 
failure  of  rubbery  materials  can  occur. 


Figure  15.  Schematic  representation  of  the  relation  between  the  adhesive 
energy  7  and  the  crack  speed  c. 


face  energy  F  ;  once  this  situation  develops,  the  interface  crack  keeps  running, 
i.e.  no  arrest  phenomenon  occurs  as  indicated  in  Figure  14. 

3.3  Viscoelastic  Effects 

The  two  delamination  configurations  were  analyzed  for  linearly  elastic  and 
isotropic  materials.  Most  encapsulating  and  adhesive  materials  however  display 
viscoelastic  behavior  when  a  wide  range  of  temperatures  is  involved. 

One  might  then  question  the  wisdom  of  using  the  results  of  theses  sample  cal¬ 
culations  presented  in  the  previous  sections,  when  time-dependent  failure  of 
rubbery  materials  is  the  subject  of  this  study  First,  it  is  worth  noting  that  we 
restricted  ourselves  to  the  fracture  behavior  of  a  crosslinked  polymer  above  its 
glass  transition  temperature  range.  Thereby  pronounced  viscous  flow  in  the  bulk 
of  the  polymer  does  not  arise  as  well  as  significant  nonlinear  viscoelastic 
response  to  deformation  which  is  more  likely  to  occur  below  the  glass  transition 
and  at  the  very  tip  of  the  interface  crack.  Some  bulk  relaxation  could  stiii 
occur,  but  its  magnitude  is  very  limited  compared  to  the  relaxation  which  a 
viscoelastic  material  experiences  when  going  through  the  glass  transition  zone 
Hence  a  corresponding  variation  of  the  elastic  properties  in  the  delamination 
sample  calculations  can  account  for  any  potential  relaxation. 

We  still  have  to  deal  with  the  viscous  dissipation  occurring  in  the  vicinity  of 
the  crack  Up.  Although  the  bulk  dissipation  can  be  calculated  more  easily,  the 
precise  contribution  of  the  viscous  processes  in  the  immediate  vicinity  of  the 
crack  tip  remains  somewhat  of  an  unknown.  In  order  to  be  able  to  incorporate 
the  importance  of  the  dissipation  on  the  unbonding  process,  one  needs  to  take  a 
closer  look  at  the  relation  between  the  adhesive  energy  7  and  the  mterfacial 
crack  speed  c  ;  such  a  graph  is  shown  schematically  in  Figure  15,  its  shape  being 
in  agreement  with  the  results  of  Chapter  2  (Figures  4  and  5).  The  deviation  of 
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(35) 


The  elastic  strain  energy  comprises  the  strain  energy  of  bending  and  the  strain 
energy  of  compression  in  the  plate  ;  the  total  elastic  energy  using  large 
deflection  theory  is  given  by  : 


U  =rrD /’ 
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n  Eh  r8 
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(  u  +  jw2  )a2  +  2i/  (u  +  ^-w2)  +  ( ^-)Z 


r  dr 


(36) 


where  a  dot  denotes  differentiation  with  respect  to  r.  The  closed  form  evalua¬ 
tion  of  equation  (36)  using  the  expressions  for  for  u  and  w  (equations  (30)  and 
(31))  is  intractable.  We  therefore  prefer  to  evaluate  the  strain  energy,  given  by 
equation  (36)  numerically,  where  £a  involves  a  thermal  strain  only.  Before 
applying  the  crack  propagation  criterion  the  energy,  stored  in  an  annulus  of 
radius  R  and  width  A  R.  which  is  about  to  debond,  has  also  to  be  accounted  for. 
This  material  volume  contains  an  amount  of  thermally  induced  energy  equal  to 

A  Eth  =  2  rr  RA  R-^—  a2  (T  -T0)2  .  (37) 

l  ~  V 

For  the  delamination  geometry  the  energy  criterion  formulated  by  equation  (15) 
can  then  be  written  as  : 


1  dU  +  1  ^  p 

2ttR  9R  2ttR  3R 


(36) 


which  is  evaluated  numerically  for  thickness  of  h  =  2.5  mm  (0.1")  and  the 
mechanical  properties  of  Solithane  113  [25].  The  results  are  presented  in  Fig¬ 
ure  14  where  the  strain  energy  release  rate  is  non-dimensionalized  by  (Eh)  If 
the  temperature  remains  constant  the  crack  starts  growing  if  the  debond  is 
sizeable  enough  for  the  strain  energy  release  rate  G  to  exceed  the  necessary  sur- 
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w  (r)  =  b 


Je(-yR)  -  J0  (yr) 

J0(7R)  "  1 


(30) 


where  b  is  the  deflection  at  the  center  of  the  plate  (r  =  0)  ;  its  amplitude  still 
needs  to  be  determined  ;  J0  denotes  the  Bessel  function  of  zero  order.  A  further 
result  of  the  perturbation  analysis  is  the  in  plane  displacement, 
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where  the  load  P  per  unit  length  of  circumference  is  : 

b2  D 

P  =  Pc  +  3.009 

R2  h2 

and  the  buckling  load  Pc  is 


(33a) 


(33b) 


The  amplitude  b  of  the  deflection  w(r)  is  now  determined  from  the  prescribed 
displacement  at  the  clamped  edges  of  the  circular  plate.  At  r  =  R  one  has  : 


u  (R)  =  e0  R  =  - 


Et 


Pc  f  3.009 


:b2 


4[J0(yR)-l 


RJ2(yR)  (33) 


from  which  b2  can  be  extracted  as  : 


1.333  ,h»2 

0  ( 1  +  0  R 

0.51 
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3  67 
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(34) 


In  equations  (33)  and  (34)  e0  carries  the  same  meaning  as  previously  defined 
through  equation  (37).  From  equation  (34)  there  follows  a  limiting  value  for  R 
below  which  no  buckling  can  occur  for  a  prescribed  specified  strain 
e0  =  -a(T  -T0)  .  a  solution  for  b  >  0  exists  only  if 
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The  second  example  involves  a  blister  type  configuration  which  may  be  the 
result  of  prior  loading  or  of  defective  manufacturing  procedures  where  portion 
of  the  coating  material  has  separated  from  its  substrate  as  indicated  in  Figure 
13.  Subsequent  heating  expands  and  buckles  the  film,  leading  to  possibly 
further  crack  propagation.  Let  a  circular  section  of  radius  R  be  detached  from 
the  underlying  adherend.  As  in  the  case  of  the  edge  delamination  T0  is  the  tem¬ 
perature  at  which  the  bond  was  effected.  In  order  to  calculate  the  strain  energy 
release  rate  G,  we  first  proceed  to  find  the  deflection  curve,  which  identifies  the 
shape  of  the  partially  debonded  coating. 

Within  the  framework  of  plate  analysis  the  problem  reduces  to  thermally 
induced  buckling  of  a  circular  plate  with  clamped  edges.  Assuming  that  the 
deflection  surface  is  a  surface  of  revolution,  and  denoting  by  <p  the  slope  of  the 
deformed  plate,  one  finds  [33]; 

V  (r)  =  A,Jt(u)  +  A2Y,(u)  (28) 

where  Jj(u)  and  Yi(u)  are  Bessel  functions  of  first  order  of  the  first  and  second 
kinds  respectively,  and 


u  =  X  r 


(29) 


where  r,  Pc  and  D  are,  respectively,  the  distance  of  any  point  from  the  center  of 
the  plate,  the  critical  (buckling)  load  per  unit  of  circumference  and  the  flexural 
rigidity  of  the  plate  Thompson  and  Hunt  [34]  present  the  results  of  a  contin¬ 
uum  perturbation  analysis  of  initial  postbuckling  behavior.  The  deflection  of  the 
circular  plate  is  given  by 
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Flgure  12.  Strain  energy  release  rate  as  a  function  of  R/h  for  the  edge  delami 
nation. 
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ture  energy  I\  However,  if  one  uses  equation  (2)  to  predict  the  temperature  at 
which  interfacial  crack  propagation  can  occur,  there  results  that  one  needs  to 
cool  the  bimaterial  assembly  (at  least)  below  -10  °C  in  order  for  debonding  to 
proceed.  This  calculation  includes  the  zero  time  glassy  response  of  the  Solithane 
113  adhesive  layer  ;  the  result  also  corresponds  more  closely  with  the  experi¬ 
mentally  observed  temperature  domains  where  adhesive  failure  occurs. 

2.  It  was  also  attempted  to  cause  debonding  by  heating  the  specimen  well  above 
ambient  conditions.  This  effort  proved  futile  as  apparently  no  high  enough 
stresses  to  cause  crack  propagation  could  build  up.  When  heating  the  bimaterial 
specimens  the  Solithane  113  expands  more  than  its  substrate  ;  this  expansion 
thus  puts  the  crack  Up  under  compression,  thereby  preventing  further  delarru- 


nation. 
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ABSTRACT 


linear  viscoelasticity  theory  is  applied  to  the  fatigue  problem  of  solids  sensi¬ 
tive  to  the  deformation  rate.  A  series  of  experiments  have  been  performed  to 
investigate  the  range  of  applicability  in  which  the  theory  is  valid  and  to  examine 
the  accuracy  of  the  theory.  The  experimental  results  are  compared  with  an 
existing  theory. 


1.  INTRODUCTION 


Compared  to  10  or  20  years  ago  there  are  many  situations  in  engineering 
today  where  polymers  are  used  for  structural  purposes.  This  use  occurs  in  the 
civil  engineering  sector;  in  water-  and  gas  piping  constructed  from  Polyvinyl¬ 
chloride;  in  agriculture  for  irrigation  and  desalination  purposes;  in  many  house¬ 
hold  goods  and  in  particular  in  the  automotive  industry.  In  the  transportation 
industry  polymers  are  widely  used  in  automobile  tires  as  well  as  in  automotive 
body  components.  In  the  aerospace  industry  polymers  are  used  to  a  large  extent 
as  sealants,  as  structural  bonding  agents,  as  well  as  bonding  together  fibers  to 
make  composite  materials.  The  composite  materials  may  be  either  of  continu¬ 
ous,  or  chopped  fiber  type.  Inasmuch  as  a  very  large  percentage  of  engineering 
application  of  any  materials  encounters  transient  or  periodic  loading,  it  is 
appropriate  that  one  give  attention  to  the  fracture  behavior  of  these  types  of 
materials  in  a  fatigue-type  of  environment. 

We  associate  with  the  terminology  'Tatigue"  the  phenomenon  that  failure  or 
fracture  of  a  component  occurs  under  a  cyclic  load  environment  at  stress  or 
load  levels  that  are  significantly  lower  than  those  that  might  have  been  sus¬ 
tained  by  the  part  if  only  a  steady  continuous  load  had  been  imposed.  In  the 
context  of  time-dependent  failure  of  polymeric  material  this  definition  may  or 
may  not  be  appropriate.  We  know,  for  example,  that  even  when  steady  loads  are 
Imposed  on  a  viscoelastic  material,  fracture  will  occur  in  a  time-delayed  manner 
due  to  the  slow  and  time-dependent  growth  of  cracks  in  such  a  material.  This 
phenomenon  has  often  been  associated  with  the  term  of  "static  fatigue".  Under 
such  circumstances  it  would  be  appropriate  by  the  classical  definition  of  fatigue 
to  add  into  that  definition  some  measure  of  a  comparative  time  length  or  life  so 
as  to  assess  whether  a  cyclic  or  repetitive  loading  leads  to  an  acceleration  over 
the  appropriate  steady-state  environment  to  which  a  part  might  be  subjected. 


By  and  large  the  study  of  fatigue  in  polymeric  materials  has  been,  to  date, 
confined  essentially  to  rigid  polymers.  In  this  connection  fatigue  in  polymers 
has  been  treated  by  virtually  the  same  methods  that  are  normally  applied  to  the 
fatigue  of  structural  metals.  Because  a  detailed  understanding  of  the  fatigue 
process  in  metals  and  also  in  the  rigid  polymers  does  not  exist  from  a  basic 
point  of  view,  one  has  generally  resorted  to  analyzing  crack  propagation  rates 
under  fatigue  type  loading  as  a  function  of  the  difference  between  maximum 
and  minimum  stress  intensity  factor.  The  resulting  test  curves  have  often  come 
to  represent  essentially  the  fatigue  material  behavior  instead  of  the  S-N  curves 
that  were  part  of  standard  fatigue  analysis  only  three  decades  ago. 

From  an  engineering  point  of  view,  this  may  be  quite  an  acceptable  approach 
for  dealing  with  the  fatigue  problem  when  the  polymer  is  definitely  in  the  rigid 
state.  However,  there  are  many  situations,  in  particular  in  aerospace  applica¬ 
tions,  where  the  material  is  exposed  to  elevated  temperatures  to  such  a  degree 
that  while  substantial  softening  is  not  achieved  the  load  duration  and  the 
fatigue  environment  may  last  long  enough  so  that  a  substantial  amount  of  dam¬ 
age  is  incurred  in  this  elevated  temperature  stage.  Inasmuch  as  at  elevated  tem¬ 
peratures  polymers  tend  to  become  significantly  viscoelastic,  the  question  arises 
as  to  how  the  fatigue  process  is  governed  in  this  kind  of  environment  when 
strong  viscoelastic  effects  are  present. 

It  appears  reasonable  to  view  the  fracture  behavior  or  failure  behavior  of 
polymeric  materials  from  the  point  of  view  of  fracture  mechanics,  that  is,  the 
failure  behavior  of  these  materials  in  the  presence  of  preexisting  cracks  In  this 
context,  it  becomes  reasonable  therefore  to  inquire  as  to  how  cracks  propagate 
in  a  material  when  strong  viscoelastic  material  behavior  is  present  and  when  the 
load  history  of  the  component  is  cyclic  in  nature.  Inasmuch  as  there  is  at  this 
time  only  an  approximate  theory  [1]  that  deals  with  crack  propagation  in 


polymers  under  cyclic  conditions  of  loading  it  is  appropriate  to  examine  this 
problem  in  some  detail  and  to  do  so  from  an  experimental  point  of  view. 

At  this  time  our  foundation  for  understanding  the  propagation  of  cracks  in 
viscoelastic  material  is  based  on  linear  viscoelastic  material  response.  For  the 
case  of  linearly  viscoelastic  behavior  this  foundation  is  outlined  essentially  in 
[2],  while  similar  developments  based  on  special  material  representation  are 
covered  in  [3]  and  [4].  These  developments  are  documented  specifically  for  pro¬ 
pagation  theory  at  constant  velocities,  where  the  underlying  assumption  is  that 
the  crack  propagation  is  applicable  instantaneously  in  situations  where  the 
stresses  at  the  tip  of  the  crack  change  continuously  with  time.  It  turns  out  that 
for  many  situations  this  underlying  assumption  is  not  a  serious  restriction  at 
all.  However,  it  appears  questionable  to  what  extent  these  relations  are  applica¬ 
ble  when  the  stresses  change  rather  rapidly  at  the  tip  of  the  moving  crack.  As 
has  been  pointed  out  in  [5],  if  the  crack-tip  stresses  change  rapidly  compared  to 
a  length  parameter  characterizing  a  thin  zone  of  fracturing  material  along  the 
crack  path,  then  these  approximations  involving  constant  rates  of  crack  propa¬ 
gation  may  no  longer  be  applicable.  The  fully  analytical  exploration  of  this  ques¬ 
tion  seems  too  cumbersome  at  this  time.  Rather,  one  would  feel  that  it  is  more 
appropriate  to  examine  this  problem  initially  experimentally.  Subsequent  to 
such  an  experimental  investigation  one  would  then  find  more  justification  or 
motivation  to  explore  particular  ranges  of  material  behavior  where  perhaps 
more  suitable  approximation  are  appropriate. 


2.  REVIEW  Of  EXISTING  LINEAR  FRACTURE  THEORY 

We  shall  summarize  here  briefly  the  developments  m  [2].  Consider  the  tip  of 
a  crack  moving  through  a  (linearly)  viscoelastic  material.  This  tip  is  character¬ 
ized  by  a  damage  or  process  zone  in  which  the  material  behaves  non-linearly.  If 
this  zone  is  sufficiently  small,  one  may  make  use  of  the  linearized  theory  of 
viscoelasticity  along  with  a  line  model  for  the  process  zone  according  to  Prandtl, 
Barenblatt  or  Dugdale.  The  process  zone  needs  to  be  incorporated  to  provide  a 
length  scale  parameter  that  couples  the  material  viscosity  to  a  crack  propaga¬ 
tion  speed.  In  Ref.  2  two  criteria  of  fracture  were  explored,  a  crack-opening  dis¬ 
placement  and  an  energy  criterion.  Both  gave  essentially  identical  results.  It 
suffices,  therefore,  to  consider  only  one,  the  energy  criterion. 

The  energy  criterion  states  that  the  rate  of  work  done  by  the  unloading  trac¬ 
tions  on  the  displacements  in  the  (line)  plastic  zone  equals  the  rate  of  fracture 
energy  absorption,  which,  in  the  simplest  form,  may  be  considered  to  be  a  con¬ 
stant  T  times  the  crack  speed  c.  Under  the  assumption  of  constant  crack  propa¬ 
gation  speed  that  criterion  leads  to  the  relation  (for  Poisson  s  ratio  v  -  1/2) 

K“®  [  —  J  =  s-E.r  (1) 

c  3 

where 

K  =  stress  intensity  factor 

a  =  length  of  (line)  cohesive  or  process  zone 

T  =  rate  independent  fracture  energy 

c  =  rate  of  crack  growth 

=  f  f~(r  ~P°)  I  ~  p(°)  I  dr  dp  =  D  (2) 

F(r)  =  VT  -  r  +  j  [  1  2  r  In  ^j==L  -  Vi  -rr] 

D(t)  =  uniaxial  creep  compliance 
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ln  this  context  we  note  that  eqn.  (1)  holds  rigorously  for  c  =  const,  and  approxi¬ 
mately  for  6  const,  provided,  (see  Ref.  2) 


iMl 

K(t) 


<  < 


c(-t) 

2a(t) 


(3) 


This  question  of  describing  crack  growth  analytically  by  an  equation  like  (1) 
was  addressed  in  a  previous  analysis  which  was  based  on  crack  growth  by  small 
but  finite  increments  [l]  but  which  did  not  limit  itself  to  the  cases  allowed  by 
[2].  The  net  result  of  that  analysis  was  that  if  viscoelastic  transients  become 
important  (i.e  if  the  equation  (3)  above  is  violated)  then  crack  propagation  per 
cycle  should  occur  more  rapidly  than  if  only  the  quasi-steady  relations  (1)  is 
involved  (6  =  const  ).  Inasmuch  as  equation  (1)  is  much  easier  to  apply  than  the 
work  in  Ref.  1  one  may  raise  the  question  whether  the  deficiency  associated  with 
violation  of  (3)  is  really  important:  After  all,  a  problem  may  exist  only  near  6=0, 
so  that  the  total  amount  of  crack  growth  accumulated  near  the  "crack  growth 
boundary”1  is  very  small  to  begin  with  and  thus  possibly  an  important  contribu¬ 
tion  to  the  toted  crack  length.  The  constraint  (3)  is  always  violated  when  a  crack 
tip  experiences  loading  that  passes  through  the  limit  for  possible  crack  propa¬ 
gation:  Below  and  at  this  limit  6  =  0  but  K(t)  *  0;  the  same  is  true  for  conditions 
above  this  limit  where  both  K  and  6  exceed  zero,  but  still  close  enough  to  the 
propagation  limit  so  that  the  order  relation  is  violated.  The  question  arises  thus 
to  what  extent  equation  (l)  can  represent  crack  growth  under  conditions  when 
the  propagation  boundary  is  passed,  particularly  passed  repeatedly  as  in  a 
cyclic  history.  It  is  possible  that  equation  (l)  would  be  more  powerful  than  the 
analytical  estimates  would  indicate.  It  appears,  however,  reasonable  to  expect 


1 .  Crack  growth  is  not  possible  if  the  stress  intensity  is  below  a  certain  value  say  K*;  above  K*  crack 
growth  occurs.  We  call  K*  the  crack  growth  boundary. 
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that  the  range  of  validity  of  (l)  is  somehow  Limited  by  the  inequality  (2)  and  that 
the  argument  of  "small  crack  growth  contribution"  means  that  the  propagation 
boundary  may  only  extend  the  range  of  validity  for  (1)  somewhat. 

In  view  of  this  situation  it  is  appropriate  to  examine  the  crack  growth 
behavior  under  transient  rather  than  steady  loading  conditions  experimentally. 
The  results  of  such  an  investigation  would  then  provide  further  input  to  analyti¬ 
cal  formulations  that  elucidate  the  behavior  of  crack  growth  under  rather  arbi¬ 
trary  transient  loading  histories.  Because  of  its  technical  importance  as  well  as 
the  ease  with  which  small  crack  growth  can  be  accumulated  over  many  cycles 
into  a  readily  measurable  quantity,  cyclic  deformation  histories  seem  to  be  well 
suited  to  this  study. 

For  later  reference  it  should  be  pointed  out  that  evaluation  of  crack  growth 
under  transient  load  histories,  if  governed  by  equation  (1),  was  carried  out  in 
considerable  detail  in  [6];  this  study  will  serve  as  a  guide  and  reference  against 
which  the  experimental  results  can  be  compared. 

We  proceed  next  to  a  description  of  the  experimental  set-up  and  analysis  of 
the  experimental  technique  (caustics).  Following  these  developments  we  turn  to 
the  recording  and  then  analysis  of  the  experimental  results. 
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It  should  be  mentioned  that  we  shall  be  interested  in  a  moving  crack,  while 
the  caustic  relation  (15)  corresponds  to  a  stationary  crack.  However,  the  crack 
will  be  moving  so  slowly  that  dynamic  effects  are  irrelevant  so  that  the  present 
analysis  is  totally  adequate  to  deal  with  the  later  test  situations. 

3.3  Experimental  Procedure 

We  shall  attempt  to  study  crack  growth  as  a  function  of  the  maximum,  stress 
intensity  during  any  cycle.  Thus  it  will  be  necessary  to  record  the  maximal  caus¬ 
tic  diameter  D  during  any  cycle.  This  is  accomplished  in  the  set-up  shown  in  Fig¬ 
ure  3.7 

A  Spectra  Physics  Model  120  Laser  shines  light  through  a  spatial  filter  and 
beam  expander  onto  the  crack  tip  of  a  specimen  supported  in  a  (servo- 
hydraulic)  test  machine  The  image  plane  is  formed  by  a  translucent  piece  of 
paper  with  mm  rulings  .  This  arrangement  allows  virtually  continual  monitoring 
of  the  caustic  size  by  means  of  a  video  camera  and  recorder.  Data  is  reduced 
later  from  the  video  tapes. 

3.4  Caustic  Data  Analysis 

In  order  to  avoid  buckling  of  the  specimen,  only  tension-tension  cycling  was 
performed  such  that  there  was  always  tension  on  the  specimen.  The  lowest  value 
of  the  caustic  was  used  to  locate  the  tip  of  the  crack  relative  to  the  boundary.  It 
turns  out  that  the  location  of  the  crack  tip  is  determined  relative  to  the  caustic 
by  the  discussions  shown  in  Figure  3.B.  The  computation  of  these  dimensions  are 
shown  in  Appendix  B  The  diameter  of  the  caustic  under  maximum  load  is 
measured  (off  the  video-play  back)  as  a  function  of  crack  length  in  the  same 
cycle.  It  will  be  noted  from  equation  (15)  that  the  stress  intensity  factor  can  be 
determined  if  all  the  parameters  characterizing  the  material  as  well  as  the 
experimental  set-up  are  known  It  turned  out  that  some  of  this  informations  was 


Loading  Cycle  Log^(HZ) 
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VISCO- ELASTIC  CAUSTIC  GROWTH  BOUNDARY 


Figure  3.0  Curve  A  corresponds  to  10  1  minutes  in  Figure  3.5 
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K,(t) 


X(t)  = 


gvg?  §• 

5. 

3f 2  h  Zo 


(14) 


Thus  by  measuring  the  "caustic  diameters"  D(t)  as  a  function  of  time,  the 
corresponding  stress  intensity  factor  K  can  be  calculated  from  (14), 

The  determination  of  K  requires  the  solution  of  an  (convolution)  integral 
equation;  this  is  not  a  trivial  matter,  in  general.  For  further  work  here  we 
therefore  make  use  of  results  in  [7]  which  establishes  an  excellent  approxi¬ 
mation  for  the  test  material  and  load  histories  used  m  this  investigation.  In 
[7]  it  was  established  that  essentially  elastic  conditions  prevail  for  times 
longer  than  10-1  min  (cf.  Figure  3.5)  at  a  temperature  of  0°C.  It  was  also 
shown  that  the  time-temperature  equivalence  principle  holds.  With  this 
information  it  is  possible  to  determine  those  conditions  on  temperature  and 
test  frequency  for  which  elastic  rather  than  fully  viscoelastic  conditions  pre¬ 
vail.  With  this  information  one  is  able  to  construct  a  curve  of  frequency  vs. 
temperature  that  "separates"  the  elastic  from  the  viscoelastic  behavior1 

This  relation  is  shown  in  Figure  3.6  along  with  the  test  conditions  (as 
points)  used  in  the  later  experiments.  It  is  clear  from  these  considerations 
that  the  caustic  data  gathered  in  the  subsequent  tests  can  be  used  to  com¬ 
pute  the  stress  intensity  factor  from  the  elastic  solution  represented  by 
equation  (15)  where  the  constant  A  (definition  following  equation  (3))  incor¬ 
porates  the  long  time  or  rubbery  Young’s  modulus  of  the  test  material. 


Ki(t)  = 


2v2tt  D(t)z 


3f  2  hZoX 


(15) 


1-  This  curve  shown  is,  of  course,  a  sharp  demarcation  instead  of  the  "region''  in  which  viscoeiastic 
behavior  is  weak.  For  this  reason  the  curve  is  adjoined  by  a  shaded  region,  in  which 
viscoelasticity  behavior  is  minimal. 


singularity  condition  (J  =  Jacobian  of  x‘.  y’with  regard  to  r,  ?) 


J=toL*L-.d!L*L  =  Q  {10) 

dr  dtp  dtp  dr  '  ' 

If  we  define  the  curve  r(?)  in  the  object  plane  from  which  the  rays  fall  on  the 
caustic  curve  in  the  image  plane  as  r  =  rs(p),  then  this  relation  yields  that 

a.  The  radius  r,  is  independent  of  p;  i.e.,  the  curve  is  a  circle  with  center  at  the 
crack  tip;  this  curve  is  called  the  'Initial  curve.” 

b.  the  radius  r,  of  the  initial  curve  is  given  by 


— 

r,=  §-^[A(t)  .  K,(t)  ]  5 


The  image  of  that  Initial  curve  (circle),  i.e.,  the  caustic  is  given  by  equation 
(9)  together  with  (11)  as 


2  3 

X’  =  r,(cos?  +  —  cosg-yi) 


(12a) 


2  3 

Y'  =  rs(sin?  +  —  sin-p>) 


(12b) 


A  plot  of  equations  (12)  is  given  in  Figure  3.4  for  a  particular  time  when  r,  is 
some  constant.  Now  define  D=2Y'max  (cf.  Fig.  3.4)  Then 

2  3 

D  =  2  max  [  r,  (  sin?  +  —  sin  )  ] 

2  3 

=  2r,max  (  sin?  +  g-  sin  ^-?  ) 

=  r,f  (13) 

where 

2  3 

f  =  2Max(sin?  +  —  sin^-?) 

J  C 


Use  of  (13)  together  with  (11)  renders  the  convolution  relation 
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As  =  h-X(<7,  +  jz)  .  (4) 

For  viscoelastic  and  optically  isotropic  materials  the  corresponding  relation  is  in 
the  form  of  a  convolution,  namely 

As  =  h  X  *  a  (5) 

where  the  star  (•)  notation  signifies  the  Stieltje's  Integral.  Thus,  the  optical 
path  change  within  the  material  is  a  function  of  the  stress  and  deformation  his¬ 
tory.  Recall  that  in  linear  elasticity  the  stresses  near  the  crack  tip  (in  the  object 
plane)  are 

_  _  Ki(l)  ,,  3tp  \  .  s 


|-(1  -  sin  |- sin  ^-)  +  ■ 

(6a) 

f-(l  4  3in f” sm  ^2~ )  + 

(6b) 

sin|- cos^- .  + 

(6c) 

It  follows  that  the  stress  invariant  a,  +  a2  -  <7,  +  o1  =  a  is 

K.Ct)  , 

°-^iZQOSz 

Therefore  equation  (1)  becomes,  with  the  help  of  Equations  (5)  and  (6) 

f  =  r  +  dZo(X(t)  •  K,(t))  7  ( ~\J~~  cos 


which,  in  Cartesian  components  is  equivalent  to 


x'  =  rcosvs  -  dz0[  X(t)  •  Kt(t)  ]  __  r  2  cos^v 

v2tt  * 


y*  =  rsimp  -dz0[  X(t)  *  K((t)  ]  — |=-r  2  sin|y?  (9b) 

We  next  observe  that  m  the  image  plane,  the  caustic  is  the  envelope  of  all  singu¬ 
lar  points  for  equations  (9a)  and  (9b).  Thus  the  caustic  can  be  expressed  by  the 
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be  deflected  to  intersect  the  image  plane  at  the  point  P"  By  the  geometry  in  Fig¬ 
ure  3.3,  we  have  then 

?  =  r  +  w(r,ip)  (1) 

where  w(r,p>)  is  a  deviation  vector.  Applying  the  theory  of  the  eikonal  one  has 

w(rt5p)  s  z0V  (As(r,y>))  (2) 

where  As(r,ip)  is  the  total  change  in  light  path  orientation  as  the  ray  passes 
through  the  material.  By  the  Maxwell-Neumann  stress  optic  law  and  Hooke  s  law 
one  has  for  As 

As  =  h  A  [  (crt  +  ct2)  ±  (cr,  —  <Ja)  ■  p  ]  (3) 

where  the  following  definitions  hold 

_  A,  +  Ag  y 

^  2  (n  ^  E 

p=  - ^2 - 

A,  +Ag-2(n-l)|- 
h  =  thickness  of  the  plate 

At,Ag  =  stress  optical  constants 
ff  1,02  =  principal  stresses 
E  =  modulus  of  elasticity 
v  =  Poisson's  ratio 
n  =  index  of  refraction 

If  At  Ag  the  materia]  is  called  optically  amsotropic,  or  birefrmgent.  The 
consequence  is  that  p  *  0  and  therefore  As  has  two  values.  Therefore  there  will 
be  two  vectors  w,  i.e.  there  will  be  two  caustics.  It  turns  out,  practically  speak¬ 
ing.  that  Solithane  113  (50/50)  is  optically  isotropic,  i.e.  A!  =  Ag,  and  p  =  0,  so 
that  only  one  caustic  is  observed.  For  optically  isotropc  materials 


Figure  3.3  Optical  paths  for  development  of  caustic. 


figure  3.2  Photo  of  a 
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crac] 


a  EXPERIMENTAL  PROCEDURES 

For  measuring  the  stress  intensity  factor  K(t)  and  the  crack  length  c(t)  in  the 
same  time  frame,  the  method  of  caustics  is  adopted.  This  method,  also  called 
"shadow  spot  method",  or  ‘the  method  of  shadow  pattern",  was  recognized  by 
Schardin  and  analyzed  by  Mannog  in  1959.  Through  the  works  of  Theocans, 
Kaithoff  and  others  (Knauss,  Rosakis,  Kim,  Ravi-Chandar)  this  method  is  not 
only  useful  for  elastic  materials,  but  also  for  viscoelastic  ones  [7], 

3.1  Physical  Description  of  the  Method  at  Caustics 

When  a  beam  of  parallel  light  is  incident  normally  on  a  stress  free  plate  of 
transparent  material,  the  beam  of  light  just  passes  through.  However,  if  there  is 
a  crack  in  the  plate,  say  under  Mode-1  loading,  the  light  rays  will  be  deflected, 
particularly  in  the  crack  tip  region.  Due  to  the  stress  field  in  the  material  the 
index  of  refraction  of  the  material  will  change.  Also,  around  the  crack  tip  the 
shape  of  the  plate  surface  is  altered.  Both  effects  will  cause  the  transmitted 
ligbf  to  deviate  from  the  original  straight  path.  By  the  Maxwell-Neumann  elasto- 
optical  law  the  refraction  index  is  decreased;  this  fact  plus  that  of  the  non- 
homogeneous  surface  deformation  makes  the  crack  tip  act  like  a  divergent 
lense  (cf.  Figure  3.1).  The  light  around  the  crack  tip  is  deflected  outward  so  that 
on  an  image  plane  the  area  around  the  crack  tip  is  devoid  of  light  and  a  dark 
spot  is  formed.  A  photograph  of  a  test  situation  is  shown  in  Figure  3.2. 

3.2  Mathematical  Description 

Consider  the  optical  arrangement  in  Figure  3.3  of  a  cracked  plate 
represented  by  the  object  plane  and  a  viewing  plane  or  image  plane  IP  is  located 
at  a  distance  from  and  parallel  to  the  object  plane.  Let  a  light  ray  be  incident 
at  point  P  on  the  object  plane  with  coordinates  P  =  P(r,  <p).  If  the  light  were  not 
deflected,  it  would  intersect  the  image  plane  at  the  point  P’.  Actually  ibe  ray  will 


not  available,  in  particular  the  stress-optic  coefficient  for  Solithane  113.  In 
terms  of  equation  (15)  this  means  that  the  parameter  X  is  not  known.  The 
determination  of  the  stress-optic  coefficient  is  not  a  trivial  matter.  Therefore  it 
was  deemed  more  time-efficient  to  ascertain  this  parameter  X  in  a  separate  cali¬ 
bration  process,  rather  than  the  stress  optic  coefficient. 

S.S  Caustic  Calibration 

To  this  end  a  typical  specimen  (cf.  Figure  3.9)  is  mounted  in  an  Instron  test 
machine  (screw-type  loading  device)  with  a  C-type  load-cell  in  the  0-20  lb  range. 
Under  a  given  load  the  caustic  is  recorded  and  the  corresponding  crack  length 
noted.  Using  the  analytical  results  for  this  geometry  (infinitesimal  elasticity 
solution  [8]),  one  notes,  as  for  example  from  equation  (15),  that  the  5/2  power 

of  the  caustic  diameter  D  is  linearly  proportional  to  the  stress  intensity  factor. 

5. 

Thus,  if  one  plots  Dz  against  the  theoretically  determined  stress  intensity  factor 
one  should  obtain  a  straight  line.  This  stress  intensity  factor  would  be  based  on 
the  conditions  prevailing  in  the  experiment  (load  and  crack  length  as  well  as  the 
rest  of  the  geometry).  As  can  be  seen  from  Figure  3.10,  the  resulting  experimen¬ 
tally  determined  points  fall  very  neatly  along  a  straight  line.  A  best  straight  line 
fic  to  this  data,  is  then  accomplished  as  shown;  this  line  determines  the  value  of 
XT5*  in  equation  (15)  and,  by  deduction,  the  long-time  value  of  the  stress-optic 
coefficient.  It  should  be  remarked  that  the  experimental  points  in  Figure  3.10 
are  obtained  by  varying  the  applied  stress  as  well  as  by  varying  the  crack  length 
for  computing  the  'theoretical"  value  of  the  stress  intensity  factor. 

For  purposes  of  completeness  let  us  record  here  the  pertinent  equation  for 
this  computation.  From  [8]  we  have  for  the  stress  intensity  factor  (Mode-I) 


where  a  is  the  average  grip  stress  and  (H,b)  are  as  defined  in  the  inset  of  Fig¬ 
ure  3.11.  The  function  iO  /}C  =  1.1215  •  F(  ^-)  is  given  in  that  figure  for  two  end 

b 

conditions  that  do  or  do  not  allow  for  rotation  of  the  grips.  In  the  present  exper¬ 
iments  no  rotation  was  allowed;  hence  only  the  appropriate  interpolated  curve 
for  h/b  =  0.67  shown  in  Figure  3.11  was  used. 

Before  proceeding  to  evaluate  the  fatigue  data  it  is  worth  while  to  mention 
why  the  photoelastic  method  was  not  used.  It  turns  out  that  with  a  specimen 
thickness  of  2.5  mm  (0.1  inch)  the  photoelastic  method  was  no  more  sensitive  or 
accurate  than  the  caustic  method.  In  this  situation  the  caustic  method  is  easier 
to  use  because  only  a  single  parameter  must  be  measured  (D)  while  in  the  pho¬ 
toelastic  method  field  data  needs  to  be  analyzed  for  a  best  fit  of  the  crack-tip 
stress  field  to  existing  photo-elastic  figures  [9). 

3.6  Data  Scatter 

A  brief  note  is  in  order  regarding  data  scatter.  It  will  be  noted  that  there  is 
very  little  data  scatter  present  in  all  the  subsequent  measurements.  This  fact  is 
a  direct  consequence  of  the  care  with  which  the  maximum  stress  intensity  was 
determined  in  each  cycle  as  well  as  the  how  the  crack  length  was  determined. 
The  latter  was  measured  under  conditions  when  the  crack  was  almost  closed 
(lowest  stress  in  any  cycle).  These  simple  rules  produced  very  consistent  data 
rubject  only  to  variations  in  material  properties  along  the  crack  path,  which 
variations  could  give  rise  to  data  scatter  measurements  of  the  crack  propaga¬ 
tion  rate.  That  this  did  not  happen  to  any  sizable  degree  testifies  only  to  the 
uniformity  of  test  material. 
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4.  RESULTS 

In  order  to  explore  the  applicability  of  the  crack  propagation  equation  (1)  it 
is  necessary  to  perform  tests  over  a  matrix  of  frequencies  and  temperature  If 
the  material  behaves  in  a  thermorheologically  simple  manner,  it  would  be 
sufficient  to  perform  these  tests  at  only  one  frequency  and  several  tempera¬ 
tures.  or  at  one  temperature  and  several  frequencies.  Although  in  previous, 
steady-state  test  the  crack  propagation  behavior  followed  essentially  ther¬ 
morheologically  simple  behavior  we  do  not  know  for  certain  what  the  effect  of 
cycle  frequency  is  on  generating  heat  at  the  crack  tip  and  thus  influence  the 
crack  growth  rate  in  this  secondary  manner.  For  this  reason  two  sets  of  data 
were  taken:  in  one  -  let  us  call  it  Set  I  -  the  temperature  varied  but  the  fre¬ 
quency  was  always  1  Hz  while  in  the  second  set  -  for  discussion  purposes 
denoted  by  Set  n  -  the  frequency  varied  and  the  temperature  remained  constant 
at  20°C.  These  test  points  are  shown  graphically  in  the  plot  of  the  viscoelastical 
caustic  boundary  (cf.  Figure  3.6). 

For  any  of  these  tests  one  records  the  length  of  the  crack  as  a  function  of  the 
number  of  cycles  as  well  as  the  corresponding  instantaneous  stress  intensity 
factor.  For  Set  I  (varied  temperature  at  1  Hz)  the  data  are  recorded  in  Figure 
4.1  and  4.2,  while  for  Set  n  (varied  frequency  at  20°)  they  are  shown  in  Figures 
4.3  and  4.4  Note  the  relatively  small  amount  of  data  scatter.  The  crack  length 
data  in  Figures  4.1  and  4.3  are  then  differentiated  numerically  by  the  five-point 
Lagrange  method  [  10,11].  Differentiation  is  accomplished  with  respect  to  the 
number  of  cycles  (as  a  continuous  variable),  though  Set  I  leads  directly  to  the 
time  derivative  because  the  frequency  is  one  Hz;  the  results  of  this 
differentiation  are  given  in  Figures  4.5  and  4.6. 

Upon  cross-plotting  the  data  from  Figures  4.2  and  4.5  one  eliminates  the 


25°C 


Figure  4.2  The  log  values  of  stress  intensity  factors  Ki  recorded  as  a  function  of 
the  number  of  cycles  in  different  temperatures 
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Figure  4.4  The  log  values  of  stress  intensity  factor  Kj  recorded  as  a  function  of 
the  number  of  cycles  under  different  loading  frequencies. 
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Figure  4. b  The  crack  propagation  speed  due  to  change  of  the  temperature. 
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cycle  number  and  derives  a  plot  of  crack  growth  rate  as  a  function  of  the  maxi¬ 
mal  stress  intensity  factor  in  each  cycle  (K,^  /Knun  ranges  between  3.0  and  3.6 
during  these  tests);  this  is  shown  for  Set  1  in  Figure  4.7.  Similarly,  if  one  cross 
plots  Figures  4.4  and  4.6  one  arrives  at  the  data  in  Figure  4.8  for  Set  II.  Thus  Fig¬ 
ures  4.7  and  4.8  constitute  the  starting  data  for  examining  the  possible  effect  of 
viscoelasticity  on  crack  propagation  under  transient  loading. 

5.  ANALYSIS  OF  EXPERIMENTAL  RESULTS 

We  notice  first,  with  respect  to  Figure  4.7,  that  the  overall  behavoir  of  crack 
growth  as  a  function  of  temperature  is  as  expected  from  the  time-temperature 
superposition  behavior.  However  two  deviations  are  rather  obvious:  If  the  data 
behaved  strictly  according  to  the  time-temperaure  superposition  principle,  the 
curves  in  Figure  4  7  should  all  have  the  same  shape.  That  is  clearly  not  the  case 
at  low  (maximal)  stress  intensity  and  progressively  less  true  as  the  temperature 
Increases 

In  this  connection  it  must  be  remembered  that  stress  cycling  occurred 
between  some  minimum  and  maximum  value  in  each  cycle.  In  these  tests  the 
minimum  value  (2.5  psi)  was  large  enough  so  that  during  most  of  the  cycle  crack 
propagation  became  possible  at  the  higher  temperature.  This  "anomaly"  of  the 
data  in  the  low-K.  high-T  region  is  thus  due  to  an  effectively  superposed  steady 
stress  intensity  factor  which  makes  a  significant  contribution  only  when  Km**  is 
small  (Kmax  /Kmin  -*  3.0)  and  when  the  temperature  is  such  as  to  allow 
significant  crack  growth  under  this  low  stress  condition. 

The  second  "anomaly  becomes  apparent  when  one  time-temperature  shifts 
the  curves  in  Figure  4.7  or  cross  plots  them  at  difference  values  of  the  stress 
intensity  factor  against  temperature.  This  plotting  is  done  in  Figures  4  9  and 
410.  In  accordance  with  the  time-temperature  behavior  indicated  by  the  dotted 


log  K  (psi 


Stress  Intensity  Factor 


Figure  4.7  Crack  propagation  rate  due  to  change  cf  temperature. 
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Figure  4.8  Crack  propagation  rate  due  to  change  of  loading  frequencies. 
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curve1  in  Figure  4.9  the  plots  in  Figure  3.10  should  exhibit  a  monotonic  curva¬ 
ture  behavior  as  a  function  of  temperature.  The  fact  that  this  does  not  occur 
but  significant  deviation  occurs  consistently  at  the  lower  temperatures  indicates 
a  systematic  variation  of  the  fatigue  crack  growth  from  standard  time- 
temperature  superposition  behavior.  In  fact  the  deviation  is  such  (the  dashed 
curves  in  Figure  4.9  correspond  to  the  dotted  one  in  Figure  4.10)  that  crack  pro¬ 
pagation  is  higher  at  these  low  temperatures  than  normal  time-temperature 
superposition  would  predict.  We  are  aware  of  the  normal  data  scatter  in  meas¬ 
urements  of  this  type,  but  believe  that  the  behavior  in  Figure  4.10  is  systematic 
at  all  stress  intensities. 

If  one  neglects  these  finer  points,  one  finds  that,  grossly  speaking  the  time- 
temperature  behavior  is  well  obeyed  in  these  fatigue  tests  as  long  as  the  tem¬ 
perature  remains  above  20°C  at  a  frequency  of  1  Hz.  This,  more  rough  com¬ 
parison  is  illustrated  in  Figure  4.11  which  shows  comparison  of  the  K-c  relation 
computed  from  equation  (1)  assuming  quasisteady  behavior  (equation  (1)  valid 
instantaneously;  dashed  curve  for  0°C,  see  Ref.  8).  The  experimental  data  for  the 
present  study  (solid  curve)  agrees  rather  well  with  the  computed  data  at  20°C 
(dotted  curve). 

Let  us  now  examine  the  behavior  at  20°C  when  the  cycle  frequency  is 
changed.  It  is  clear  from  Figure  4.0  that  a  standard  plot  of  dc  /dN  against  K 
makes  the  cyclic  growth  rate  strongly  dependent  on  frequency.  It  is  natural  to 
ask,  therefore,  whether  the  characterization  on  a  "per-cycle”  basis  is  more 
appropriate  than  on  a  "per-unit-time"  basis.2  After  all,  we  have  seen  in  Section  2 
that  the  crack  growth  should,  at  least  approximately,  be  predictable  by  a  time- 
differential  equation.  Thus,  if  one  converts  the  ordinate  in  Figure  3  8  from 

1 .  This  figure  is  an  accumulation  of  data  from  different  tests  on  the  same  material  [ 2-4 ]. 

Z.  Standard  fatigue  data  a  presented  on  a  ''per  cycle"  basis. 


(iftch/kec) 
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Flcure  4.1 1  Comparison  of  computed  and  measured  crack  propagation  rate. 


dc/dN  to  dc/dt  one  finds  in  Figure  4.12  a  much  closer  agreement,  although, 
due  primarily  to  the  low  data  scatter  in  Figures  4.3  and  4.4,  one  can  distinguish 
a  systematic  separation  of  the  curves  with  increasing  frequency.  Specifically  we 
note  that  higher  frequencies  lead  to  faster  crack  propagation  rates  than  lower 
frequencies.  For  'low"  frequencies  the  curves  all  seem  to  collapse  into  a  simple 
master  curve  which  should  be  equal  to  that  computed  from  the  quasi-static  and 
monotonically-based  crack  propagation  equation  (l). 

With  reference  to  the  previous  test  series  'Set  I"  were  low  temperature  pro¬ 
duced  crack  growth  rates  that  were  higher  than  those  predicted  by  the  quasi¬ 
static  theory  the  same  is  true  here  for  higher  frequency.  Since  higher  frequency 
has  the  same  effect  on  viscoelastic  behavior  as  lower  temperature,  we  conclude 
that  under  increased  viscoelastic  material  response  in  the  nearrelaxed  range  of 
material  behavior  crack  propagation  is  accelerated  by  a  cyclic  load  history.  This 
behavior  is  commensurate  with  the  concept  of  fatigue. 

That  the  curves  in  Figure  4.13  do  not  all  collapse  onto  a  single  master  curve  - 
thus  attributing  the  separation  of  the  curves  possibly  to  data  scatter  -  is  evident 
when  one  cross  plots  the  data  in  Figure  4.13  for  constant  frequencies  as  in  Fig¬ 
ure  4. 12.  If  all  the  curves  where  essentially  statistical  variants  of  each  other  the 
dotted  curves  in  Figure  4.13  should  have  all  zero  slope.  It  is  evident  that  that 
would  hardly  be  an  appropriate  interpretation  and  that  the  systematic  and  pro¬ 
nounced  slope  is  significant,  thus  demonstrating  the  systematic  variation  of  the 
crack  growth  rate  with  frequency. 


Log  ^  Shifting  factor  of  crack  speed  in  log  scale 


Figure  4.13  Shifting  factor  of  crack  9peed  due  to  change  of  loading  frequency. 
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APPENDICES 


APPENDIX  A  ENVIRONMENTAL  CONTROL 

In  this  Appendix  we  describe  briefly  the  experimental  set-up  with  regard  to 
temperature  and  frequency  control,  This  is  appropriate  because  relatively  small 
changes  in  temperature  may  affect  the  crack  propagation  rate  measurably  and 
singificantly.  The  schematic  of  the  interacting  components  are  shown  in  Figure 
A-l. 

THE  LOADING  SYSTEM 

The  servo-hydraulic  loading  device  is  home  built  but  essentially  from  com¬ 
ponent  that  normally  makes  up  MTS  equipment  employing  an  MTS  frame  ratio 
at  22  klb  However,  the  maximum  load  applied  was  always  less  than  10  lbs  ,  so 
the  frame  can  be  treated  as  rigid  body 

A  hydraulic  power  system  with  servo  controler  monitors  and  controls  the 
loading  history  through  a  function  generator  with  a  DC  offsetting  sine  wave  of 
fixed  frequency  and  maximum  amplitude.  Thus  the  loading  system  provided 
essentially  amplitude  (sinusoidal)  displacement  control  to  the  ends  of  the  speci¬ 
men.  Because  of  the  limitations  on  the  available  cooling  system  the  whole  load¬ 
ing  system  can  run  continuously  for  up  to  16  hours  with  steady  put  The 
environmental  chamber  (Figure  A-2)  is  a  double-wall  wooden  container  with 
dimensions  that  just  (it  between  the  MTS'  frame.  The  testing  sample,  Sohthane 
113  strip,  is  located  in  the  chamber  in  a  steady  temperature  environment.  With 
the  help  of  a  fan  air  circulates  from  a  conditioning  chamber  into  the  environ¬ 
mental  chamber  (Figure  A-2)  which  by  itself  is  an  adjustable  temperature 
chamber  The  temperature  capability  of  this  facility  is  in  the  range  100°C  to 
150°C  (in  the  low  temperature  range  cooling  is  accomplished  with  the  aid  of 
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liquid  Nitrogen). 

Two  Nickel-Chromium  vs.  Copper-Nickel  thermo  couples  are  placed  on  the 
surfaces  of  the  Solithane  113  specimen  and  the  temperature,  as  the  output  vol¬ 
tages  from  the  thermo-couples,  is  recorded  continuously  on  HP  7200B,  a  strip- 
chart  recorder.  For  the  range  of  test  temperatures.  12°C  ^  T  35°C,  the  tem¬ 
perature  control  system  can  maintain  temperature  fluctuations  on  the  speci¬ 
men  surfaces  to  within  ±  0.02  mV  of  the  thermocouple  voltage,  i.e.  average  tem¬ 
perature  error  is  ±  0.3  -  0.5°C.  This  range  corresponds  to  approximately  a 
variation  in  temperature  of  t  0.4°C. 


k  =  !  _  I  +  I  !  =  3rs 


± 

3r,  ~  9 


the  position  of  crack  tip  is  at  ^-L  .  And  the  position  of  the  maximal  distance  of 

the  caustic  from  the  x-axis,  i.e.  Y’  =  -jj-  occurs  at  a  distance  x  from  the  left  of 

4 


the  caustics  (see  Figure  3.8). 


_  jra  +  0.103r, 

L  =  3r, 


=  0.479 
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